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Abstract 



I study the post-inflation oscillation and decay of light coherent scalar field condensates 
that may develop during an inflationary phase of the universe. In particular, the 
light scalars studied are a composition of the scalar particles of a supersymmetric 
theory which correspond to the flat directions of the theory's scalar potential. Some 
toy models that possess supersymmetric flat directions are presented and numerical 
solutions for the evolution of the scalar fields are obtained. Both analytic and numeric 
results suggest that such condensates, if they existed in the early universe, can decay 
through a rapid and nonperturbative process long before these condensates could 
significantly affect the thermalization of the universe. 
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Chapter 1 
Introduction 



Inflation and Supersymmetry are central topics in modern cosmology and modern 
particle physics. There is indirect evidence in support of inflation [TJ, and there is 
theoretical and phenomenological motivation for supersymmetry |3J. This thesis is 
devoted to some problems in early universe cosmology assuming both inflation and 
low energy supersymmetry are a part of nature. The discussion will be of the dynamics 
of coherent scalar fields in the early universe and their possible effects on Cosmology. 
One coherent scalar field to be reviewed is the inflaton which is the proposed field 
responsible for inflation. The other scalars, which are the focus of the work, are those 
of a supersymmetric theory such as the Minimal Supersymmetric Standard Model 
(MSSM). These fields may have been prepared in a coherent state through the process 
of inflation. Our goal is to obtain an order of magnitude estimate for the decay (or 
decoherence) time of these coherent scalar fields after the conclusion of the inflation. 

It is worthwhile to briefly consider the larger context of fundamental scalar fields in 
nature. Scalar fields are associated with many strange effects in particle physics and 
cosmology. For example, the Higgs Mechanism, which is believed to be responsible for 
Electroweak symmetry breaking, is mediated by a fundamental scalar Higgs field(s). 
Also, inflation is believed to be due to the evolution of a coherent scalar field called 
the inflaton. Similarly, the present observed acceleration of the universe [8J may be 
caused by the dynamics of a scalar field called quintessence [9]- While composite 
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scalar fields (bound states) have been observed, not one fundamental scalar field has 
ever been observed. All the observed fields of the standard model of particle physics 
are spin 1/2 fermions or vector gauge bosons. The yet to be observed Higgs boson 
is in fact the only fundamental scalar field in the standard model. One possibility is 
that fundamental scalars simply do not exist in nature. For instance, there exist well 
motivated models of particle physics which allow for Electroweak symmetry breaking 
without a fundamental scalar Higgs field [10J. However, there appears to be no a priori 
reason to disallow fundamental scalars, and on the contrary, nature may be signaling 
their existence in the above examples. 

Scalar fields also possess peculiar classical and quantum mechanical properties. 
For instance, without some protective mechanism, the effective masses of scalar fields 
are not stable to quantum corrections and they will run up to the fundamental scale 
of ones theory which is typically the Planck scale M P ~ 10 19 GeV or the Grand 
Unification Scale Mqut ~ 10 16 GeV. This sensitivity to the ultraviolet scale is one 
aspect of what is commonly known as the Hierarchy problem. The second aspect to 
the Hierarchy problem is why these scales should be so much larger than the scale of 
Electroweak symmetry breaking which is of order TeV. We do not attempt to explain 
the Hierarchy of scales here, but this Hierarchy, which amounts to thirteen or sixteen 
orders of magnitude, will in fact play a role in our models and analyses. 

A second peculiar property of scalar fields is that they may obtain a uniform 
classical Vacuum Expectation Value (VEV) over all the observable universe. The above 
three examples of the Higgs mechanism, inflation, and quintessence are all essentially 
consequences of this property that scalar fields may obtain a VEV. Other fields such 
as the gravitational field and gauge vector fields may obtain Vacuum Expectation 
Values, but for these to be uniform over our observable universe implies a preferred 
direction to space, or could be simply reinterpreted as a scalar component to these 
fields (in the case of cosmological perturbations for instance). Vector and tensor 
fields can be classical on local scales of course, as in a laser or in the gravitational 
fields of stars and planets. Only a scalar field can possess a Lorentz invariant vacuum 
expectation value which would be essentially unobservable except for its modification 
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of fundamental properties of nature. Specifically, in the Higgs mechanism, the Lorentz 
invariant scalar Higgs VEV couples to the leptons and quarks of the standard model in 
such a way to give them the effective masses that we measure in the laboratory. The 
proposed inflaton evolved in time and thus broke Lorentz invariance; however, it was 
a uniform VEV over a large patch of space which resulted in a largely uniform energy 
density and temperature in the early universe. The remnants of this uniform bath 
are observable today as the nearly perfect thermal spectrum of the cosmic microwave 
background The concept of a vacuum expectation value of a scalar field is so 
fundamental to modern particle theories, that when it arises in a model, it is simply 
referred to as "the" vacuum of the model. 

The first property mentioned above of the instability of the scalar field's mass to 
quantum corrections can be corrected by making ones theory supersymmetric which 
cancels the instability and allows for a hierarchy of scales to be phenomenologically 
viable. Note that supersymmetry does not explain why the hierarchy is there however. 
The second property mentioned is a generic property of bosonic fields that they may 
coalesce into a Bose-Einstein condensate!]] To understand how such VEVs can 
develop in nature, it is unavoidable to look to cosmology for guidance. 

The theory of inflation proposes that there was a stage during the early universe in 
which space stretched at an accelerated rate. If the inflation was exponential, the line 
element would be written ds 2 = dt 2 — e 2Ht dx 2 . This inflation ended when the agent 
which was responsible for the acceleration had spent enough of its energy into the 
expansion that the exponential expansion could no longer continue. If one neglects 
quantum effects during inflation, then all other species of particles are diluted to 
negligible amounts by the inflation, with the only remaining constituent to the energy 
density being the inflaton. The model of inflation solved a few outstanding problems 
of cosmology at the time it was conceived [12], and other useful models of inflation 
were soon developed [121 IE]. The model of chaotic inflation [14] is the one which 
will be applied here. 

1 Yet a third peculiar property of scalar fields discussed in Section [4] is that scalar fields can have 
an equation of state with negative pressure. 
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A first ingredient to our study is the fact that inflation can magnify the quantum 
fluctuations of a field into an observable classical wave with non-negligible energy 
density. It may be shown in particular that this happens to scalar fields whose mass 
is sufficiently small [15j. This effect on the quantum fluctuations of a field is actually 
a generic property of all inflation models because it is an effect of the accelerated 
expansion itself. Specifically, all that is required is a small enough event horizon during 
inflation for the fluctuations to propagate past and become frozen [16J. The event 
horizon during inflation is approximately the inverse Hubble parameter H^ 1 which can 
be as small as 10~ 30 meters in many inflation models. It is not so unbelievable to 
imagine a quantum fluctuation propagating over such a short distance. Finally, it is 
noted that this effect of inflation on the quantum fluctuations of a field is believed to be 
the source of the cosmological perturbations [17J which has left a signature in the large 
scale structure of the universe [18J and in the Cosmic Microwave Background [19J. 
While the study of the cosmological perturbations is a central theme in cosmology, it 
will not be the focus of this work. 

The second ingredient to our study is that supersymmetry and the MSSM in par- 
ticular provides many such light scalar fields. If inflation lasted a long enough time, 
then some set of these scalar fields would have obtained non-negligible VEVs by the 
conclusion of inflation through the above effect on their quantum fluctuations. Such 
scalar field VEVs may then have implications for cosmology. For instance, one benefi- 
cial feature of these scalar field VEVs is their potential to contribute to baryogenesis 
[20] . There may also be adverse effects from light scalar field VEVs in the early 
universe, and these will be discussed shortly. 

The scalar fields we are interested in are the scalar particles of the MSSM. In 
general, these fields are coupled to one another, and if during inflation any one of these 
begin to obtain a VEV, it will induce an effective mass in the other fields it couples to. 
It may even induce an effective mass to itself through quartic self-interactions of the 
MSSM. To see this, one may expand the field as <p = $ + 5(f> where $ is the VEV (to 
be treated classically) and 5(f> are the quantum fluctuations. Considering a potential 
such as V = m 2 \(f)\ 2 + A|0| 4 and expanding in 5(f), the mass term of the perturbations 
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may be extracted. Schematically, one obtains, 

m 2 |0| 2 + A|0| 4 -> (m 2 + A$ 2 )<50 2 (1.1) 

Then assuming the bare mass m is small (specifically m < H), and assuming inflation 
lasts a sufficiently long time, the VEV $ will accumulate during inflation, and the above 
effective mass of the perturbations will grow. In some patches of the universe, the 
effective mass may grow to exceed m and become comparable to H. After this instant, 
the growth of the VEV in these patches will stop. In this way, the interaction terms 
of the scalar potential in the MSSM can suppress the growth of VEVs during inflation 
to be orders of magnitude below the Planck mass. However, in the MSSM and other 
supersymmetric theories, there are special combinations of the scalar fields for which 
the renormalizable potential is exactly zero! These are known as supersymmetric flat 
directions. The flat directions are tilted slightly by supersymmetry breaking masses 
m ~ TeV, but these masses are small compared to the Hubble parameter during 
inflation which can be 10 13 GeV in many models. In this case, assuming there are no 
other higher order Planck suppressed interactions in the scalar potential, and assuming 
inflation lasts a relatively long time, the scalar field VEV can become larger than Mp. 
Because the mass is small, the energy density of these VEVs should still remain below 
that of the inflaton during inflation. However, after inflation has ended and after the 
inflaton condensate has decayed, there are scenarios in which the energy density of 
the light scalar fields can in fact overtake that of the inflaton's decay products. One 
potential adverse effect of such scenarios is that if the light scalars possess a baryonic 
or leptonic charge, as in the Affleck-Dine mechanism [20], they could yield an overly 
large baryon asymmetry [501 I2TJ- Additionally if the scalars dominate the inflaton, their 
fluctuations will be the cosmological perturbations, and then ones modeling must be 
sensitive to the associated phenomenological constraints from the CMB (see [22] and 
references there in). 

The above problem of a flat direction field coming to dominate the energy density 
of the universe is one example of a generic problem in reconciling many extensions of 
the standard model with inflationary cosmology. It is known as the moduli problem. 
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The problem is that weakly interacting light fields with large enough energy density 
may persist in the universe for an extended period of time and may upset the successful 
predictions of the big bang cosmology. For instance if the fields decay late, they may 
generate too much entropy and result in an unacceptably small baryon asymmetry 
[23 121] or they may upset the the predictions of nucleosynthesis [24] . Alternatively, if 
these fields are stable, they may provide too much dark matter to be consistent with 
observation such as in the case of the gravitino problem to be discussed shortly. 

However, there can be benefits of having flat direction VEVs. One benefit is the 
mechanism of baryogenesis mentioned earlier [20]. Another possible benefit is a lower 
reheating temperature [25], but this benefit may be difficult to realize [26, 27], and 
why this is the case will be reviewed shortly. In general, the cosmological effects of 
flat directions depend strongly on the evolution of the inflaton and flat direction fields 
during the post-inflation phase of the universe but before the thermal phase. This is 
the phase of thermalization or "reheating" as it is commonly known and we shall briefly 
discuss this. After inflation has ended, the inflaton dominates the energy density of the 
universe, and it will be mostly uniform over large patches of space, but it will oscillate 
in time and otherwise behave as a condensate of massive zero momentum particles in 
any particular patch. It will decay by some unknown mechanism into the quanta of the 
standard model or of MSSM or of some other extension of the standard model. The 
reheating temperature mentioned above is simply the temperature achieved by the 
inflaton's decay products once they have thermalized. The reheating temperature can 
also be loosely thought of as the highest temperature that the universe ever reached^ 

An upper bound for the reheating temperature is determined by the energy density 
of the inflaton at the conclusion of inflation if one assumes that decay and thermal- 
ization happen instantly at this time. The energy density at this time is related to the 
Hubble parameter through the Friedmann equation, and measurements of the CMB 
put a bound on the Hubble parameter at the end of inflation of about Hj < 10 15 GeV. 

2 The thermalization stage is called reheating because in the first models of inflation [FJ] the 
universe was initially in a thermal phase which was followed by the inflationary phase. The "reheating" 
was then a return to thermal equilibrium after the inflationary phase. 
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From this, a bound on the reheating temperature T RH < 10 16 GeV is obtained (see 
Section [5J. At the opposite extreme, the reheating temperature must be greater than a 
MeV because nucleosythesis happens at this temperature and requires a thermal bath. 
Thus, there is approximately 19 orders of magnitude of uncertainty in the reheating 
temperature. However, we note that in most models, the reheating temperature is 
significantly lower than the above upper bound to account for a finite decay time of 
the inflaton. To summarize, both Hi and Trh are two important benchmarks for a 
cosmological model as both quantities set the energy scale for the physical processes 
which may occur. For instance, one successful feature of inflation is that it can dispose 
of dangerous heavy relic particles or topological defects such as magnetic monopoles 
which are prevalent in Grand Unified theories and also generated in large quantities at 
GUT scale temperatures [12J. Such particles are disposed of by the effect of dilution. 
Inflation thus solves the monopole problem, but there are other heavy relics which can 
be generated thermally and perhaps non-thermally after inflation has ended. If the 
reheating temperature is higher than about 10 6 — 10 9 GeV then one such hypothetical 
heavy particle, the gravitino, can be produced in large enough quantities to either 
upset the predictions of nucleosynthesis, or if the gravitino is stable, provide too much 
dark matter [281 ESI [30] . Thus, there is incentive for cosmologists to construct models 
with naturally lower reheating temperatures. A simple and well motivated mechanism 
to lower the reheating temperature is if the inflaton decays through Planck suppressed 
interactions, here referred to as a gravitational decay. The reheating temperature is 
then restricted to be below approximately 10 9 GeV which is still not quite enough to 
prevent the above mentioned gravitino problem [28, 29, 30J. The reheating temper- 
ature may be further reduced if the interactions which thermalize the inflaton decay 
products are somehow suppressed. These interactions would presumably be gauge 
interactions, and a simple mechanism for suppressing them is if the gauge bosons 
which mediate the interactions have acquired large masses [20, 25J. In fact, the flat 
direction VEVs just discussed can provide this mechanism because the VEVs typically 
break many or all of the standard model gauge symmetries [HI [32]- The effective 
vacuum established will then consist of heavy gauge vector bosons among other heavy 
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states. This would indeed suppress the thermalizing interactions and it could lower 
the reheating temperature as mentioned above. However, to obtain a significant sup- 
pression, it requires the flat directions to remain coherent over a long time and in 
such cases the flat direction's energy density may come to dominate over that of the 
inflaton [26, 33J. Flat direction domination can be problematic as discussed above. 

The suppression of thermalization depends on the flat direction VEVs remaining 
coherent over a relatively long time, and this has spurred a recent discussion over the 
time scale for the decay of flat direction VEVs £3123123121133135113313811321 of 
which this thesis is a part. A perturbative estimate of the decay rate of the flat direction 
VEVs support the hypothesis of a long lived flat direction [25], but this estimate clearly 
depends on the VEVs remaining coherent. The VEVs will decay eventually, putting 
the universe back into the gauge symmetric state, but until the decay, the vacuum is 
in the broken state discussed above. However the vacuum is also evolving - so strictly 
speaking it is not a vacuum - and it has been pointed out in [26] that the evolution 
may in fact be non-adiabatic. The time scale of the evolution is long compared to the 



fast (heavy) scales such as the oscillation of the inflatonF] In this sense it is a good 



approximation to treat the background VEVs as a vacuum. However, the flat direction 
VEVs typically possess more than one complex degree of freedom, and the relative 
evolution of these light degrees of freedom can be a source of non-adiabatic evolution. 
The corresponding light quanta to the flat direction degrees of freedom are thus not 
necessarily adiabatically invariant. The argument along the above lines is proposed 
and developed in [261 EZ], and the results indicate that non-adiabatic evolution may 
contribute to flat direction decay. Much of the new material presented in this thesis 
is based on the work of [27J and also extends and clarifies some of these results. In 
the remainder of the introduction, this work is described. 

The decay of the flat directions can equivalently be described as the onset of spatial 
instability to the flat direction condensates. Like the inflaton, the flat direction VEVs 
can be treated as condensates of massive zero momentum scalar particles, but with 

3 For a uniform field such as the inflaton or the flat direction condensates, the mass of the field 
is also the natural frequency of classical oscillation of the field. Thus, a heavy mass means a fast 
oscillation. 
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masses approximately 10 10 times smaller than that of the inflaton. The decay of the 
condensates proceeds by "kicking" quanta from the zero momentum modes into higher 
momentum states through the effect of the non-adiabatically evolving background. In 
position space, this simply corresponds to the development of spatial instability of the 
condensate. 

If the condensate is perfectly uniform over space and if all fields are treated clas- 
sically, a spatial instability will not develop. However, one expects there to be small 
noise fluctuations on top of the uniform field. Even assuming the absence of classical 
noise, there is quantum noise, and the dynamics of the uniform field is not necessarily 
stable to either source of noiseQ Specifically, both classical and quantum noise 
can be amplified by the effect of the evolving background which couples to the noise 
perturbations through effective and time dependent mass terms. The instability is 
thus produced by parametric resonance. 

The formalism here used to quantify the above instability is that of quantum 
fields in external backgrounds, and it is commonly used in quantum field theory and 
cosmology [40]. It is typically applied to handle the case of a single field with a 
time dependent frequency. The formalism was extended to handle multiple scalar and 
fermionic fields in [41] and it is further developed for scalars here. As an example 
of these new developments, consider in the single field analysis, the quantity uj/uj 2 . 
Specifically, one indication of nonadiabatic evolution is when the frequency of a state 
changes in a time comparable to or faster than one period of oscillation for the state, 
or uj/uj 2 > 1. This measure of non-adiabatic evolution is here generalized to the case 
of multiple scalar fields where one must not only keep track of the time evolution of the 
eigenfrequencies, but also the time evolution of the eigenvectors. In particular, if any 
eigenvector has changed during a time comparable to or greater than the oscillation 
period of any state of the system, the evolution can be non-adiabatic. These quantities 
can thus be arranged in a matrix. The matrix expression is relatively simple and is 
presented here and in [27], but to our knowledge has not been reported elsewhere. As 

4 Alternatively, if the condensate is coupled to other fields which can have classical or quantum 
noise, spatial instability can develop in the these other fields. An example is provided in Section H 
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in the single field case, the elements of this matrix provide a convenient preliminary 
measure of non-adiabatic evolution. They are applied in this thesis for this specific 
purpose. 

An alternative method for studying the spatial instability of the flat directions, 
which contrasts to the above formalism, is the lattice simulation of the classical field 
theory [42] . Such classical simulations are possible when the condensates have large 
occupation numbers. For instance, lattice methods have been performed successfully 
to study nonperturbative decay of the inflaton condensate [43]. Lattice methods 
may in principle also be applied to study flat directions but there are no published 
examples to date. However as these simulations are very demanding with respect to 
computer memory, the types of models that can be studied are thus limitted. The 
methods used here are more amenable to analytic arguments, are less computationally 
demanding than lattice methods, and these methods also incorporate quantum effects 
by construction. Lattice methods are not applied in this work, but have not been ruled 
out for future work. 

One limitation to our formalism is that it is performed at quadratic order in the 
perturbations. Interaction terms are not included, and also the production of quanta 
is not reflected by a corresponding depletion of the evolving flat direction condensate. 
Our models strictly do not conserve energy in this respect and thus become unreliable 
when the energy density of the produced quanta becomes comparable to the energy 
density in the flat direction VEVs. However, the formalism is known to adequately 
model the initial phase of exponential particle production. For instance, a comparison 
of the two methods for a toy inflation model has been found in qualitative agreement 
[44] . Our goal is to obtain an order of magnitude estimate of the VEVs decoher- 
ence (decay) time, and by determining the initial growth exponent, this goal may be 
achieved. The effects of interactions will be discussed at the conclusion of Section [5] 

Another important aspect to studying the decay of flat directions which is devel- 
oped here is the handling of the gauge symmetry. Due to the complexity involved in 
studying flat directions with the full SU(3) x SU(2) x U(l) standard model gauge 
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group, only U(l) gauge models are here presented^] Despite this limitation, our 
analysis correctly handles the gauge fixing of the perturbations and in particular the 
vector longitudinal modes (Goldstone mode) which do not necessarily decouple from 
the other scalar perturbations. Choosing the unitary gauge is the most convenient 
choice and once this gauge is fixed, the constraints of gauge current conservation 
become transparent. We show how the constraints of current conservation lead to 
non-trivial consequences for the evolution of the flat direction VEVs. 

We apply the above gauge fixing to two simple U(l) models involving both a 
single flat direction and also multiple flat directions [26J. The multiple flat direction 
model will be shown to yield non-adiabatic evolution. One expects the perturbations 
to the flat direction VEVs to also be the states involved in the non-adiabatic evolution 
because these are the light states of the vacuum and so evolve on the same time scale 
as the vacuum. To understand why the model can lead to a parametric resonance, one 
first tabulates the degrees of freedom. The multiple flat direction model consists of 
the gauge vector and four charged complex scalars. The time-like component of the 
vector is non-dynamical so this removes one degree of freedom, and the gauge fixing 
removes a second degree of freedom. There remains ten propagating physical degrees 
of freedom. Three of these degrees of freedom will belong to the massive vector and 
one will belong to a heavy Higgs state of the same mass as the vector. There remains 
six light degrees of freedom, and these necessarily correspond to three flat direction 
degrees of freedom of the model. However not all of the flat directions necessarily 
obtain VEVs, and in the solutions we consider, the gauge current constraint mentioned 
above allows only two of the three flat directions to obtain a VEV. The remaining flat 
direction which did not obtain a VEV mix with the Higgs perturbation and this system 
of three real fields mixes strongly and results in parametric resonance. 

It will be shown that all three states including the heavy Higgs will be produced in 
the resonance. It will also be found that the production of the heavy states compared 
to that of the two light quanta is consistent with the principle of equipartition of the 

5 Note however, that in |27j the method is generalized to SU(N) models which behave as N 
copies of the U(l) model. 
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energy. The equipartition here will not be that of a thermal bath, but rather a property 
of the strong coupling between the modes, and more akin to the equipartition seen 
in some classical systems of strongly coupled oscillators [45] The equipartition 
observed in the numeric results is also obtained by direct analysis of the Heisenberg 
equations of motion from which a scaling property of these equations is obtained. In 
particular, this scaling allows us to rescale our numerical solutions, which are limited 
by machine precision to a small hierarchy of approximately nine orders of magnitude 
to the phenomenologically relevant case in which the hierarchy is up to sixteen orders 
of magnitude. 

The numeric results of this thesis will suggest that a system of multiple flat di- 
rections can decay nonperturbatively for a broad range of parameters via parametric 
resonance. Numerical estimates of the decay time will be part of these results. In 
addition to the numeric results, some formal results stand on their own. In particular 
the scaling arguments, the adiabaticity parameters, the reformulation of the Heisen- 
berg equations of motion (not mentioned above), and the gauge fixing procedure are 
all new results, and potentially useful tools in further study. 

The plan of the thesis is as follows: In Section [2] we lay out the formalism for 
identifying and cataloguing of MSSM flat directions. This includes a motivation of 
the scalar potential (the D and F terms) of the MSSM, the formalism of gauge 
invariant monomials for identifying flat directions, and the lifting of flat directions by 
the nonrenormalizable superpotential. This section concludes with a brief discussion 
of supersymmetry breaking, and the spectrum of soft sparticle masses. 

In Section [3] we discuss the formalism for describing the quantum evolution of 
a system of coupled scalar fields in a time dependent background, or the Heisenberg 
equations of motion. In particular we continue the program begun in |4T] by presenting 
some new formal tools and discussing the solution of the Heisenberg equations when 
there is a hierarchy of mass scales. We also draw the connection between the evolution 
and an equipartition of energy. 

In Section [4] a brief introduction to classical inflation is presented followed by a 
description of the mechanism by which quantum fluctuations of light scalar field such 
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as the SUSY flat directions will be converted into a VEV through inflation. We discuss 
the post-inflation evolution of the VEVs and the likelihood of obtaining large VEVs of 
order the Planck Mass. 

Section [5] contains discussion of different scenarios for reheating but is focused on 
a standard case in which the bulk of the produced entropy is from the inflaton. In 
our discussion of reheating, we emphasize nonperturbative dynamics or "preheating" 
scenarios as they are known in the literature. The calculation of the reheating temper- 
ature and the baryon asymmetry is discussed. The potential effects of flat directions 
on the reheating process is also discussed 

In Section [6] we explain the procedure for gauge fixing a generic model of U(l) 
flat directions to the unitary gauge. We then apply the procedure to two simple U(l) 
symmetric models of flat directions, and we calculate the formal expressions which 
describe the non-perturbative decay of one of the models. 

In Section [7]the multiple flat direction model is solved numerically and the results 
for the production of quanta is presented. A useful scaling relation is described which 
allows us to rescale our results from an energy scale dictated by limitations of the 
computer to the much smaller Electroweak scale masses that are phenomenologically 
relavant. 

We conclude in Section [8] with some discussion of the implications of our results, 
and in particular which features may be generically applicable to the MSSM. We also 
indicate future directions of study. 

In the Appendix, a brief discussion of the notation of the thesis is provided, as well 
as some calculations too lengthy to show in the main text. 
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Chapter 2 

Supersymmetry and Flat 
Directions of the Scalar Potential 

In this section, we present a brief introduction to supersym metric field theory in order 
to motivate the scalar potential in these theories, and the "flat directions" which 
are the minima of the scalar potential. We describe the methods for identifying flat 
directions of a supersymmetric theory which involve gauge invariant monomials [47J. 
We then discuss the terms of the superpotential which are expected to lift the flat 
directions [HJ[22], providing a natural cutoff for the VEV a flat direction may obtain 
during inflation. We conclude with a brief discussion of supersymmetry breaking and 
the spectrum of soft sparticle masses which also appear in the scalar potential and 
which characterize the classical evolution of a flat direction during and after inflation. 
Much of the background material related to supersymmetry was gathered from the 
review [4]. 



14 



2.1 The Lagrangian of a Supersymmetric Field 
Theory and its Scalar Potential 



Below we introduce supersymmetric field theory and the scalar potential in these 
theories. 

A supersymmetry transformation will transform a bosonic field into a fermionic 
one and vice-versa while preserving the form of the terms in the Lagrangian. For 
instance, the canonical kinetic terms for a complex spin-0 scalar 0, and that of a two 
component spin-1/2 fermion ijj in the Lagrangian are[] 

d^d^cp + i^a^ip (2.1) 

A supersymmetry transformation should have the effect {0, "0} — > {0, "0} such that 
the form of the above kinetic terms will have not changed after the transformation. 
A simpler symmetry which one may compare to supersymmetry, is the symmetry of 
the following Lagrangian involving two complex scalar fields 0i and 2 , 

d^fa + d^fa - m 2 (|0i| 2 + |0 2 | 2 ) (2.2) 

This Lagrangian is invariant under SU(2) transformations of the doublet (0i, fa) which 
rotates the scalar fields 0! and 2 into each-other. The transformation matrix can 
be written e Jer where e = (ei,e 2 ,e 3 ) is a vector specifying the rotation and r = 
|(ci)C2,cr3) are the group generators. The infinitesimal form of this transformation 
is, 

h ) V 02 / 

For example, the standard model Higgs field transforms in this way. A supersymmetry 
transformation is slightly more involved than an SU(2) transformation however as the 
generators of the supersymmetry transformation are anti-commuting spinors. Under an 



1 The notation here is ct m = (<7q,— (?) and er M = (cto,<?) where ctq = ( g i ) ar| d or are the 



p I ■ , ■ , 1 \ / -i \ (10 

Pauh spin matrices o\ = [ ^ Q I , er 2 = I . Q ,(T 3 = q _ 
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infinitesimal supersymmetry transformation a scalar <p and a fermion ip will transform 
via an anti-commuting two component spinor e as follows, 



6</>^<f> + ^ a , i^i + i(a"4^ (2.4) 
where a = 1,2 indexes the two components of the spinor. This infinitesimal SUSY 



transformation is to be compared with the infinitesimal SU(2) transformation (2.3) 
shown above. There is a further subtlety in the above SUSY transformation that 
for the transformation to be consistent, it is necessary that the number of bosonic 
degrees of freedom in <\> match the number of fermionic degrees of freedom in ip. At 
a cursory level there is a problem since the spinor is described by two complex values 
while the scalar is described by one complex value. Of course, there are only two 
on-shell degrees of freedom in both the scalar and spinor, but the other two spinor 
degrees of freedom can go off-shell when calculating a quantum amplitude^ and we 
will want our transformation to be applicable both on-shell and off-shell. In order to 
match both the on-shell and off-shell degrees of freedom, a complex auxiliary field F is 
introduced. This field will not be dynamical. It will be rather a Lagrange multiplier in 
the classical sense, but it will make up the difference of off-shell degrees of freedom for 
the scalar so that the fermion may consistently "rotate" into the spin-0 boson. The 
auxiliary F field is necessary for the supersymmetry transformation to be consistent. 
The combination of a complex spin-0 scalar <p, a two component spinor ip, and the 
auxiliary field F is known in the literature as a "chiral multiplet", and it transforms 
consistently into itself under supersymmetry transformations. The chiral multiplets 
have the following form in the Lagrangian, 

^chiral = d^d^i + i^d^i + F*F { + £ chiral int (2.5) 

which is simply the canonical kinetic terms shown previously, plus a quadratic part in 
the the auxiliary field F, plus interactions between multiplets encoded in £ c hiraiint to 
be specified shortly. Every standard model fermion (a quark or lepton) is expected 



2 by off-shell, we mean that the Heisenberg uncertainty relation allows the relation p^p^ = m 2 
to be violated for finite periods of time (ie for a virtual particle). 
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to belong to a "chiral multiplet" such as above, and so is expected to have a spin-0 
counterpart (a squark or slepton). Additionally, the as yet to be discovered Higgs 
sector fields are also expected to belong to chiral multiplets, and have their respective 
spin-1/2 counterparts (Note the MSSM requires two Higgs doublets versus one for 
the Standard Model). 

The interactions, £ c hiraiint shown above must preserve supersymmetry as well as 
gauge invariance, and these two constraints are enough to restrict its form. Without 
going into the details of the derivation, £ c hiraiint may be encoded into a complex 
quantity referred to as the superpotential W. The interaction Lagrangian and the 
superpotential are then defined, 

^chiral int = " \ (W^j + W^^ + C.C. (2.6) 

W = b a <$> a + m a6 $ a $ 6 + A a6c $ a $ fe $ c (2.7) 

SW 5 2 W 
w % = — W lJ 



5$i ' SQiSQj 
where the standard convention of writing W in superfield notation is used. The 
superfield $ a is a compact notation to represent both the scalar <p a as well as its 
fermionic superpartner ip a . Note also that W is not in truth a potential because it is in 
general complex valued. To summarize, W is invariant under both supersymmetry and 
the gauge symmetries, and it encodes the information required to construct £ c hiraimt 
as shown above. 

The above defined scalar fields and their superpartners must be charged under 
the standard model gauge group. Our SUSY models then require gauge vector par- 
ticles A a and each of these will possess a fermionic superpartner, A" (a gaugino) 
for supersymmetry to be preserved. For the same reasons as we required an auxiliary 
field F for the chiral multiplet, we will require a real scalar auxiliary field D for the 
gauge vector multiplet. The counting of the off-shell and on-shell degrees of freedom 
in this case are as follows. The vector and spin-1/2 particles both have two on-shell 
degrees of freedom. Off-shell, the vector has three degrees of freedom as the gauge 
symmetry always allows us to eliminate one of the degrees of freedom. We must then 
add one degree of freedom to make up the difference, and this is provided by a real 
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scalar auxiliary field D. The vector gauge multiplet thus contains a real vector A a , a 
spin-1/2 fermion A and real scalars D a which appear in the Lagrangian as follows 



-£total 

r 

•'—gauge 



-^chiral ^gauge ^additional int 

--F°F^ a - i\ ]a a»V„\ a + -D a D a 
^ p,v 2 

d.Al - d v Al + gf abc A b Al 



(2.8) 
(2.9) 



abc A b \ c 



d,X a + gf abc A b \ 



where f abc are the gauge structure constants for the specific vector gauge particle 
contained in the multiplet. The interactions of these gauge particles is partly specified 



by replacing gauge covariant derivatives for the regular derivatives in (2.5). However, 
there are additional interactions between the gauge particles and the chiral multiplets 
which are allowed by supersymmetry and gauge symmetry. These are parametrized by 

Additional int above and have the form, 



C 



additional int 



-V2g [(0*T»A a + c.c] + g{<fT a <p)D a 



(2.10) 



where g is the relevant gauge coupling. To be complete, the supersymmetry trans- 
formation for the fields after inclusion of both chiral and gauge vector multiplets is 



generalized [4] from (2.4) to 



5F 



SD a 



iip^cP^i + e a Fi 
ieW^i + V2g(T a <f))^^ a 
— (eVA a - A ta c^e) 

2^^ + Tf^ 



V2 



(eVZ> M A° - P M A ta a^e) 



The transformation is quite complicated, but this should not obscure the general idea 
which is very simple. To summarize, the Lagrangian of our supersymmetric field 
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theories is defined by the Eqs. 2.5 - 2.10 The scalar potential of the candidate theory 
is determined by solving the equations of motion for the auxiliary fields F a and D a 
(which are Lagrange Multipliers) and substituting these back into the Lagrangian. The 
result of this procedure is the following, 

V{W) = ^Zn<^0*) 2 + ^|F a (0)| 2 (2.11) 

a a 

D a = g a Y,^T a ct> (2.12) 

a,b 

dW 

= W a (2-13) 

where the first line is the scalar potential, and the latter two lines are the equations 
of motion for the auxiliary fields. Note that the scalar potential ultimately arises as 
a constraint from the demands of both supersymmetry and gauge invariance since (i) 
the auxiliary fields were required to satisfy the demand of supersymmetry and (ii) the 
interaction Lagrangians were required to satisfy both demands. Note also that the 
scalar potential is strictly greater than or equal to zero. The above expressions for the 
scalar potential, and the auxiliary fields will be referred to frequently in the subsequent 
sections. 

With the scalar potential now motivated, we return briefly to a larger scope and 
to the analogy between supersymmetry and the SU(2) symmetry. One requirement 
that both supersymmetry and the SU(2) symmetry demand is that the fields which 
are transforming into each-other must have the same mass. However, in the same way 
that the SU(2) symmetry may be broken spontaneously, the supersymmetry may also 
be broken spontaneously so that it is effectively hidden from us. One consequence of 
supersymmetry breaking is that the masses of a particle and its superpartner generally 
become different. We conclude that if the fields in nature are truly supersymmetric, 
then supersymmetry is necessarily broken, because experimenters have never observed 
a fundamental boson with the same mass and charge as the electron. If such a particle 
exists, its mass must be larger than the mass scale probed by existing colliders. The 
topic of supersymmetry breaking will be raised again when we discuss the the MSSM 
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and the parametrization of supersymmetry breaking in the the MSSM's effective La- 
grangian. The scalar potential above will then be augmented by soft SUSY breaking 
contributions. 



In the next section we identify the roots to the scalar potential (2.11), which are 
the minimum field configurations called flat directions. 



2.2 D-Flat Directions and How to Identify Them 

The Lagrangian of supersymmetric field theories just discussed in general possesses 



directions in field space for which the scalar potential (2.11) is zero. The fields which 
compose these flat directions may thus in principle be excited to large classical field 
strengths at no cost to the potential energy. As we saw in Section [4} inflation provides 
a mechanism by which these large VEVs for flat directions may develop. Of course, 
there will be soft supersymmetry breaking additions to the effective scalar potential 
as well as non-renormalizable terms which may lift the flat directions. Both types 
of terms will limit the growth of the flat direction VEVs during inflation with the 
limiting effects of the non-renormalizable terms becoming more relevant at larger 
field strengths (as these terms possess higher powers of the VEV amplitude). Also, 
there may be supergravity corrections in the effective potential which, if present could 
strongly limit the growth of a VEV [4S1 [32] ■ The presence of SUGRA corrections are 
model dependent however, and there exist well-motivated models of supergravity such 
as "no-scale" supergravity |49J in which such corrections are absent at tree-level and 
thus su pressed [5U] , 

In this section, we temporarily put aside the above mentioned subtleties of lifting 
a flat direction, since our first task is to identify the flat directions. To solve for 



all the roots of the scalar potential (2.11) is a formidable task, but one may use 



gauge symmetry to simplify the problem. Specifically, to solve for these special field 
configurations, one may use the theorem of [47J which relates solutions of D a ((p a ) = 
to the gauge invariant polynomial of the fields <p a composing the flat direction. Having 
identified these solutions, one may then impose the F-flatness constraint F a = 0, to 
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further isolate the flat direction manifold. The theorem of [47J is summarized next. 
To begin, define the set of $ a 7^ to be a particular solution of the D-flat 



constraint, D a = which, using (2.11), is written 



a 



for all a 



(2.14) 



In the following we refer to <p a as fields, although the proof only requires that they 
are complex quantities which transform under the gauge symmetries. Recall that an 
infinitesimal gauge transformation may be written 



iT a e a )<f) 



(2.15) 



where T a are the gauge generators of the group, and e a are infinitesimal parameters. 
Assuming one can construct a gauge invariant polynomial I((j> a ) of the same set of 
the (p a , and performing an infinitesimal gauge transformation on the 4> a , the function 
/ transforms as, 



7(0 o + iTV> a ) 
dl 



The gauge invariance of / then requires, 

91 T°6„ 



d<p a 



for all a 



(2.16) 
(2.17) 



(2.18) 



This statement for the gauge invariance of / becomes the statement of D-flatness 



(2.14) if one chooses a = $ a and if one can additionally show that 

dl 

C$! for all a 



w 



(2.19) 



here C is a nonzero complex factor independent of a. The authors of [47] point out 



that this additional condition (2.19) is typically obtained. 

As an example of the above correspondence between flat directions and gauge 
invariant polynomials, consider a toy model with a U(l) symmetry consisting of two 
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complex scalar fields <p a and 



potential (2.11) is, 



b b with charges +1 and —1. The D-term in the scalar 

1 



-g\\4>a\ A -\M 



2\2 



(2.20) 

and the flat direction is simply {$ a = fe a+e , $6 = fe^ 9 } where / is real and 9 
is a complex phase which can be eliminated via gauge transformation. The gauge 
invariant polynomial corresponding to this flat direction is the monomial, I2 



and the condition (2.19) may be checked, 

dh 



0(j)a 

dh 



P 
<$>a<$>b 
P 



(2.21) 
(2.22) 



where it is apparent that the common complex factor in this case is C = We note 
that the MSSM monomials studied in the following sections will satisfy the constraint 



(2.19) in a similar way as above. 

To summarize, by constructing gauge invariant polynomials I(4> a ) with the property 



(2.19), we implicitly find solutions to the D-flat constraint. Next, the F-flat constraints 
may be imposed. To be precise, a D-flat direction is lifted, when there is a nonzero 
F a = |^ in the scalar potential (2.11) which can happen when there exists a term 
in the superpotential (2.7) with one or ess fields that do not obtain a VEV in the flat 
direction. It is helpful to notice that the terms of the superpotential are themselves 
gauge invariant monomials and there are some flat directions which are lifted just by 
their corresponding monomial appearing in the superpotential. 

One comment is in order with regard to the flat direction VEVs and to the larger 
goal of studying the dynamics of these flat directions which is perhaps obvious, but we 
will mention none-the-less. The VEVs $ a can take on a continuum of values - they are 
not isolated solutions to the D-flatness constraint - and the solution space is typically 
a n-dimensional surface with n > 1. Hence, there is no obstacle to constructing space- 
time dependent solutions $ a (x M ). The price to pay of course is an increase in the 
energy of the system from the kinetic and gradient terms for $ a which would appear 
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in the Lagrangian. This suggests that a flat direction VEV with an initially uniform 
value such as would be generated during inflation can decay via spatial instability - 
one need only study the dynamics to determine if this realized in any particular model. 
This is done in later sections, but we now return to the topic of identifying the flat 
directions. 

The above procedure for determining the F and D-Flat directions is systematically 
carried out for the MSSM fields in [31J and in [22]. In the next section, we outline the 
procedure and we present a few examples. 

2.3 Identifying D-Flat Directions of the MSSM 

The correspondence between flat directions and gauge invariant polynomials can be 
exploited to catalog the possible flat directions in the MSSM by computing all the 
monomials invariant under the standard model gauge group SU (3) x SU(2)l x U(1)y- 
This has been done in [31] and [32]. The procedure of j3T] for constructing monomials 
of the MSSM fields is summarized as follows: First, the color indices of all the fields are 
contracted to form 377(3) singlets in a basis of monomials denoted B 3 . The SU(2)l 
indices of the fields in B 3 are then contracted to form all the possible S77(3) x SU(2) l 
singlets in a basis B 32 - Finally this basis B of all S77(3) x SU(2) L x £7(1) singlets are 
formed by combining the polynomials in B 32 into hypercharge zero combinations. It 
can then be shown that not all of monomials constructed in this way are independent. 
There are nonlinear relations between some monomials, and the authors of [21] are 
able to reduce the set of monomials with these relations. 

We will explain this in slightly more detail. The notation used is that of [31] in 
which Greek letters a, ft, 7 refer to SU(2)z indices, latin letters a, b, c, ... refer to color 
indices, and latin letters i,j,k,... refer to family (generation) indices. The uppercase 
latin letters I, J, K, ... are also used below and their meaning should be clear from the 
context they are written. 

The basis B 3 is constructed by contracting the color indices of the squark fields. 
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The SU(3) invariant combinations are, 

(QiQj) 
(QiQjQk) 

(qiQjQk) 



= QiQaj (2.23) 
= QiqW^abc (2.24) 
= qaiq b jqc K e abc (2.25) 



where I,J,K,... = 1...6 index the quark, and the symmetry under interchange of 
/, J,K... is not specified. For the squark SU(2) L doublets, the following two classes 
of monomials are possible, 

{QiQjQkT = Qi a Q] b QT^ bc e^ not all family indices equal (2.26) 

{QQQt M) = QTQfQl c tab C e ijk (2.27) 

where the first of the above transforms as an SU(2) L doublet and the second of the 
above transforms in the spin-3/2 representation. The SU(2)l singlets for the basis 
B 32 are constructed from the doublets via, 

(<pnpj) = tfifie^ (2.28) 

Finally, the basis B is constructed by combining the monomials of B 32 which appear 
with the hypercharges Y = {—2, —1, 0, 1, 2} into monomials of the form 

£o , Zit-i , , , (2.29) 

The following relations between monomials are also applicable to reducing the above 
bases B 3 , B 32 and B into smaller bases, 

q^qMiqu) = q a K {qiqLqM) + qHqmiqu) + q a M {qKqLqi) 
{qiqjqK){qLqMqN) = {qiq~L){qjq~M){qKqN) ± (permutations) 

{.¥i¥j)(.VkVl) = (<Pi<Pk)(<Pj<Pl) + (<Pi<Pl)(<Pk<Pj) 

(tlt J -i)(t?&) = (^-i)(£f^i) 
(6^X6^) = £1X6^1) 
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LLe 


= L a L^e Qf3 e 


QdL 


= Q"d a L l3 e af3 


QuH u 


= Q a a u a H^ 


udd 





where the first three above have been obtained from the properties of Levi-cevita 
symbol. It is instructive to look at a few simple examples of monomials (with the 
gauge index structure also shown), 

(2.30) 
(2.31) 
(2.32) 
(2.33) 

where the family indices have been suppressed. Also note these particular monomials 
will appear again in the next section when R-parity is discussed. Tables of the numerous 
other monomials can be found in [311 [32]- In the next section, we look at the lifting 
of the D-flat directions from the F-terms and from the soft supersymmetry breaking 
terms of the MSSM. 

2.4 The F-terms of the Scalar Potential and Lift- 
ing of D-Flat Directions 

The D-flat directions are lifted by three main contributions; the renormalizable F- 
terms, the non-renormalizable F-terms, and soft supersymmetry breaking terms of the 
Lagrangian [2UE2]- We will discuss the renormalizable F-terms first. The renormal- 



izable F-terms are determined from the superpotential (2.7) via (2.13). The MSSM 
superpotential has the form, 

Wmssm = uy u QH u - dy d QH d - ey e LH d + fiH u H d (2.34) 

where the y's above are 3 x 3 dimensionless yukawa matrices which operate in the 
family space and which are assumed to have magnitude of order unity. The gauge 
invariant contraction of indices is done in the same way as in the previous section, 
and in fact the above terms in the superpotential are nothing else but gauge invariant 
monomials. The fields H u and H d are the two Higgs SU(2) L doublets of the MSSM, 
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and the /i parameter has dimensions of mass and it is taken to be order of magnitude 
TeV [4]. There are some terms that could have been included into the MSSM su- 
perpotential, but are disallowed because they violate baryon number conservation or 
lepton number conservation. These terms are, 

W omitted = ^LiLjEk + X' ijk LiQjd k + ^UH U + ^X'^mA (2.35) 

The first three terms above will lead to lepton number violation and the last term 
will lead to baryon number violation. One disastrous consequence of such terms is a 
rapid decay of the proton. Rather than imposing B and L conservation separately, all 
of these dangerous terms above can be eliminated with the assumption of a discrete 
symmetry called R-parity, 

P R = (_ 1 )3(B-L)+2. (2>36) 

where s = ±| is the spin quantum number of the field involved. From the above 
definition, it can then be shown that all of the standard model particles (including the 
Higgs Bosons) have R-parity P R = +1, while all of the superpartners have R-parity 
P R = — 1. That R-parity is conserved by the MSSM means that all the terms in the 
MSSM Lagrangian should have an R-parity of P R = +1. One notable consequence 
of assuming R-parity is that the lightest supersymmetric particle is required to be 
absolutely stable, and thus a promising dark matter candidate^] 

To see how the D-flat directions may become lifted by the F-terms, consider the 



four flat directions specified previously in (2.33). The LLe, QdL, and udd monomials 
of (2.33) would clearly be lifted by the corresponding terms in (2.35). However R- 
parity prevents these terms, and thus allows these three directions to be both D-flat 
and F-flat. R-parity does not protect the QuH u direction however, so this direction is 
lifted. 

Note that for those flat directions that are not lifted by the F-constraints such as 
LLe, QdL, and udd, their flatness will be preserved at higher loop corrections in the 
renormalization group by the SUSY non-renormalization theorems. One must then 
look at non-renormalizable terms in the superpotential to determine when a particular 

3 For a more comprehensive discussion of R-parity, see [4| 
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D-flat direction is lifted. The non-renormalizable superpotential is composed of gauge 
invariant monomials of order n = 4 and higher, and the lifting is determined in the 
same manner as is outlined in Section 12.21 The dimensionful coefficients to these 
terms are assumed to be in inverse powers of the GUT scale or some other ultraviolet 
cutoff M. The potential can thus be written as a series expansion 



(2.37) 



V($) =m 2 |$| 2 + K>3(|$|) , Wl$l) = $^|An| 2 M 

71=4 



M 



where m is a soft supersymmetry breaking mass to be discussed in the following 
section, $ is the VEV of the flat direction, and n indexes the corresponding nonzero 
term in the superpotential. The lifting of all the MSSM flat directions by terms in 
the superpotential is discussed in detail in [31, 32J, and in [3T], the corresponding 
complex dimensionality of the flat direction subspace is also listed. The flat direction 
VEV obtained during inflation is partly controlled by the leading nonzero term in the 
superpotential which is specified here and in [31] by the index n. The main result of 
this section which will be applied later is, 

The majority of the MSSM flat direction degrees of freedom 

(2.38) 

remain flat at least until n = 3, 4, 5, 6 

There are only two complex degrees of freedom remaining which could remain flat until 
n = 7, 9. We note that while n is a dummy index running over integers n > 3 above, 
for the remaining chapters we will typically use n to mean the index corresponding to 



the lowest non-vanishing term in the series (2.37). We emphasize that the existence 
of these terms which can lift the flat direction are model dependent. Based on known 
symmetries these terms can be present, but this does not mean that they actually are 
present. 
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2.5 SUSY Breaking and the Spectrum of Soft 
Sparticle Masses 

The effects of SUSY breaking on a model of flat directions appear in the parameters 
of the effective scalar potential and if one further assumes supergravity, the effects 
may also appear in the Kahler potential, or equivalently the kinetic terms [221 EE]. In 
our models, we only consider effects on the scalar potential, and in particular the mass 
terms. However, we briefly discuss the supergravity case as well as the implications 
of the spontaneous supersymmetry breaking dynamics occurring during inflation or 
reheating. 

In general, global supersymmetry is broken when the scalar potential obtains a 
VEV, and this can only happen when the equations F = and D = cannot 
be simultaneously satisfied for any value of the fields. The F and D terms of the 
scalar potential then become supplemented by supersymmetry breaking corrections 
with a characteristic energy scale m susy ~ TeV which is the energy scale suggested 
by particle phenomenology and cosmology. Models in which the MSSM fields are 
directly responsible for SUSY breaking are difficult to reconcile with phenomenology. 
Instead the belief is that supersymmetry is broken by "hidden" sector fields [4J which 
communicate the breaking to the MSSM "visible" sector indirectly through a flavor 
blind interaction such as gravity (Planck-scale Mediated Supersymmetry Breaking) 
or through the standard model gauge interactions (Gauge Mediated Supersymmetry 
Breaking). The most obvious effect on the "visible" sector is the lifting of the sparticle 
mass spectrum to the TeV scale. The focus here will in fact be on the mass spectrum 
of the squarks, sleptons and scalar Higgs fields since these particles comprise the 
possible MSSM flat directions. 



A generic superpotential (2.7) can have mass parameters, and the MSSM (2.34) 
contains one such term which is the fxH u Hd term (/z is also assumed to be of the TeV 
scale). A new breaking potential £ soft is added to the model in which only one particle 
from any supersymmetric boson-fermion pair need appear^] In this potential, there 

4 The standard model partners do not appear in Cl^l ar since these masses come only from the 
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will be mass terms for the squarks, sleptons, Higgs fields and gauginos, as well as 
cubic interaction terms [4], but we restrict our attention to the scalar mass terms 
which appear with the following parametrization, 

+m 2 Hu H* u H u + m 2 Hd H*H d + (bH u H d + c.c) 

where m? are 3 x 3 Hermitian matrices in the family space. When one speaks of 
the spectrum of masses, it is the above parameters which specify this. There are 
some mass degeneracies for the particle states which are immediate consequences of 
gauge symmetry. Namely that the components of any SU(2) doublet are degenerate 
in mass and similarly for any SU(3) triplet. This fact is implied by the above notation. 
Additionally, it is expected that the overall mass scale for the above terms should be the 
TeV scale. Also there is phenomenological evidence related to CP violation, Lepton 
number conservation, and flavor changing neutral currents which further restricts the 
values of the above mass terms [4|. Without going into detail, we note a simple 
paradigm for satisfying these phenomenological constraints which is to assume each 
of the above matrices m? are proportional to the identity matrix. This leaves 6 
parameters to specify the matrices completely. This paradigm is known in the literature 
as "soft supersymmetry-breaking universality". However, we note that in stronger 
versions of universality all the mass matrices are degenerate to a single mass scale 
m . The universality paradigm is connected to the expectation discussed above that 
supersymmetry is broken by a hidden sector and communicated to the MSSM sector 
by a flavour independent interaction. The flavour independence of the interaction 
would naturally lead to degeneracies imposed by universality [4|. 

This is not the whole story however because the mass parameters are energy scale 
dependent, and they "run" with renormalization group equations. The degeneracy is 
only exact at some universal scale which is usually taken to be the GUT scale. The 
renormalization group equations for the MSSM are then used to solve for the mass 
parameters at the typical interaction energy scale Q of ones experiment. In our case, 
Higgs VEV after Electroweak symmetry breaking. 
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the energy scale of the flat direction dynamics is Q ~ TeV as the flat directions mass 
and decay products would be of this scale. The renormalization group equations result 
in corrections to the masses which are on the order \n(Q/Qo) where the input scale 
Q ~ M GUT is assumed. In this sense the corrections are "soft" compared to a theory 
without supersymmetry which typically has quadratic corrections with the energy scale 
Q, and thus enormous mass corrections [4j. The notation "soft" for our susy breaking 
potential £ so f t thus refers to this property. If the renormalization group equations are 
solved assuming some form of universality, typically one finds the first two flavors of 
any family in the squarks and sleptons remain nearly degenerate, but the third flavor 
of a family receives corrections within a factor of ten or so of this degeneracy [4]. 

To summarize, the assumptions we make on the spectrum of scalar masses in our 
later models is modest, that they are all near the TeV scale but not all degenerate. 
It is noted however that universality is not necessary in the MSSM, and the spectrum 
may in fact have less degeneracy but probably not more. 

There can be other new effective interactions to the scalar potential which appear 
due to supersymmetry breaking and which can be mediated by new heavy quanta of 
mass at or slightly below the Grand Unification Scale M GUT ~ 10 16 GeV. The class of 
terms we are interested in are those which violate the difference of Baryon and Lepton 
number B — L. Such interactions are not part of the MSSM, but they are present in 
some Grand Unified models such as SO(10) models (for example, see [51J). We do not 
attempt a discussion of Grand Unified models which are far beyond the scope of this 
thesis. However we will motivate the nature and origin of these effective interactions 
by again considering our toy model. Our toy model will in fact not contain all the 
desired elements we wish (it will not allow for CP violation), but it will be possible to 
proceed. Specifically, we will determine the effective interactions applied in our later 
numerical simulations of Section 

The toy model introduced in Section |Z2] is composed of two complex scalar fields 
0i and 02 charged oppositely under a U(l) symmetry and with masses of the scale 
m susy . In addition to the local U(l) gauge symmetry, the Lagrangian and in particular 
the D-term, is invariant under two separate global gauge transformations 0i — > e lAl 0i 
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and <p 2 — > e tA2 4> 2 - The model will thus conserve separately the numbers of 0i and 
02 quanta. One can loosely think of these as individual lepton flavor numbers for 
instance. Now consider adding to this theory a new neutral heavy vector X which 
violates the number conservation of 0x and 2 with a decay such as X — > 0102 or 
X — > 0*02; ■ The heavy particle would also mediate a process such as 0102 - ► 0102 
which violates the 0i number and 02 number but conserves the gauge charge. Since 
these X particles are very heavy, they would be created as virtual quanta to mediate 
scattering such as above, but never appear on-shell as an initial or final state quantaj^] 
The quanta would be unobservable and the above interaction would appear in all 
respects as due to a quartic interaction of the form A0 2 02 + c.c. In particular the new 
quartic interaction would be interpreted as a correction to the classical Lagrangian 
and thus a correction to the classical equation of motion and dynamics. 

Our next task is to estimate the magnitude of the coupling A based on the tree- 
level interaction. The scattering 0102 — > 0102 is modified by the presence of the 
condensate (0i) = (02) = $ and in particular the propagating heavy particle will 
acquire an effective mass m x ~ max [|$| 2 , m x ] where m x ~ M GUT is the mass in 
the absence of the flat direction VEV. It is assumed that the VEV is large, |$| > m x 
so rn 2 XcB ~ |$| 2 . Thus the effective coupling acquires a factor l/m x from the X 
propagator and a factor m 2 usy which is the typical momentum transfer at the energy 
scale we are interested in. To summarize, our effective couplings have the following 
form and magnitude, 

'si 

li> 



Til 

Xfyl + A*0f 0* 2 2 , A~^f (2.39) 



Again, this interaction is U(l) gauge invariant but violates the 0i and 02 number 
conservation. We emphasize that in realistic models, the leading order CP violating 
effects appear at one loop order or higher [20J and the magnitude of the coupling can 
be model dependent. However, as a first estimate (2.39) is sufficient. 

The above interactions can be used to model the classical generation of a baryon 



5 This effective four scalar vertex may be compared to the 4-fermion vertex one obtains after 
integrating out the W boson from the electro-weak Lagrangian. 
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asymmetry or what is commonly known as Affleck-Dine Baryogenesis [20J. Recall 
the three criterion for Baryogenesis are known to be (i) baryon number violating 
interactions, (ii) C and CP violation and (iii) out of thermal-equilibrium evolution 
[52] . The first criterion is present in the above model, and the third criterion is also 
present since the fields will be in coherent states and thus far from thermal equilibrium. 
The C and CP violation is not as obvious. One must consider the phase of the flat 
direction VEV $ = |$|e s as well as the phase of the complex coupling A = |A|e e . 
Either one of these phases may be removed at some fixed time by a field redefinition 
$ = e iA °$*, but the relative phase may not be removed, and this is the source of C 
and CP violation [20]. Our tree-level analysis is insufficient to establish this relative 
phase, and in fact, our toy model which only possesses two fields and a U(l) symmetry 
is too restrictive to allow for C and CP violation at all. Hence we simply insert this 
phase by hand. A careful discussion here of of the origin of the C and CP violating 
phase would require a long digression, but it is sufficient to note that these effects 
are present in the standard model, and can also be present in the MSSM and Grand 
Unified models. Our goal was to motivate the origin and magnitude of the above 



quartic coupling (2.39), and this has been done. 

To conclude, we return to the topic of the source of supersymmetry breaking, 
and its dynamics. The potential C B0 ft may be time dependent in the context of 
cosmology, and one should be sensitive to this fact |7J. For instance if supersymmetry 
breaking in the hidden sector is somehow tied to inflation, it is not guaranteed that 
m susy ~ TeV is the correct mass scale for the scalar particles during inflation or the 
early stages of reheating. Similarly, there may be competing mechanisms which break 
supersymmetry. An example of a competing susy breaking mechanism is provided in 
[3"2l [45] in the context of supergravity. These authors point out that the inflaton 
itself breaks supersymmetry by virtue of it having a VEV, and this breaking can be 
communicated to other fields through the Kahler potential. In this situation, there 
may be corrections to the masses of the fields of order the Hubble parameter H with 
the sign of these corrections being a free parameter. Such corrections, if present 
would presumably be the dominant contribution to the flat directions potential, and 
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thus a dominant contribution to the evolution of the flat direction during inflation and 
reheating [321118]. However there exist well-motivated models of supergravity such as 
"no-scale" supergravity [49J in which such mass corrections are absent at tree-level 
and thus supressed [50J. In the following we simply assume that the flat directions 
are lifted only by £ soft and by possible non-renormalizable terms of the superpotential 
mentioned in the previous subsection. Similarly, we assume £ so ft is time independent 
during the later stages of inflation and during reheating. 
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Chapter 3 



Quantization of Scalar Fields in 
External Backgrounds 

Here we introduce the formalism to describe quantum fields in external backgrounds. 
Some well known examples in which this formalism has been successfully applied are 
Hawking radiation [53J, the cosmological perturbations, [IT] Schwinger pair production 
[54] . and resonant preheating of the inflaton field after inflation [56]. In all these 
examples, there is a bosonic field which has acquired a vacuum expectation value. 
The VEV evolves classically and acts as a source for the other fields to which it is 
coupled. The coupling then may lead to a time evolution of the parameters of the 
Lagrangian for the other fields which is non-adiabatic and which may result in the 
production of quanta. 

This physical picture and the accompanying formalism are the main tools we apply 
to study the evolution of the flat direction field(s) during and after inflation. The 
cases we study in fact involve multiple scalar fields which mix and the mixing can 
provide more possibilities for non-adiabatic evolution [4TJ £>5]. In this section, we 
determine the differential equations which describe the quantum evolution of these 
multiple fields [41], and we explain how to extract the particle production from the 
solutions. We discuss how to manipulate a model's Lagrangian to extract the driving 
terms which appear in these differential equations. It is then possible to arrange the 
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driving terms into a simple quantitative measure of non-adiabaticity which generalizes 
the measure of non-adiabaticity commonly used in the single field case. To conclude, 
we develop a series solution to the equations, and determine the scaling properties 
of these solutions under changes in the energy scales of the model. The material 
here is thus focused on the differential equations and the techniques for determining, 
analyzing and interpreting the equations and solutions. The techniques presented are 
applied in later sections. 

3.1 Diagonalizing a System of Multiple Scalar 
fields, Bogolyubov Transformations, and the 
Heisenberg Equations of Motion 

We are interested in studying actions consisting of multiple scalar fields such as the 
following, 

S = J d A x^ (0 i)jU <f> U(t -faM^j) (3.1) 

where the mass matrix M 2 is time dependent and non-diagonal. Note that the actions 
we consider in later sections acquire the above time dependent form only after splitting 
the action into a classical time dependent background and quantum fluctuations away 
from this background. The fields fa will then represent the quantum fluctuations. For 
convenience, the above action is transformed to momentum space and also written in 
matrix notation with the fields arranged into iV-vectors X = {</>i,0 2 , , ,4>n}, 

S = J dr] d 3 k^ (x' f X' - X f fi 2 x) , fij = M 2 + k 2 5 tj (3.2) 

Additionally, the time coordinate is 77 which will later be identified as the conformal 



time, and time derivatives are denoted with a prime. The system (3.2) is studied using 
the canonical quantization and the Bogolyubov transformations to relate initial time 
Schroedinger operators to the Heisenberg operators. With the solutions to the Heisen- 
berg equations of motion, expectation values of observables may then be determined. 
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Specifically, we calculate the production of quanta, and we determine properties of the 
derived equations which permit a large production of quanta. The reader is referred 
to appendix [G] or to reference [41 J for further details of this formalism. 



To motivate the discussion of the multi-field case (3.2), consider first a simpler 
one-dimensional system, a harmonic oscillator with a time dependent frequency. The 
action is, 

'S'harmonicoscillator J dTj — i^x' 2 - u(r]) 2 x 2 j (3.3) 

In the Heisenberg picture, operators evolve with the same equations as their corre- 
sponding classical equations of motion. Thus the operators x and p, satisfy, 

x = p 

p' = -uj 2 x (3.4) 

The position and momentum operators x and p may be converted into the complex 
basis a and a) , 



x 



1 (a + a f ) 



2u 



p = z y2^-^) ( 3 - 5 ) 

which are the familiar raising and lowering operators that satisfy the commutation 
relations [a, a*] = 1. In this basis, the Hamiltonian H = \{p 2 + co 2 x 2 ) is written, 

H = u (tfa + 0=^(^+0 ( 3 - 6 ) 

where n = a) a is a number operator h\n) = n\n) which returns the principle quantum 
number n with n = being the ground state. Now, if the frequency uj is constant, 
then the state of the system will exhibit a trivial time dependence - if it begins 
in any particular energy eigenstate, it remains in this eigenstate. Transitions may 
occur when time dependence is introduced to the Hamiltonian, either by addition 
of new terms or by giving the frequency to time dependence. The latter case of a 
time dependent frequency is the one considered. The time evolution of the operators 
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must be determined from the Heisenberg equations of motion (3.4). To solve for the 
evolution, first parameterize the operators with a Bogolyubov transformation which 
will relate the raising and lowering operators a(rj) and 0^(77) at some time 77 to the 
initial Schroedinger operators at time r] \^ The Schroedinger operators are defined, 
a = a(rjo) and a) = aJ(r] ), and the Bogolyubov transformation is then, 

a(rj) = a(r))a + (3(r])(r 

a\ v ) = a*(rj)a! + /3*(r?)a . (3.7) 
The reason for this transformation is that the Heisenberg operators a and a) lead to 



the canonical form (3.6) in which a and aJ are the annihilation/creation operators 
of the physical eigenstates at any given time. The Bogolyubov transformation thus 
allows us to compute expectation values of observables at a given time once a solution 
for a and (3 has been determined. 

Recalling that the canonical commutation relations must be ensured at all times, 



one substitutes (3.7) into [a,aJ] = 1 and then requires [a, a)} = 1 to determine the 



following constraint on the parameters a and (5, 

\a\ 2 - \(3\ 2 = 1 . (3.8) 

However, performing the same substitution on [a, a] = does not lead to a constraint 
because the relation is satisfied trivially. The evolution of the a and (3 are subsequently 



determined by substituting (3.5) into the Heisenberg equations of motion (3.4) which 
yields, 

a = —iujol H a 

f3' = iu(3+—a (3.9) 
2uj 



It can be verified that these differential equations will preserve the constraint (3.8) 



It is thus only necessary to enforce the constraint on the initial values 01(770) ar| d 



1 Since our problems involve a time dependent Hamiltonian, the Bogolyubov transformations 
are time dependent. A space dependent Hamiltonian would have called for a space dependent 
transformation. In this sense, the Bogolyubov transformation is a general tool (see appendix |g|). 
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f3(r]o). The ground state of the system is defined by a(r/ ) |0) = 0, where a(r/ ) = a 
which means (3(r]o) = and a(r]o) = 1. The equations (3.9) can then be solved 
analytically or numerically, and with a solution for a(r)) and (3{rf), the time evolution 
of the operator is known. For example, putting the initial state as the ground state 
|0), the expection value of the number operator, h = a)a at some time is determined, 

(n(rj)) = (0\(a*a j + (3*a)(aa + f3a j )\0) 

= 13*13 (3.10) 



Note finally from (3.9) and (3.10), that if to' is zero, then (n) = const, and the state 
of the system will not change, which is the expected result. However, if > c<^, 
then the time dependence becomes nontrivial, and the system can make transitions. 
In other words, energy can be pumped into and out of the system by the effect of 
the rapidly changing frequency lo. The quantity ^- will be referred to as an adiabatic 
parameter since it quantifies the likelihood of transitions between the eigenstates of 

71. 



We proceed with the multi-field system (3.2) with the goal of reproducing the 



results of [41] . and we refer to the harmonic oscillator example where appropriate. 
The major difference between this system and the harmonic oscillator is the presence 
of multiple frequencies and thus multiple energy eigenstates for the system to populate. 
Additionally, the mass matrix is time dependent and not necesarily diagonal. One may 
diagonalize fi 2 with a time dependent rotation matrix C{rj), 

C T ( V )tf(r ] )C(r ] )=uj 2 (r ] ) (3.11) 

where u is the diagonal eigenfrequency matrix. Note the transformation C is time 
dependent, so we cannot generally perform this change of basis on the Lagrangian 
without introducing new kinetic structures, but if the mass matrix is evolving slowly 
enough, the column vectors of C will define the physical eigenstates of the system. 
Now define conjugate momenta in the nondiagonal basis IT = -§^r_ = X[, and define 
the fields and conjugage momenta in the diagonal basis X and fl as follows [4TJ . 

x = cx , n = cn 
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with 



d 3 k 
(2tt) 3 / 2 L 
d 3 k 



ikx 



%j{ri)aj + e ikx h*Ar))a 



(2tt) 3 / 2 L 

which, in the case of positive eigenvalues cut > 0, are equivalently written, 



(3.12) 



ft,: 



d 3 k 



(2vr) 3 / 2 ^2^ 
d 3 k . fuji 
'V ~2~ 



ikx 



di(k,r]) + e ikx a\(k,T]) 



e ikx di(k } ri) 



al(k,ri) 



(3.13) 



where the Bogolyubov transformation is now, 



di(k,r]) = a ij (k,r])a j (k) + P*Ak,r])aUk) 



a](k,r}) = a^kirjjaJ^k) + /3 ji (k,r])a j (k) 



(3.14) 



and for consistency between representations (3.12) and (3.13), the following relations 
will hold 



h 
h 



1 



'2uj 
iuj 



> + /3) 
-a + P) 



Note X and II satisfy canonical commutation relations [Xi(x),Uj(y)] = 5(x — y)5ij, 
and one may check using the orthogonality of C, that the nondiagonal fields sat- 
isfy the same commutation relations, [Xi(x),Hj(y)] = [CikX k , Cj m fl m ] = 5(x — 
y)8kmCikCj m = 5{x — y)5ij. The constraints on a and (3 coming from the commuta- 
tion relations [a*, at] = and [a,, a,j] = are respectively, 



aa' 



{3*{3 7 



a/3 1 - /? a 1 



1 





(3.15) 
(3.16) 
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These constraints reduce to the constraint (3.8) when N = 1. The Hamiltonian has 
the same form in either basis, 

H = J d 3 k^ (n f n + X j tt 2 X) (3.17) 

= J d 3 k^(&IL + X^uj 2 X^ (3.18) 
though in the diagonal basis, after normal ordering, has the well known form, 

H — / d 3 k uJihi with hi = a\(k)ai(k) (3.19) 

i J 

where the occupation number operator hi{k) has been defined. The diagonal basis 
Xi thus corresponds to the physical eigenstates of the system at any time during the 
evolution, ie. eigenstates of h. The Heisenberg equations of motion for the fields are, 

X' TT- 9C 

n' = -Q 2 X (3.20) 
and when expressed in terms of a and (3 take the form, 

a 1 = ( — iui - I)a+ - j \ (3 

(3' = (iuj-I)P+(^-jja (3.21) 

I,J = -(v^r^i^rv^J) , T = C T C (3.22) 
2 V Juj Juj J 



where the antisymmetric T and / matrices and the symmetric J matrix have also been 
defined [41J. Written using index notation, these equations are equivalently written, 

C^ij i^i^ij lik^-kj 2u~-@ij JikPkj 

flij — i^iPij — Iikflkj + 2i~.&ij — JikOtkj no summation on i, j 

The Gamma matrix also may be written in terms of the eigenvectors as follows, 

Tij =C i -C' j (3.23) 
40 



and thus the elements contain the information for the rate of change of the 
eigenvectors, 



As in the SHO example (3.4), the equations (3.21) preserve the constraints (3.15) 



and (3.16) at all times if the constraints are satisfied at some arbitrary time. Us- 
ing (3.14), the occupation number of the i'th state is determined, 

(fii(r))) ■■ 



{fl \a\(k,r])ai(k,r])\Q ) 
((3*f3 T ).. , no summation on i 



(3.25) 
(3.26) 



where |f2 ) is the initial vacuum, and where the Schroedinger operators ctj(/c) annihilate 
this initial vacuum. Note that the occupation number is not a good quantum number 
when lu 2 is negative. This is most easily seen by trying to compute the expectation 
value of the Hamiltonian which becomes imaginary. In this case, the modes are 
tachyonic, and a particle interpretation of the states is not appropriate. However the 
field operators are still well defined when written in terms of the mode functions h and 
h [4DJ introduced in (3.12). In the absence of a particle interpretation, observables 
such as a correlation function (4>(x)(j)(y)) may still be evaluated in the mode function 
decomposition (3.12). A calculation of this sort is performed in Section [4] for just 
such a circumstance. We will use the representation a and (3 whenever the particle 
interpretation is available. 



Finally, we mention that the Heisenberg equations (3.21) may be reduced to a 



smaller apparently equivalent set of equations by applying the algebraic relations (3.15- 



3.16). This reduction is shown in (3.5). 



3.2 Expressing Lagrangians in Canonical Form 



Typically the Lagrangians one wishes to study are not initially in the form (3.2), 
but must be put into this form. Our models are no exception, so we here present 
some useful techniques for manipulating a Lagrangian to extract the matrix Y and 
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the eigenvalues tOi which specify the Heisenberg equations of motion (3.21). Working 



in momentum space, the form of Lagrangian we are typically faced with contains 
additional mixed terms specified by a known antisymmetric matrix , 



L= l - (x V X' + X V KX + X ] KX' - X ] tt 2 X 



(3.27) 



where X is a vector composed of an arbitrary number of real scalar fields, (l 2 is a 
known symmetric matrix, and both (l 2 and K depend on time. Note that the matrix 
VL 2 is not the frequency matrix due to the presence of a nonzero K. To remove the 
mix terms proportional to K, one considers a vector Y which is related to X by a time 
dependent rotation, Y = TLX with the condition that the representation Y puts the 
Lagrangian in the canonical form. An analytic solution for 1Z is not always possible, 
but since 1Z is a basis dependent quantity, one expects it should not appear in physical 
results. We thus do not attempt to solve for 1Z. By expanding out the kinetic term 
\Y' T Y and comparing to (3.27), one may write (3.27) in the canonical form,[^] 



c = -y' 1 y' - ^Y T n(n 2 + k t k)tz t y 



(3.28) 



where one reads from above that the frequency matrix is fi 2 = 1Z(VL 2 + K T K)1Z T . 
Since the eigenvalues of VL 2 are invariant to a rotation of Q? , one can determine the 
the diagonal frequency matrixes by diagonalizing the known matrix (£l 2 + K T K) with 
a rotation U. Specifically, 



TZ[yi 2 + K T K)1Z T = 1ZUuj 2 U t 1Z t 



2 if K is not anti-symmetric, it can be made so by addition of suitable total derivative terms. 
Similarly, if the kinetic term posesses a time dependent matrix factor this may be diagonalized and 



made of the form (3.27) 

3 A relation between the kinetic term for Y and that of X is obtained by expanding, 



l Y lT Y' 
2 



1 



-X' 1 X' + X' T K T K'X + ~X T 1l' x K'X 



1 



\x' T X' + X' T KX 



1. 



-Y 1 KK 1 KK 1 Y 



where we have identified 1Z 1 7Z' = K, from which it also follows that 1Z' T 7Z' = K 1 K 
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It is also clear from the above that the previously defined C matrix (3.11) will have 
the specific form C = 1ZU. The matrix T is then computed, 



r = c T c = u T n T (nuy = u T n T n'u + u T u' 

= U T KU + U T U' (3.29) 
and thus both uj and Y matrices are determined without an explicit form for the rotation 



1Z. To summarize, given the Lagrangian in a form (3.27) with the two matrices K 
and Vt 2 specified, one may extract the relevant quantities for the evolution (3.21), 
namely the eigenfrequency matix uj and the matrix Y describing the rate of change of 
the eigenvectors. 



After having put the Lagrangian into one of the above forms (3.2) or (3.28), 
it may be that some of the fields which compose the vector X and which have 
nonzero entries in K and Vt 2 may in fact decouple. A decoupling may not always be 
obvious from examination of K or Vt 2 , but a computation of the / and J matrices can 
often indicate a hidden supression or decoupling. It may then be possible to perform 
additional rotations on X to decouple the fields and thus reduce the Lagrangian 
further. If the rotation is time-dependent, the transformation will induce terms in K 
and Vt 2 as discussed above, but the induced terms may be chosen to yield cancellations, 
decoupling, or other simplifications of these matrices. 

One specific case of interest is when the Lagrangian can be expanded in a small 
parameter. In our models the small parameter will be the ratio of two mass scales 
||. We will see later that the strong mixing (unsupressed in the expansion ||) in 
the Lagrangian dominates the Heisenberg Equations, and the weak mixing (mixing 
supressed by powers of jj) in fact contributes neglibly to the dynamics and may 
therefore be neglected. The fields which mix only weakly should then be decoupled 
with a small rotation, a rotation which deviates from the identity by small terms 
in positive powers of The strongly mixed parts of the Lagrangian will remain 
unchanged by such a rotation, but the transformation may be chosen to decouple the 
weakly mixed fields in the Lagrangian to a corresponding order in The decoupled 
fields may then be discarded. By this procedure, the remaining Lagrangian should only 
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contain the minimal field content, those fields which mix strongly and thus may lead 
to a non-adiabatic evolution. The reader is referred to Section [6] and Appendix [B] for 
specific examples employing this procedure. We now move on to further discuss the 
form of the Heisenberg equations and its interpretation and solution. 



3.3 Quantifying Non-adiabatic Evolution 



To gain intuition for the Heisenberg Equations (3.21), notice that there are two qual 



itatively different parts to the equations; a "unitary" part represented by the anti- 
Hermitian matrices (±iu — I) and a "nonunitary" part represented by the symmetric 
matrix — J). One can check that the unitary part preserves the total occupation 
number Tr[/?*/? T ], but the matrix (3* fi T is not preserved in general which indicates 
a possible conversion from one species of particles to another from this evolution. 
The non-unitary part will boost (or contract) both aa^ and (3*(3 T in (3.15) while 



keeping the difference invariant - preserving (3.15)-(3.16). This part of the evolu 



tion will change the total occupation number in the system. Note however, that the 
unitary part can not be neglected as the mixing between species effects the evolu- 
tion and so has physical implications. We make one final rewriting of the Heisenberg 



equations (3.21). By extracting factors of y/u, the equations are written 



a 



1 j4= 

to Juj 



■jj \ -I -=1—=) JujOl + Ju ( -£-r 

Juj y'UJ J \2uj Z 

U I / - / ' ) \fu)fi + Juj ( f ., •/ I \/^V! 

oj yio ) \2uj z \/uj v/ uj I 



up 

j= J V^a (3.30) 



The equations are now written entirely in terms of the following matrices, 

, A = 



1 1 

l±i—I—= 
uj Juj 



uj 



UJ" 



1 "4= 

UJ \ UJ 



where the matrix A has been defined. One recognizes the diagonal part of A as adi- 
abatic parameters from the single-field analysis. The off-diagonal parts of this matrix 
may also be interpreted as adiabatic parameters as follows. One expects nonadiabatic 
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evolution when either the oscillation frequency of the physical state or the field com- 
position of the state (its eigenvector) has changed in a time comparable to or faster 
than the period of one oscillation of the state — . The latter condition on the field 

r Ul 

composition is clearly relevant only if there are two or more fields from which the state 
is composed. Specifically, one may guess nonadiabatic evolution when \C-\ > uJi or 
specifically when, 



\ 



p x 2 
L i 



J 



This condition is not quite correct however, since in the case of degenerate eigenstates, 
the eigenvectors have no preferred directions in their subspace, and so the rate of 
change of these eigenvectors in this subspace is not physically meaningful. This feature 
is taken into account by the A matrix (as well as the J matrix) which has the form, 

A = 

Aij — t / \ V no summation (3.32) 




and it handles both limiting cases, 



lim A^ = , i 7^ j 



lim Aij = ( — ] , i^j 



*7 

LOj»UJi \ Ul, 



so for nondegenerate frequencies the condition (3.31) is asymptotically achieved. We 
thus take the elements of A as our adiabatic parameters. If any one element of the 
matrix is greater than unity \A^\ > 1, we expect nonadiabatic evolution. We note 
finally that large A^ are a necessary condition for net production of quanta, but not a 
sufficient condition. Thus the matrix elements A^ may be treated as useful indicators 
of particle production but not conclusive indicators. 
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3.4 Mass Hierarchy, Equipartition, Averaging, 
A Leading Order Solution, and Scaling 



In this subsection, we discuss the solutions of the Heisenberg equations of motion when 
the states of the theory have a hierarchy of two mass scales, the scale m ~ TeV of 
the flat direction, and the scale of the flat direction VEV ((pi) < Mp. For notational 
convenience we will denote the heavy scale simply by M and the relevant perturbation 
parameter is then ||. The quanta of the theory will also have masses of these two 
scales, but the nonperturbative dynamics, and thus the particle production will depend 
on the slow oscillations of frequency ~ m. If one wished, it is possible to consistently 
remove the heavy fields of the theory leaving behind an effective theory of the light 
states. However, this assumes the heavy states are always produced as virtual particles 
- there never being enough energy available to produce the heavy particle. This is not 
a good assumption. Our system of flat direction VEVs are extremely far from thermal 
equilibrium, so there is no notion of temperature to clearly indicate whether the heavy 
degrees of freedom will be "frozen" or "unfrozen." Additionally if nonperturbative 
effects are important, there can be processes with many light particles annihilating to 
create a heavy particle. We thus include the heavy states in the theory, though we 
use the hierarchy of scales as a perturbation parameter to approximate the dynamics 
of the full theory. Specifically, the approximation will allow us to numerically solve 
our system at a small value of the ratio || and to consistently rescale the solutions 
to correspond to a smaller value of this ratio. This is useful because when obtaining 
numerical solutions, there are limitations of machine precision which prevent one from 
taking the ratio || at the value of physical interest which is ^ < ~ 10~ 14 . 

In the following, we present a leading order solution to the Heisenberg equations of 



motion (3.21) valid during a parametric resonance of the heavy and light states, and 
valid after an averaging has been applied to the fast oscillations The above men- 
tioned scaling relations will follow once the leading order solution to these equations 

4 By fast oscillations we mean those with frequency ~ M. The slow oscillations have frequency 

~ m. 
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are identified. The first step is to determine the conditions on the equations (3.21) 
to yield a parametric resonance, and specifically conditions which take into account 
the expansion in fj. Recall, that after fixing the initial conditions a = 1 and (3 = 0, 
the equations are uniquely specified by the driving matrices I,J,uj and — which are 
in turn dependent on the matrices Y and to. The driving terms can be organized into 
a series in half integer powers of (tt). We will only require the leading order terms 
of these expansions to determine the leading order solutions of the equations. For 
instance, the heavy states have frequencies uja ~ M which are zeroth order, and the 
light states have frequencies uj a ~ m which are first order. Also, for the models we 
consider, both V and uj will evolve only on the slow time scale For instance, 

the rotations C which bring the system to diagonal form will evolve on this time 
scale. Thus C ~ m assuming the CV, are not further suppressed in any way, and so 
Yij ~ m. Now denoting the heavy fields with upper case indices A,B,C and light 
fields with lower case indeces a,b,c, the leading order driving terms which appear in 
the equations are tabulated as follows, 

uj a ~ \/m 2 + k 2 uja ~ M 

< ~ ^ ~ m (3.33) 



hb, Jab ~ rn I A b, Jab ~ m Im, JAb ~ VmM 
where the dependence on k is also shown assuming k < m. We also note by exam- 



ination of the definition for in (3.22), that a degeneracy in eigenvalues, uoi = Uj 
will lead to further suppression in J^, greater than the order listed above. 

We may also make a heuristic guess of the leading order solutions {«,/?}. An 
important consideration is that we are seeking a series solution in the specific case 
of non-adiabatic evolution when the occupation numbers of the states rtj = (P*jPij) 
are expected to be large. Assuming there is no suppression in the driving terms, all 
the light states should have comparable occupation numbers n a ~ n b at any given 
time. We do not expect the heavy states to have comparable occupation numbers 



5 In principle, V and ui could also evolve on the fast scale M , but such cases would correspond 
to a coherent classical motion of the heavy Higgs-like excit ations which are orthogonal to the flat 



directions. This possibility is discussed again in Section 



5.7 
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to the light states. Assuming there is sufficient mixing between the heavy and light 
states, we do however expect the energy density of any heavy degree of freedom to be 
comparable to that of a light degree of freedom, and thus loa^a ~ u a n a which then 
implies 

OL a Pa 1 [Ua~ fm 



(3 34) 

ola Pa] V^a V M 
This is essentially an assumption of "equipartition" though it does not depend on 
the presence of a thermal bath. Rather the "equipartition" is due to strong cou- 
pling between the states combined with non-adiabatic evolution for at least one of 
the states (45] . The heuristic guess of equipartition also allows one to make some 
statements about how the solutions rescale under a change in the parameters m and 
M. Specifically, so long as the non-adiabatic evolution does not vanish as one varies 



these parameters, we still expect the relation (3.34) to hold for the solutions. How- 
ever, the argument assumes there is a large production of quanta for the light states, 
and says nothing about this amount. To determine this one must directly examine the 
Heisenberg equations of motion which is done next. 

The analysis of the Heisenberg equations is facilitated by separating the two ex- 
treme time scales set by m and M. As shown above, the driving terms of the Heisen- 
berg equations all evolve on the slow time scale, though there are fast oscillations 
of frequency M introduced to the dynamics due to the presence of terms ujaola and 
ujaPa in the equations. Note that the slow oscillations will appear as constants on the 
time scale of the fast oscillations. Similarly, we expect the slow oscillations to be only 
sensitive to the average value of the fast evolving quantities. We should be able to 
formally average the fast oscillations from the Heisenberg equations leaving behind the 
slow dynamics which is relevant to the non-adiabatic evolution. The fast oscillations 
are averaged over by integrating the equations over exactly one or a few periods || to 
obtain a difference equation. The driving terms can be treated as constants over this 
very small integration region and they may thus be extracted from the integration. 
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Writing the Heisenberg equations in matrix notation, 

this equation is of the form X' = AX where A is an approximately constant ma- 
trix over short time intervals and X is a vector. The general solution is X{r]f) = 
e^f-Viixirii). From this one constructs a difference equation over a very small time 
interval of time At; which corresponds one fast oscillation period. The exponential 
may be approximated over this small time e AAv ~ 1 + AAr). The difference equation 
is X (j] + Arj) — X(rf) = AX(i])Ar] which for small time At] will simply be a differ- 
ential, except the functions X have been averaged over a small window of time. By 
this argument, we re-write our differential equations, 



d I a 




(3.36) 



where we have introduced the time averaged quantities a and (3. The equations are 
identical in form to our starting equations because the driving terms were treated as 
constants over the short time interval. Apparently nothing has changed except the 
notation, but the reason this averaging step has any meaning is that now a! ~ ma 
and P' ~ 111/3 where beforehand we had for some components a' Ai ~ Ma^ and 
(3' Ai ~ M(3ai- This subtle difference allows us to construct the leading order behavior 



of the solutions to equations (3.36) in the small parameter ||. The equations are 
first written using index notation again with upper-case A, B,C corresponding to the 
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heavy modes, 

O-'aB = —i^A®-AB — lAC&CB — IacO-cB + TT^PaB ~ JacPcB — JAcPcB 

ZU) A 

j3' AB = tUlA^AB — IacPcB — lAcficB + ~^~~OtAB ~ JaC&CB ~ Jac&cB 

uj' - 

&Ab = —iUA&Ab — lAC&Cb — lAc&cb + ^ @Ab ~ J Acficb ~ J Acficb 

Zlo a 



0Ab = ^A^Ab — lAcPcb ~ ^Ac^cb + ~ OlAb ~ J AC&Cb ~ J AcOicb 

Lu)A 



Oi' a B — —iUaOiaB ~ I a C a CB ~ lac&cB + W^PaB ~ JacPcB ~ JacPcB 

Zu a 



fi'aB — i^aftaB ~ IaCaCB ~ lacftcB + TT~ a aB ~ JaC&CB ~ JacOicB 

ZUJ a 

-/ - - - - oj' a - 

OL a b = —iu a Ol a b — IaC&Cb — lac&cb + ^-fiab ~ JacPcb ~ Jacficb 
fi'ab = ^ofiab — lacficb ~ lacftcb + ^-&ab ~ JaC&Cb ~ Jac&cb 



In the first four equations above we notice using (3.33) that for instance 
and so the left-hand-sides may in fact be neglected in the series approximation. From 
the second pair of equations, we may then infer there exists a solution with relative 
order {cxa, Pa} oc \J 1 ^ I {oia- l f3 a } which was guessed by our heuristic argument of 
equipartition (3.34). There are also many terms above which are subleading and may 
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be discarded. The resulting leading order equations are obtained to be, 

= -iuj A a AB 

= iuj A (3 AB 

= — iu A a Ab — I Ac a cb — J r Ac f3 c b 

= iuj A f3 Ab - I Ac f3 cb - J Ac a cb 

&'aB = —lac&cB + ~^~^aB ~ JacficB 

(3'aB — —lacPcB + W^-O-aB ~ Jac®-cB 



O-'ab — — ^a^aft ~ I a C a Cb ~ lac&cb + ^-fiab ~ JacPcb ~ JacPcb 

2u a 



P'ab — i^aPab ~ IaCaCb ~ lacftcb + -^~ a ab ~ JaC&Cb ~ Jac&cb (3.37) 

We should check that the leading order solutions to the equations do not exhibit fast 
oscillations, and in fact, this is guaranteed through the heavy modes {a A i,fi A i} being 
algebraically solved in terms of the light modes in the second pair of equations. In 
this sense, the heavy modes have been removed from the dynamics, though the heavy 
quanta are still allowed to be produced. We also learn from the above that {a AB , (3 A b} 
are zero to leading order. 

Finally, we may show the scaling of the leading order solution to the averaged 
equations a, j3 under changes in the parameters which was our original purpose. Sup- 
posing a solution has been determined for the Heisenberg Equations at some small 
value of the parameter jj. Then the scales are changed as follows m — > fim and 
M — > 7M for which the ratio || — > ^ is still sufficiently small. Using the scaling 
properties of the driving terms (3.33), one may infer rescaled solutions to the leading 



order equations (3.37) in terms of the known solutions. Table 3.1 illustrates the rela- 
tion. Notice in particular that because the light mode solutions a a j,(3 b j are invariant 
under the rescaling, their initial conditions are preserved, but this is not true for the 
heavy modes a A j,p A j. Having averaged over the fast oscillations, the history of the 
heavy states has been erased. This is consistent with the dynamics of the modes being 
attracted to an equipartition of energy. 
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known sol. 


rescaled sol. 


Scales 


m , M 


fim , 7M 


Solutions 


&ij i ftij i 


&aj 1 fiaj j "^a 


Functional dependence 

e.g. 6c k {r]) -> a knew (r] new ) 


7] , k 


n - n h — k 
'Inew .. 1 n new ,. 

H 1 



Table 3.1: With known solutions to the leading order Heisenberg Equations (3.37) 



at the scales m and M, one may infer rescaled solutions at the scales fim and 7M 
assuming both || and ^ are sufficiently small. The table shows how the solutions 
and their functional dependence change under the rescaling. The upper case index 
"A" corresponds to a heavy mode, the lower case "a" correspond to a light mode, 
and i,j correspond to either heavy or light. 
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To summarize, the solutions to Heisenberg equations have been shown to exhibit a 
clear scaling at leading order in the expansion ^. The scaling of the solutions is shown 
in Table 3T The scaling has been determined from a heuristic analysis applying the 
equipartition principle and from a direct analysis of the differential equations. The 
scaling will additionally be verified numerically in later sections over a range of ratios 
H and finally the scaling will be applied to extrapolate numerical solutions for small 
values of || to the very small values of ~ 1CT 14 of physical interest. A byproduct 
of obtaining the scaling relations was the physical prediction of equipartition of the 
energy during the non-adiabatic evolution. 



3.5 Reformulating the Heisenberg Equations 



In this section we present a nontrivial re-characterization of the Heisenberg Equations 



of motion as written in (3.21) in which the constraints (3.15 3.16) are built-in. The 



equations and the constraints again are, 

(-iw-I)a+(£- J)P 



a 



0' 



aa^ - f3*f3 T = 1 
a/?* - (3*a T = 



a 



(3.38) 
(3.39) 



where the dependent variables a and (3 are square matrices, uj is the diagonal eigenfre- 
quency matrix, / is antisymmetric and J is symmetric. We attempt to re-write these 
equations in quadratic variables n = (3* 1 '. Recall this is a hermitian matrix which 



yields the occupation numbers on the diagonal (3.26). We compute, 

n' = (3*' f + (3* f 



'-iuo - I)(3* 



UJ' 

2uo 



J a* 



(3 T + (3* 



(3 T (iuj + I) + a T [£--J 



'—iuj - I)n + n(iuj + I) + ( ^~ 



J) a *^ + (3*a T ^ 



\2oj 
J 



(3.40) 



53 



Examining this equation, it seems appropriate to also define the elements of the matrix 
a*p T = (3a^ as dependent variables (notice the constraint (3.38) is applied in 



m 



this definition). We compute the time derivative of m to check if the system of 
equations closes, 



Tl 



m = a p + a p 
(iuo - I)a* + 



J P* 



r (iuj + I) + a T ( — - J 



(iu - I)^^ + a*f3 A {iuj + 1) + 



\2u 



J\n+(n* + l) 



\2uj 
2uj 



J 



(3-41) 



where the constraint (3.38) was applied in the last step to substitute a*a = (n* + l). 
The system of equations (3.40 3.41 ) does in fact close with the new variables n and m. 
The original variables a and (3 are N x N complex matrices subject to two algebraic 



constraints (3.38 3.39). Analogously, the new variables n and m are iV x N complex 
matrices subject to structurally equivalent algebraic constraints n) = n and m T = m 
which makes n hermitian and m symmetric. The number of independent variables 
of the old and new representation are thus seen to match, and the two systems of 
equations are apparently equivalent. However, one must be careful because our new 
variables are quadratic expressions, and it is still not obvious whether some physical 
information has not been lost because of the nonlinear transformation. It may be 
possible to write a and [3 explicitly in terms of m and n and thus show a strict 
equivalence, but this computation is not performed here. 

It is worthwhile to express the equations in a Real representation using Real sym- 
metric matrices n R = Re(n) and m R = Re(m) and the Real anti-symmetric matrices 
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nj = Im(n) and mj = Im(m), and the system is then written in the form, 



{ a/ T \ ( 



n' R = n R I - In R + ojui - niu + f — - J j m R + m R ( — - J 
n, = n/J — In/ — cun^ + n R uj + J m; — mj I J 



V2cj / \2u; 
m' R = -umi - miuj - Im R + m R I + - J^j n R + n R J 



m'j = ujm R + m fi c<j — Irrij + mjl + .J^jrij + nj ^- J^j (3.42) 

One notable feature of these equations is that the initial conditions one usually consid- 
ers of vanishing occupation number are extremely simple - they are n = m = 0. There 
is only one term above in the third set of equations which provides a seed to drive the 
evolution of the system — j), and this quantity has already been identified in the 
previous sections as being responsible for particle production. 

We make a few final notes regarding the above reformulation of the Heisenberg 



Equations of motion. One point to be made is that the constraints (3.38 3.39) are 
essentially enforcing canonical commutation relations (which was discussed in the pre- 
vious section), and apparently the reformulated equations now have the commutation 
relations automatically enforced! To build upon this formulation, one might attempt 
to construct analogous equations for fermions in which the anti-commutation relations 
are incorporated. This is in the same spirit of the program begun in [41] in which the 
equations (3.21) were originally presented in a multi-field context. Also, the dimension 
of the system was never specified, so the equations are generally applicable. We note 
that the single field case of the above equations has been noted and applied already 
in |55]. 

Finally, while this new formulation may in fact be more efficient for numerical 



computation than the original formulation (3.21), it was derived only after the nu- 
merical results had been obtained and verified via (3.21). Conseqently, these extra 
computations were not required and so not performed. The material of this section 



55 



may thus be considered tangential to the thrust of the thesis. 
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Chapter 4 



Inflation and the Evolution of 
Scalar Fields in the Early 
Universe 

The goal of this section is to introduce inflationary cosmology, and to explain the 
process by which inflation converts the quantum fluctuations of a light scalar field 
such as a flat direction into observable classical perturbations of the field [15J. To 
an observer present at some time after the conclusion of inflation, the spectrum of 
perturbations for the field 4> which are larger than the observer's Hubble radius will sum 
up to appear as a homogeneous classical field or a vacuum expectation value (VEV) 
defined (4>) pa tch = Because a quantum process has generated the field, different 
observers in different patches of the universe will observe different realizations of $. 
We will consider patch sizes which correspond to the present observable universe, and 
we determine the variance of the distribution of $ for an ensemble of such patches. 
The computed variance will then be a measure for the typical amplitude $ obtained in 
any arbitrary patch, and thus a best estimate for the value in our observable universe. 

In computing the variance, we also take into account limitations to its growth due 
to a finite duration of inflation as well as due to the effects of higher order terms in 
the fields potential which will cutoff the growth of the VEV [33J. The main result 
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of this section is a characterization of the post-inflation initial conditions for the flat 
direction fields (fa) prior to their oscillation and subsequent decay during reheating 
which will be studied in detail in Sections [5]j7j A critical input to the generation of 
the above perturbations is the Hubble parameter during inflation Hj. We conclude by 
briefly discussing the COBE normalization of the inflaton's potential, and a bound on 
Hi set by non-observation of a gravitational wave signature in the Cosmic Microwave 
Background. 

4.1 The Classical Theory of Inflation 

We first review the basics of inflation, beginning with a list of the problems in the big 
bang cosmology that inflation resolves^ 

Horizon problem: The Cosmological horizon defines the maximum distance from 
any observer from which particles travelling to the observer are reaching the 
observer at some time t. It is thus the size of a causally connected patch in 
the universe at this time. The horizon problem is that approximately ~ 30,000 
patches that were causally disconnected at the time of the last scattering of the 
CMB photons, appear to have the same temperature to better than one part 
in 10~ 4 as measured by us today. There is no reason apriori that these regions 
should have the same temperature. Note: The Cosmological horizon is not to 
be confused with the event horizon of an observer in an accelerated reference 
frame. The event horizon marks a region of space about an observer from which 
an emitted light signal at a time t will never reach the observer. 

Flatness problem: The spatial curvature contribution to the Friedmann equation 
is measured by us today to be < 0.02, compared to the contribution from 
matter, radiation and dark energy [19]. This is consistent with a flat or nearly 
flat universe. A nearly flat solution is an unstable solution of the non-inflationary 

1 The list of problems that inflation resolves can be found in many references on cosmology. The 
list presented here was compiled from the reviews [67 68 69l I70j 
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cosmology. This means that even assuming the curvature contribution is as large 
as 1CT 2 today, it had to be smaller in the past. For instance at the Electroweak 
epoch, the contribution would need to be < 1CT 26 . Why is the curvature so 
small? One might postulate that the curvature is exactly zero by some symmetry, 
but a dynamic mechanism which drives this term to small values might also be 
a solution. 

Entropy problem: The entropy of our current observable patch of the universe can 
be estimated to be proportional to the number of CMB photons in our patch 
~ 10 88 . The entropy in this volume should remain essentially constant as one 
evolves the universe to early times. A legitimate question is "where did the 
entropy come from?" Did it come from one (or a few) explosive event(s) or is 
it coming in a more continuous process? In either case, one is challenged to 
provide a workable mechanism to create this entropy. One notes that both the 
horizon and flatness problems may be cast into the form of the entropy problem, 
and in this sense the entropy problem is not necessarily a different problem]^] 

Isotropy, homogeneity problems: Why is it that for each direction one looks, 
the universe appears mostly the same, and for different points in the universe 
the conditions also appear mostly to be the same? 

Problem of primordial inhomogeneity's: The CMB and large scale structures 
tell us that the early universe had density perturbations y ~ at the time 
of recombination and earlier. Where did these small perturbations come from? 

Gravitino Problem: The MSSM is a globally supersymmetric theory, but super- 
symmetry may also be promoted to a local symmetry in which case the SUSY 
algebra includes general coordinate transformations. Local supersymmetric the- 
ories thus include the graviton as well as the graviton's fermionic superpartner, 
the gravitino. The gravitino acquires mass during the breaking of local super- 
symmetry through a super-Higgs mechanism. The gravitino should also only 

2 In Guth's original paper [12J the words "entropy problem" never appear but his argument very 
clearly draws the connection between entropy and the horizon and flatness problems. 
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couple to matter gravitationally, so its lifetime may be very long [£] . Addition- 
ally, gravitinos will be produced thermally in the early universe. The gravitino 
problem is that an overproduction of gravitinos in the early universe leads to se- 
vere phenomenological consequences such as upsetting nucleosynthesis, diluting 
the baryon number of the universe, or in the case of a stable gravitino, leading 
to an overly large dark matter component. In order to avoid the gravitino prob- 
lem, the reheating temperature must be less then approximately 10 6 — 10 9 GeV 
to suppress the production of these particles where this range depends on the 
gravitino mass and the dominant decay channel of the gravitino [291 130], 

Monopole (topological defect) problem: Topological defects are stable field 
configurations of many grand unified field theories (GUTs), and defects such 
as magnetic monopoles are expected to form if the universe cools through a 
temperature of about 10 16 GeV which is the approximate unification scale of 
these theories. For Grand Unified theories to be viable, one needs to explain 
why these topological defects have not been observed^ 

The entropy problem is related to the assumption of evolving adiabatically from 
the present state of the universe to very early times. In Guth's paper [12] in which he 
introduces inflation, he connects the assumption of adiabatic expansion to the flatness 
and horizon problems, and he quantifies these problems as fine tuning problems of 
the initial conditions. He then points out, "...both problems could disappear if the 
assumption of adiabaticity were grossly incorrect" [12] ■ He then proposes a period 
of accelerated expansion as the solution. Going beyond the context of Guth's paper, 
here is how the above problems in cosmology are resolved: The acceleration, if it 
lasts long enough, has one effect of transforming a small causally connected region 
into a larger region of numerous apparently causally disconnected regions after the 
acceleration has ended. This solves the horizon problem, the acceleration also has 
the effect of driving any initial curvature contribution in the Friedmann equation to 

3 Though if the Gravitino is the LSP its lifetime would be necessarily infinite 

4 For an introduction to magnetic monopoles, see [6J. 
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very small values by the end of the acceleration period, thus solving the flatness 
problem. The acceleration will additionally dilute the particle density of all particle 
species, including topological defects such as monopoles or exotic species such as 
gravitinos. This solves the monopole problem and partly solves the gravitino problem^] 
Then, if the acceleration is isotropic, it will lead to a homogeneous and isotropic 
universe. However, the acceleration also has the effect of red shifting all fluctuations 
of the fields in the universe well outside any observers Cosmological horizon after the 
acceleration phase, so that the fields effectively reside in their ground states - the 
entropy per comoving volume has thus been driven to nearly zero. ^ We clearly 
do not want to wind up with an empty universe. Additionally, the consistency of the 
Friedmann equation requires some agent to be responsible for the acceleration. There 
are problems with identifying the agent of the acceleration with Einstein's cosmological 
constant; one cannot get the acceleration phase to finish, nor will there be any large 
entropy production. In Guth's paper, he models the agent as a Lorentz invariant false 
vacuum state. Other researchers later consider the agent as a scalar field [T31 E]. In 
either case, the entropy problem may be resolved via the large entropy generated by 
the decay of the false vacuum or by the decay of scalar field once the inflation has 
ended. Guth's contribution was to recognize that the monopole, flatness, entropy and 
horizon problems could in principle be resolved by this concept of inflation [T2]. It was 
soon after recognized that inflation may additionally resolve the problem of primordial 
inhomogeneities. 

We would like to see explicitly how the accelerated expansion is obtained. A 
first step is to write down Einstein's equations with the Friedmann Robertson Walker 
metric which is a good description of our universe during every stage of its evolution, 
including the inflation stage. The FRW metric has the form, 

ds 2 = dt 2 - R(t) 2 ( + r 2 dn 2 ) (4.1) 

\1 — KT / 

5 The gravitino problem is only partly resolved because the gravitinos may then be produced 
after inflation during reheating or afterwards if the reheating temperature is high enough. 

6 Note the current acceleration of the universe if it continues will have the same effect leading to 
a nearly zero entropy density for radiation for any individual patch. The radiation entropy will never 
be exactly zero as there is a finite entropy associated with the de Sitter event horizon |71l [72] 
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where the spatial part of the metric is both homogeneous and isotropic (or maximally 
symmetric), and where the k term specifies the spatial curvature of the universe. 
Substituting the FRW metric into the left-hand-side of Einstein's equations G^ v = 
8%GT^, and substituting the energy momentum tensor of a perfect fluid into the 
right-hand-side, one obtains the Friedmann Equation, 

* 2 =(l) 2 = f G '-F 

where p is the energy density of the fluid, and p would be its pressure. As mentioned 
above, the curvature contribution to the Friedmann equation, which is the term pro- 
portional to k above, is negligible today and even more negligible at earlier epochs 
under the assumption of adiabatic evolution, so we set k — 0, and do not mention 
the curvature term again. The FRW metric and the Friedmann equation we consider 
for the universe are then, 

ds 2 = dt 2 -R 2 dx 2 (4.2) 

The Einstein Equations additionally result in an equation for the acceleration of the 
scale factor, 

R -AtxG, . 

^ = ^— (P + 3p) (4.4) 

and there is also the statement of the conservation of energy and momentum for a 
perfect fluid in an expanding universe, 

d(pR 3 ) = -pd (R 3 ) (4.5) 

which is equivalently the first law of thermodynamics for the fluid. It can be shown 
that only two of the three equations ( 4.3|4.4|4.5 ) are independent. From (4.5), if we 



assume an equation of state p = wp, then one may solve for the energy density p in 
terms of the scale factor, 

R \ -3(1+™) 



P = Po ( ^ ) (4-6) 
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where p and R are the energy density and scale factor at some epoch (usually the 
current epoch). Now we explore the dynamics for different equations of state for the 
dominating fluid. The equation of state for a matter fluid is w = 0, as it may be treated 
as non-interacting pressure-less dust. The equation (4.6) then tells us p m atter oc -R -3 
which is interpreted rather simply as the expanding universe simply dilutes the particle 
number density by an inverse factor of the volume R~ 3 . For a relativistic particle 
species whose particles satisfy E « k, the equation of state is w = 1/3. The particle 
number density is diluted as it is for matter, and the wave-number k is also red-shifted 
by R^ 1 to give the final result p ra diation oc i?~ 4 . 



For both matter and radiation, the acceleration equation (4.4), implies a deceler- 
ating universe. We then notice the curious, yet remarkable fact that for a fluid with an 
equation of state w < —1/3, the acceleration equation gives a positive acceleration. 
Additionally, in the specific case w = —1, Eq. (4.6) tells us the energy density is con- 
stant. For a constant energy density, Friedmann's equation (4.3) requires the Hubble 
parameter to be constant, and the solution is an exponential expansion R oc e Ht . 

One simple way to achieve the equation of state w = — 1 is with Einsteins Cosmo- 
logical constant which gives an exactly constant vacuum energy p vac in Friedmann's 
equations. This mechanism is not viable for inflation since it is not clear how one 
would end the inflationary stage and begin a radiation dominated era necessary for 
the hot big bang. [] The choice that Guth made was a Lorentz invariant false 
vacuum state [12J which would mimic a cosmological constant, and yield the equation 
of state w = —1. However, he recognized there were problems with this picture with 
regard to the conclusion of the inflationary stage. The true vacuum would develop 
through bubble nucleation from the false vacuum, and most of the latent heat for the 
decay would be concentrated in the bubble walls. The entropy of the universe would 
thus need to be created through collisions of the bubble walls. One problem was that 
assuming enough collisions occurs, they produce too large density perturbations, and 
a second larger problem is that in order to obtain enough inflation to solve the horizon 



7 Additionally, there are no cosmological density perturbations generated in a pure de Sitter 
background. 
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and flatness problems, one requires the false vacuum decay to be so slow that the col- 
lisions would in fact rarely occur [12]. Shortly after Guth's paper, these problems were 
resolved in the general context of a scalar field evolving on a potential and achieving 
w ~ — 1 for the initial stage of the motion [Ol H]. We explain this latter case next 
as it is simpler and also the commonly accepted paradigm. Specifically, the inflation 
model we assume is the chaotic inflation model of Linde's [14]. We first write the 
action for a single real scalar field, the proposed inflaton, which is minimally coupled 
to gravity, 



S = Id X\f^g 



(4.7) 



1 

T 

Note the inflaton field will be denoted by \ nere ar| d in the remaining sections. The 
classical equations of motion for the inflaton are obtained to be, 

X-V 2 x + 3H X = -^- H=^ (4.8) 

dx R 

One may also show that the energy and momentum tensor for the scalar field is, 

t; = xr - Qx,« x a - v(x)j (4.9) 

which may be put in the form of the energy momentum tensor for a perfect fluid 

P=lx,«x, a +v( x ) 

t; = ( P + p )ujj v - P 5; with { P = lx, a x, a -v(x) (4.10) 

U>1 = 75fe VaNd When X, a X,a> 

with energy density p and pressure p and where defines the local four-velocity of 
the fluid. In the case where the scalar field has only time dependence, X-,aX-° = X 2 < 
we notice from the expressions of the pressure and energy density, that the equation 
of state p = —p is approximately obtained if the condition \\ 2 V(x) is satisfied. 
In this case, the Friedmann equation approximates, 

8ttG 1 

~3 w/ /v - 'x=xz 

We want the equation of state w ~ — 1 to be satisfied for as long a time as possible 
to solve the above mentioned problems in cosmology, so we impose another condition 



Hj^—V{ X i) when *'| « V{ X i) (4.11) 
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that the second derivative of the field also be small, or looking at the equation of 



motion (4.8), the condition is \ ^ H\ which then implies 3Hx = —V. Both 
conditions on x and on x can then be compactly written using two dimensionless 
parameters known as the slow-roll parameters e and r] defined as follows, 

x 2 ^v - ^l M2 p(y) 2 ( 4 - 12 ) 

X^H X - V = M^(jPj (4.13) 

with inflation being obtained when e C 1 and \r]\ <C 1. The convenience of these 
parameters is that they indicate when inflation is obtained given only information 
about the potential, and where the field is sitting on the potential. 

The above represents a basic picture of inflation. Before moving on to discuss 
quantum fluctuations during inflation, we return to quantify inflation's solution to the 
horizon problem and also to make some comments on the duration of inflation. In the 
context of inflation, the horizon problem is solved by putting the largest observable 
scales today within the event horizon during inflation so that these regions are in causal 
contact. The cosmological horizon length today is approximately H^ ay ~ 10 26 m ~ 
10 4 Mpc. Observers today are thus able to see fluctuations up to this approximate 
length but not greater. This length L > Htoday ma y ^ e ev °l ve d backwards in time in 
proportion to the scale factor]^] The minimum number of e-folds of the scale factor 
during inflation required to put this length scale within the event horizon is obtained 
from the expression, 

Htoday _ RRH R*i (4 14) 

U* DP \ ' J 

tij titoday J^RH 

We define Rrh as the scale factor at the start of a radiation dominated universe 
which is here assumed to also be the end of inflation (though see Section [5]). We 
define R} as the scale factor during inflation when our scale of interest has exitted the 

jR* AT 

horizon. We then parametrize ■^- J - = e efolds where N e f Q i ds has been defined which 



8 Fluctuations of this length and larger are frozen from the time they exitted the horizon during 
inflation to the time they re-enter the horizon today. They change only by being stretched with the 
scale factor. 
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is the number of e-folds to the end of inflation. For the remaining ratio above, we 



can use « where T today fa 10~ 13 GeV is the measured temperature of the 

CMB today. Taking the logarithm of our equation and reorganizing the fractions to 
put the unknown quantities T RH and H} together, one obtains, 

N efolds ~ 66- log f 7 ^-) (4.15) 



H* 

where we have used H to day ~ 2.33- 10~ 42 GeV. For the number of e-folds of observable 
inflation, we will hereafter use N e f Q i ds m 60. 

A final ingredient to our later analysis is the duration of inflation. One may 
determine this time by rewriting the slow roll approximation 3Hx = —V in a form 
that may be integrated, 

tfinal rXfinal ^JJ 

dt= ~Y rdx 
Assuming a quadratic potential, one obtains, 

_ g {Hjnit — H final) ^ 

ml 

A 

One must specify the initial and final Hubble parameters which are fixed by the initial 
and final values of the field. It is reasonable to choose Xinit ~ y ^-Mp since begin- 
ning with fields larger than this results in a situation where the quantum fluctuations 
of the inflaton dominate over the classical evolution, and inflation is then "eternal" 
[731 . With this choice, we have ~ ~ \ —, and using H ini t > H Hnal we 

1 a £1 final Xf V m * 

approximate the duration of inflation to be, 



m x V m x 



This concludes our development of the classical theory of inflation, and we now con- 
sider quantum fluctuations. 
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4.2 Quantum Fluctuations of Light Scalar Fields 
During Inflation 



Neglecting quantum fluctuations, the effect of inflation on any other field is to stretch 
the profile of the field and dilute its energy density to negligible amounts. However, 
the story is quite different when one considers the quantum fluctuations of the field. 
Specifically, a large quantum fluctuation of a field may exist for a short period of time 
in accordance with the uncertainty relations, and for fluctuations whose wavelength 
are far inside the event horizon, the field will see space as approximately Minkowskian. 
However the inflation will stretch this fluctuation, and once the perturbation's wave- 
length has stretched past the event horizon, the fluctuation can no longer propagate 
- it is frozen. Hence, if the horizon scale is small enough - which is to say if the 
Hubble parameter is big enough - then the mode will be frozen at a relatively larger 
value. This effect is known to occur for scalar fields with small mass, for the scalar 
perturbations of the inflaton field (which have a small effective mass), and also for 
the transverse part of the metric tensor perturbations. The result of the inflaton 
perturbations are the cosmological density perturbations which seed the structure for- 
mation of the universe and leave their imprint on the CMB. The results of the tensor 
perturbations are primordial gravitational waves which may also leave an imprint on 
the CMB. These latter two topics are not the focus of this dissertation, but they will 
be mentioned again briefly when we discuss phenomenological bounds placed on in- 
flation from CMB observations. The case we now consider is the first case mentioned 
above of a light test scalar field during inflation. 

We wish to study the evolution of the scalar field <f>(x, t) during and after inflation 
for some patch of the universe. If the patch width L is taken to be the scale of 
the cosmological horizon at some time after inflation, then an observer in this patch 
will have no knowledge of long wavelength fluctuations larger than his/her horizon. 

9 It is the equation of motion of the field which governs whether the field survives with non- 
negligible value at the event horizon. The equation of motion for the metric perturbations can be 
put in the same form as the equations for a scalar field. [75j 
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The sum of all these long wavelength fluctuations will appear to the observer as a 
uniform field which is a Vacuum Expectation Value (VEV) for the field. Additionally, 
the VEV is in general different for different patches owing to the random quantum 
fluctuations which seeded these VEVs. Sampling over the whole ensemble of patches, 
one determines a mean of zero (0) = 0, but the variance is not zero ((f) 2 ) ^ 0, and 
the variance will be the relevant measure of the amplitude of fluctuations^ In our 
problem, the patch of interest is the observable universe today, but our computation 
of the variance will be done at a time during inflation. Knowing that any fluctuation 
is frozen after it exits the event horizon during inflation and remains frozen until 
this fluctuation later returns within the Cosmological horizon after inflation, we may 
simply perform our sum at the instant during inflation when the scales of interest have 
exitted the event horizon. The largest observable fluctuations today were determined 
to have exitted the horizon during inflation approximately 60 e-folds before the end of 



inflation (4.15). Thus, our computation will be done at this time t 60 , and the value 
of L at this time will simply be the inverse Hubble parameter L = H$q. Finally, 
as should be clear from the above, when we determine the variance, the sub-horizon 
wavelengths will not be included in the sum. The variance will then be denoted 
(<t> 2 ) kL < 1 indicating that only wavelengths ^ > L are being summed. 

To perform the computation we change to conformal time dt = Rdr/, and we 
redefine the scalar field = ^ and work with the rescaled field ip. We express 
ip in the momentum space decomposition outlined in the previous section (and in 
AppendixjG]), which again is, 

/d 3 k 
j^j- 2 {e ikx h k (rj)a + e-*~hl(r,)cf) (4.18) 

where in our conformal coordinates, the momentum is related to the physical momen- 
tum by kphys = ]f ■ Using the above mode expansion, the variance is determined to 

10 Higher correlation functions could be appreciable if their are large cubic or higher order terms 
in the potential V (</)). For the moment, is assumed these terms are suppressed, so all statistical 
information is cont ained in the variance - the statistics are Gaussian. Non-Gaussianity will be 
discussed in Section 4.3 
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be (see Appendix [C]), 



d 3 k 



\h k 



(4.19) 



(2tt)- 

One must now determine the time evolution of the mode functions h^{rj) which follow 



from the field equations. The field equations for are the same as in (4.8), though 



now written in terms of rescaled fields cp and in conformal time units as follows 

R" 



tp" - v V + m 2 R 2 - 



R 



<p = 



(4.20) 



upon substitution of ( |4.18[ ), the mode functions h k (rj) will satisfy the equation, 

k 2 + m 2 R 2 (r]) 



K + u 2 ( V )h k = , uj 2 ( V ) 



~R 



(4.21) 



which is the equation for an oscillator with time dependent frequency as was studied 
in the previous section. To model the inflation in the slow roll regime, the metric is 
approximately de Sitter, 

-1 



ds 2 = R 2 (drj 2 - dx 2 



R 



(4.22) 



where Hj is the Hubble parameter during inflation which is approximately constant. 



The equations for the modes (4.21) with the above metric is, 

,2 



hl + 



k 2 + 



m 



1 

rf 



h, = 



(4.23) 



This may be put in the form of Bessel's equation, and the solutions h k {rf) thus deter- 
mined (see AppendixjC]). The solutions must be canonically normalized at early times 
(or large k) where the mode functions are well approximated by simple plane waves. 



Specifically, the equation (4.23) for the modes at sub-horizon scales, asymptotes to 
; = with solutions e ±tkr> . Referring to (4.18), we want h k {r]) — 1 



K 



k 2 hu 



early times for our fields to be canonically norma 
to be, 



-ikr) 



2k 



at 



ized. The normalization is obtained 



hk{v) 



,(-krf) +iY n (-kr})) , n 2 



9 
4 



m 
H 2 



(4.24) 
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These solutions are considered in the limit of small mass m < Hj where n 2 > 0, and 
large mass m > Hj where n 2 < 0. The variance (4.19) is converted back in terms 
of the physical momenta k 



phys 



R< 



and in terms of the original fields 



and the integration is performed over the dimensionless variable x 
wavelength limit x < 1. The results are, 

m 



<0 2 > 



<Hi 



2 (2n-5) r(TO )2 



f 1 dx (iffj) 
JO x x 



6|ra|(27r) 3 



for m -C Hi 
for m 3> Hi 



= <p/R, 
in the long 



(4.25) 



Notice the result (4.25) is scaling with the Hubble parameter as H] which is a measure 
of the vacuum energy density during inflation, Hi = ^pp^. Thus, the higher the 
energy scale of inflation, the larger the variance in our test field. However, in the case 
of a large mass, the variance (4.25) is suppressed by a factor - « — . The reason for 



the suppression is that a massive field can only propagate a distance of order — , so the 
field will have been damped to small values long before it ever reaches the de Sitter 



event horizon. In the case of a small mass, the integration in (4.25) is logarithmically 
divergent in the limit m — > 0, but with finite and small m, yields the result, 



^ W < H, 



Hj 



\47r 3 / 2 m 

To understand the source of the logarithmic divergence as m - 
the time evolution of the VEV during inflation. The results (4.25) were obtained 



(4.26) 

0, one must consider 



assuming a de Sitter metric which implies an infinite duration of inflation during 
which the mode functions hk(rj) are fixed points of the de Sitter evolution, and thus 
the variance at a given length scale is invariant in timej^] If one takes into account 
that there was an initial time t in which $ is zero over all the patches in the universe, 
the time dependence of the variance is seen to have the form [15] , 



3Hj 
•!>7r 2 m 2 



1 — cxp 



2m 2 



(t - 1 ) 



(4.27) 



which, at later times (t — t ) ^> asymptotes to the results (4.25 4.26) differing 



only by a factor of order one. In the case of early times t <C ^ the variance is 
11 See Appendix[c]and also the discussion in reference [16] 
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proportional to (4> 2 ) oc Hf(t — t ) which is also valid in the limit of zero mass. The 
time evolution is in fact that of a random walk for the field VEV. In the case inflation 
lasts a sufficiently long time, the logarithmic divergence in the variance would emerge 
as a result of a random walk over an indefinite number of time steps. However in this 
case, the VEV may instead be cutoff by higher order terms in the fields potential V($>) 
as discussed in [33J. Specifically, a reasonable guess for the flat direction potential is 



the series expansion defined in (2.37) which again is, 

oo 

\/($)=m 2 |$| 2 + K>3(|$|) , K>3(|$|) = ^|A n | 2 M 



n=4 



Mf 



2(n-l) 



(4.28) 



where the heavy scale in the expansion is assumed to be the Planck mass and the co- 
efficients are assumed A„ ~ 1. If the variance becomes large enough that the higher 
order terms in the potential are comparable to the mass term, then the variance 



results (4.25 4.27) become a poor description of the distribution. Higher order corre- 
lation functions become large, imparting potentially large deviations from gaussianity. 
We do not attempt to calculate these higher order n-point functions. However, the 
behavior of the distribution's width as a function of time may still be approximated. 
Assuming the flat direction VEVs begin at some time during inflation with negligible 
Vacuum Expectation Value $(to) = 0, the distribution will be a Gaussian and its width 



(variance) will grow according to the early time behavior of (4.27), and it will continue 
to grow until either the asymptotic value is reached or until the lowest non-vanishing 



higher order term in the potential (4.28) becomes comparable to the mass term. In 
the latter case, the distribution begins to acquire a nongaussian component due to 



these higher order terms and the evolution ceases to follow (4.27). One may guess the 



width of the distribution will remain essentially unchanged after this time, but with 
the shape deviating further away from the gaussian. By this argument, the width of 
the distribution is approximately, 



width 



m 



J-) 



(4.29) 



M P \nM Pi 

which was obtained by equating the lowest non-vanishing higher order term in the 
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potential to the mass term, and dropping factors of order one. An alternative ap- 
proach for estimating the width of the distribution is to split the field into a classical 
background and a fluctuation, = $ + 50 where the fluctuation 5<p is quantized. 
In this argument (outlined in the Introduction), one would expect the growth of $ 
to stop once the effective mass of the perturbations becomes the same order as the 
Hubble parameter, m c s ~ Hi so that the quantum fluctuations cannot propagate past 
the event horizon and thus the accumulation of the VEV will stop. Specifically, the 



equation for the fluctuations 5(f are the same as (4.20), but with an effective mass, 
where the mass m has been neglected. One then approximates m e g 



d 2 V 



using the leading non-vanishing term in the series expansion (4.28). This argument 
leads to the following bound for $, 



(fjhard 
bound 

M P 



nM F 



i 

n — 2 



(4.30) 



which is also bounded below M P but larger than (4.29) for a give value of n (see 



Figure 4.1 ). Either of the above estimates (4.29) or (4.30) for the width are acceptable 
for our purposes, though the second estimate is established when the mass of the field 
can be neglected against the higher order terms so it implies the distribution has 
evolved longer in time compared to the first case. Thus it is a larger bound. The 
latter bound may also be written in terms of the energy densities of the flat direction 
and inflaton as follows, 



P<I> bound 



Hi 



2 

n-2 



(4.31) 



where the quadratic term of the potential has been neglected in obtaining this. From 
this expression it is clear that p$ is bounded by p x which is expected. For instance, 
taking Hi ~ 10 _6 Mp, the ratio of the energy densities is maximized when n ~ 20 
10 -4 . V\ Note that since the quantum fluctuations cannot propagate 



where — ~ 

Px 

past the horizon once $bound ' s reacr, ed, it is essentially a hard bound beyond which 



The approximation in (4.31 1 gives 



in the limit of large n. This is expected, since in this 



limit, one expects the quadratic term of the potential to domin ate, a nd th us obtain p$ = 
Note the quadratic term, although neglected in obtaining (4.30) and ( 4.31 ) , is still present 
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larger values of the VEV are impossible to obtain. Finally, note that neither of the 
above two bounds on the distribution says anything about the shape of the distribution 
which is presumably very different than a gaussian in both cases. 



13 



There is one final consideration to be made which is that the Hubble parameter 
is not a constant during inflation, but decreases slowly as the inflaton field rolls down 
the potential. It is thus inappropriate to use a single value for Hj in the expressions for 



the variance (4.19 4.27). The calculations may be adjusted to take this into account, 



and if this is done, the infrared divergence noted above will return as the spectrum 
will have obtained a "red" tilt [69J. One must therefore impose an infrared cutoff 
to the corresponding momentum integral which physically can correspond to a finite 
duration of inflation. This calculation is not performed here. However, let us make 
some final crude estimates to take into account the duration of inflation. From the 
result ( |4.27| , one sees that a time At asymptotic > ^ is necessary for the asymptotic 
value ( |4.26[ ) to develop. Assuming zero temperature, the duration of inflation was 
determined for a quadratic potential in (4.16) to be At in fl at ion ~ — - Taking 



Hj to be the Hubble parameter at the conclusion of inflation, and comparing these 
two times, 



At 



asymptotic 



At; 



inflation 



V 2 

m 
10 17 



where the Hubble parameter during inflation has been fixed by the COBE normalization 
to Hi ~ 10~ 6 M P assuming a potential V(x) = km^X 2 ( see Section 4.5). The 



above estimate suggests that inflation may not last long enough for the asymptotic 



variance (4.26) to be reached. The variance will thus approximately grow as HfAt 



13 One example of a strongly nongaussian distribution is one which develops a minimum at the 
origin and a peak at some nonzero value of |$|. If the peak is sufficiently narrow, then in fact the 
field has been driven to an attracting surface (or fixed points if it is a real field) 



73 



through the whole evolution reaching the following value at the end of inflation. 

&RMS 



HfAt 



inflation 



El 

Mr 



Mr 



5/4 



10" 



(4.32) 



At the conclusion of Section 4.1 we determined that the Hubble parameter at the 
start of inflation was a factor . — larger than at the conclusion of inflation. If we 
were to take this into account carefully as discussed above, one would expect the 
above estimate of &rms to increase. Also note that this approximation assumes that 
the variance for $ at the start of inflation is fixed at zero, but this is not guaranteed 
especially in the context of eternal inflation [73J. Vacuum expectation values which 



have grown near to the two bounds (4.29) and (4.30) are thus not impossible. 



4.3 Estimating the Distribution of VEVs 



The previous discussion was to explain the mechanism of the growth of flat direction 
VEVs. We gathered that two limiting factors could be (a) the duration of inflation 



and (b) the higher order terms of the flat directions potential (4.28). It was finally 



noted that our estimates of the variance (4.19 4.27) did not take into account that 
the Hubble parameter changes from the beginning of inflation to the end of inflation 
making these results incomplete. However, our estimates for Q W idth and Abound are 
still applicable. 



Hereafter, we assume that higher order terms V n> s of the potential (4.28) are 



indeed present and that the distribution of VEVs has evolved so that the bounds (4.29) 



and (4.30) are appropriate. For instance, these bounds may be easily obtained if 
inflation is eternal. Both bounds are shown in Figure |4~T| as a function of the lowest 
non-vanishing term in the superpotential indexed by n (see ( 2.38[ )). 

Note that to keep the discussion in the previous section as simple as possible, the 
calculations were performed assuming a real scalar field. The MSSM flat directions 
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Figure 4.1: Estimates of the width (lower curve) and absolute upper bound (upper 
curve) for the distribution of flat direction VEVs as a function of the lowest nonzero 
term in the superpotential indexed by n (see (2.37) and (4.28)). The width and 
upper bound curves are obtained from formulas (4.29) and (4.30). In the above 
m ~ 10" 16 M P ~ TeV and H ~ KT 6 M P which is fixed by the COBE normalization 
for a quadratic inflaton potential (See Section 4.5) 



however are complex scalars. The above results (4.29) and (4.30) are still valid for a 
complex field. However a complex scalars may also be charged under some symmetry. 
In the case of the MSSM flat directions, this charge is typically the difference in the 
baryon and lepton numbers, B — L and the symmetry is the conservation of the B — L 
charge. Our estimates of the variance would apply to the modulus of the complex field. 
One expects the phase of the field to be randomized as well, though the distribution 
need not be uniform on the circle if C and/or CP symmetries are explicitly violated 
in the theory. 

Another consideration is that there are many MSSM flat directions, and if some set 
of them obtains large VEVs during inflation, then this may exclude another complemen- 
tary set of flat directions from obtaining VEVs as a consequence of the large induced 
mass. Thus, when the quantum fluctuations of the flat directions are generated during 
inflation, there will be competition between mutually exclusive flat directions. If infla- 
tion lasts long enough, some set of the flat directions may have obtained large VEVs 
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at the exclusion of a complementary set As the process is random, some patches 



of the universe may be dominated by one set of VEVs and another patch may be 
dominated by a complementary set, and these may smoothly interpolate with patches 
with small or no VEVs at all. A comprehensive analysis of the resulting distributions 
of VEVs would be a difficult task, and far beyond the scope of this work. However 
for our purposes, it is sufficient for us to have made an estimate for the approximate 



width of these distributions as quantified in Figure 4.1 



What is relevant to our models is not the distribution of VEVs at the conclusion 
of inflation, but the distribution at a later time when H < m, and in particular the 
distribution's tail. Our models will describe the evolution just after Ru^m- The motion 
from the end of inflation to this instant may in fact alter the shape of the distribution's 
tail. These issues are discussed next. 



4.4 Post-Inflation Classical Evolution of Flat Di- 
rection and Inflaton VEVs 

In the previous two sections, the focus was on the statistical distribution of VEVs 
obtained by a flat direction by the conclusion of inflation. In particular, both a width 



and absolute upper bound to the distribution was estimated in (4.29) and (4.30) and 



shown in Figure |4.1| Because this thesis is focused on the fate of very large flat 
direction VEVs, we will consider only the range of initial VEVs from this distribution 
falling between these two bounds and we study their subsequent classical evolution 
after inflation has ended. In general, the distribution of VEVs valid at the end of 
inflation will be transformed to a slightly different distribution at the instant when 
H ~ m. After this instant, the classical evolution has a relatively simple form. This 
state of the evolution is also the stage relevant to our models. We then discuss the 
early generation and subsequent conservation of the B — L charge belonging to the 
flat direction fields. Discussion of the quantum effects (decay) on the inflaton and flat 



14 since the effective masses of the complementary set have become large enough to suppress the 
propagation of their fluctuations through the horizon 
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directions after inflation is postponed to Section [5] which builds upon the description 
here. Elements of the following discussion can also be found in |20| [331 E21 174] , 
Recall from Section @]that inflation ends when the inflaton has evolved to a location 



on its potential when the slow-roll conditions (4.13) are no longer valid. Thus, the 
inflaton field begins to move in the sense x 2 ~ V(x), and begins oscillations according 
to the equation of motion (4.8) and the Friedmann equation. The spatial gradient 
term is negligible compared to the other terms initially, so the inflaton remains a 
homogeneous field over the horizon volume as its oscillations begin. The solutions are 
damped oscillations (see Appendix ID]) about x = with frequency m x ~ 10~ 6 M P , 



M P ■ ( + \ , /n ( M P 



From this, one calculates the energy density of the inflaton in terms of the scale factor 
R and one sees the evolution to be that of a matter dominated universe, 



p x » m x M> f ^ j (4.34) 

where Ro is the scale factor at the onset of inflaton oscillations or equivalently, the end 
of slow roll inflation, and where the prefactor has been dropped in order to simplify 
later calculations. The behavior of the energy density as a function of the scale factor 



"' is shown in Figure 



Ro 



4.2 



Note, the inflaton energy density p x will dominate over the energy density of the 
flat direction initially and thus drive the Hubble parameter. The precise evolution of 
the flat direction is determined through the flat direction's equations of motion and 
the Friedmann equation which are, 



$ + 3#$ + % = (4.35) 

(i) 2 = ¥(l $ l 2 + vw + Px) 

where we have assumed one complex degree of freedom for the flat direction $, and 
where V($) has the proposed form (4.28). To develop an understanding for how 
the flat direction evolves, consider two different scenarios for $ which lie on the 
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Figure 4.2: Evolution of the energy density of the inflaton (solid curve) and a flat 
direction VEV (long dashes) at the conclusion of inflation. The horizontal axis is 
Logarithmic in the scale factor where R is the scale factor at the onset of inflaton 
oscillations. The flat direction (long dashes) has been given an initial amplitude 
|$o| ~ lCT 2 Mp which corresponds to the the lower bound of Figure 4.1 with n = 10. 
Also shown is a case assuming the same potential but with the VEV having reached 
the upper bound of Figure |4.1[ In this latter case, the early motion of the flat 
direction (short dashed curve) is dominated by the higher order terms V n >3 in the 
potential (4.28), and so the initial energy density can be orders of magnitude higher. 
However, assuming the field is in slow-roll, one may argue that p$ should decrease 
slightly faster than p x as shown for example in the short dashed curve. As $ decreases, 
the terms V n> 3 in the potential (4.28) become comparable to the mass term, and the 
evolution is expected to smoothly transition between the two dashed curves (transition 
not shown). Note the analysis of this last case is somewhat heuristic and not by any 
means definite. Also shown are Rd x , the scale factor at the instant of the inflaton's 
gravitational decay and R eq , the scale factor at the instant the flat direction would 
overtake the inflaton (F.l) assuming $ has not already decayed. For a comparable 
analysis on which this figure is based see 
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two curves of Figure 4.2 The scenario with $ on the lower curve is simpler, so 
we begin with this. Recall from the previous section that in this scenario the mass 
term is comparable to the contribution of the higher order terms of the potential, 
m 2 $ 2 ~ V n> ^. At the conclusion of inflation, the field will be frozen, but will begin to 
move when the condition ~ H 2 is satisfied which is then the condition m ~ H . 
Thus the flat direction VEV is frozen until the scale factor R m as shown in Figure |4~2 



Another way to state this is that the field value is nearly frozen until the age of 
the universe H^ 1 is comparable to the period of one oscillation of the field mT 1 . 



The inflaton will still be executing damped oscillations (4.33) at this instant. The 
subsequent motion of $ will be discussed shortly. 



Next consider the scenario with $ on the larger bound (4.30) and thus on the 



higher curve of Figure |47lj The energy density of the flat direction in this scenario 
can be many orders of magnitude above the first scenario, but still sub-dominant 



to the inflaton energy density (see Figure 4.2). Also in this scenario, the VEV will 
begin to move immediately at the conclusion of inflation, since ~ H 2 is satisfied 
immediately. The exact motion of $ at the onset of inflaton oscillations depends on 
the nonlinearities of the potential. When the field moves, it is reasonable to assume 
that the energy density will decrease as fast or faster than the inflaton. For instance, 
if one assumes the field is in slow-roll (neglecting $ in the equation of motion), then 



this may be shown. A situation such as this is depicted in Figure |4T2] In this case, 
the energy density would decrease as fast or fast than the inflaton until the mass term 
becomes comparable to the higher order terms of the potential m 2 |$| 2 ~ V n> z($>). 



At this instant, the field would reside on the lower curve of Figure pTT} The evolution 
would then smoothly make a transition to the case of a mass dominated potential 
which was described above. During this whole phase of the evolution the assumption 
is that the field is in slow roll. What we draw from this analysis is that initial field 



values larger than the width (4.29) may evolve to approximately the same field value 
and energy density at the instant H ~ m. 



After this instant H ~ m, the equation of motion (4.35) is approximately linear 
in $, and the motion will be that of damped linear oscillator with a time dependent 
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damping term H 15 The solutions are similar to that of the inflaton, 



p \3/2 /p3 



$ « ( J $ m [A- cos(mt + <J r ) + zA 4 sin(mt + $)] + O I ) (4.36) 



which evolves with a matter-like energy density, 



- ™^ i^f) (4-37) 

and with A r , A, <5 r , Si being dimensionless coefficients of the evolution. Three of these 
coefficients are fixed by initial conditions and the fourth is fixed by the normalization, 

A 2 r cos 2 5 r + A 2 sin 2 5, = 1 (4.38) 



This evolution is shown for a sample set of parameters in Figure 4.3 where one sees 
the field executes a elliptic-like trajectory about the origin while spiralling to smaller 
values due to Hubble expansion. During this time of damped linear oscillations, 

the ratio of the flat directions energy density to that of the inflaton is fixed to the 
following ratio (assuming the inflaton and/or flat direction have not already decayed), 

2 

(4.39) 



Px \ M p, 

and we see that so long as $ m < M p , the inflaton will dominate the total energy 
density. If this condition is not satisfied, then the flat direction will overtake the 
inflaton even before the flat direction begins oscillations. As we have just discussed, 
the higher order terms in the potential, if present, will keep the flat direction VEV 
bounded below the Planck mass 



16 



One other element to the classical evolution of the flat direction is the charge 
associated with the motion of the flat direction. The charge of the flat direction is not 



5 The Hubble parameter does in fact depend on $ as shown in (4.351, but only very weakly. 
Strictly speaking, this is a nonlinear system, but it is approximately linear since H is driven mostly 
by the inflaton energy density p-^_ 



16 If however the potential ( |4.28| ) contains only the soft mass parameters and exactly zero higher 
order terms, and if inflation lasts a long enough time, then in principle one may obtain <f> m > Mp, 
and the flat direction may overtake the inflaton. This could result in a second stage of inflation 
driven by the flat direction VEV. However, this situation cannot lead to large enough cosmological 
perturbations, and is thus not viable. 
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Figure 4.3: Evolution of the complex scalar field $ for a sample choice of parameters. 



The parameters chosen are Ai 
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charge asymmetry as computed by formula (4.40) is e b—l 
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a standard model gauge charge but rather it will typically be the difference of lepton 
number and baryon number B — L^\ The field $ is then by definition a gauge singlet 
field. As was mentioned in the previous subsection 4.3 the phase of the flat direction 
VEV is also randomized during inflation, though the distribution would presumably 
depend on the specific breaking of C and CP symmetries of ones model. Additionally 



if there exist B — L violating interactions such as those discussed in Section [275] which 
involve the flat direction fields, then a B — L charge can be generated dynamically 
through the flat directions classical equations of motion. This is essentially the well 
known baryogenesis mechanism of Affleck and Dine [20J. If inflation has driven the flat 
direction VEVs close enough the lower curve of Figure [4~Tj then such nonlinear terms if 
present can in principle contribute a phenomenologically relevant B — L charge. Once 
the field has reduced in amplitude past H ~ m, the charge will then be conserved. 



17 Note however that in many Grand Unified models, the discrete B — L symmetry may be 
promoted to a continuous gauge symmetry U(1)b-l Kl- 
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The charge in this regime is determined from our solutions (4.36) as follows, 



n B -L = 

{m&JAtArCceQt 



S r ) 



and the number charge per co-moving volume (n B _ L R?) is thus conserved. Notice, the 
charge can take positive and negative values which will correspond to the handedness 
of the rotation of its trajectory in the complex plane. Also notice that given the 
constraint (4.38), the charge density is extremized when (5j — S r ) = ±ir and A+ = 
A r = 1, and in the case Si — 5 r = ±7r/2, the charge is zero. Defining the ratio of the 
B — L number charge density to the total particle density as the parameter e B -i, one 
determines, 

" L '~ L A i A r coB{5 i -S r ) (4.40) 



£b-l 



where this ratio was determined at R m . The quantity £ B ~l is the B — L number 
per particle of the condensate, and it is bounded between ±1. In the absence of 
B — L violation, the net charge should be zero and the ellipse will have collapsed to 
a line passing through the origin. In the above we have parametrized the motion, but 
we have not explained how the charge asymmetry is developed through the nonlinear 
terms. An example of such terms was introduced in 2.39| To see how a charge is 
developed, redefine $ such that its initial phase is zero and thus $ initially lies on the 



real axis. As discussed after (2.39) this field redefinition implies A is explicitly complex 
(not on the real axis). Also assuming $ = initially, the equation of motion for the 
imaginary part of <E> at the initial time reduces to 



Im 



-lm[A]$ 



3 

Re 



(4.41) 



and thus the field will begin to rotate off the real axis into the complex plane. Whether 
the rotation is clockwise or counter-clockwise depends on the phase of A. In our later 
numerical simulations we will instead take A to be Real valued and the initial phase 
to be nonzero which is an equivalent situation. The amount of baryon asymmetry 



(computed later in Section 5.2) is dependent largely on the initial VEV $ as we ll as 
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on the existence of interactions such as (2.39). At later times when the higher order 



interactions such as these become sub-dominant to the mass term, the built up charge 
will be conserved to a good approximation |20J. The B — L asymmetry may also be 
smaller or larger depending on how the inflaton and flat directions proceed to evolve. 
For instance if the flat direction energy density can overtake that of the inflaton, the 
baryon asymmetry can be order one. This situation would require very large sources 
of later entropy to dilute the B — L charge density to be phenomenologically viable. 
Alternatively, if the inflaton remains the dominant fluid, the B — L asymmetry obtains 
a suppression inversely proportional the amount of entropy released by inflaton decay. 
These issues will be discussed again in Section [5j 

4.5 COBE Normalization and the Energy Scale 
of Inflation 

Before continuing with our discussion of the evolution of the flat direction fields, we 
make a brief digression to discuss two very important inputs to our study from obser- 
vational cosmology. These are phenomenological constraints on the Hubble parameter 
during inflation, Hj obtained from the Cosmic Microwave background observations. 



As we have seen from the previous section in the results (4.19 4.27) the Hubble pa 



rameter or equivalently the energy scale of inflatior^ controls the amplitude of the 
perturbations one obtains for a light scalar test field present during inflation. It can 
be shown that the amplitudes of the perturbations of the scalar and tensor modes 
of the gravitational field (primordial density perturbations and gravity waves) are also 
controlled by the Hubble parameter in a similar way [171175] . A large value for Hi thus 
implies large cosmological perturbations which may be in conflict with observations of 
the CMB. The cosmological perturbations are not the focus of this thesis, but we shall 
briefly describe a constraint on Hi related to the tensor modes as well as a constraint 
between Hi and inflaton's potential known as the COBE normalization, related to the 



IS 



since during slow-roll inflation one has H] sa (3^2- J Px w ' tn Px ~ V(xi)- 
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scalar modes. 

The tensor modes are like the scalar test fields just discussed in that the only 
physical parameter which appears in their spectrum is the Hubble Parameter [171 175] . 
In this sense, a signature of primordial gravitational waves in the CMB is a direct 
measure of the energy scale of inflation. If one supposes the gravitational waves are 
a sizeable contribution to the CMB anisotropy, then the measured anisotropy places 
an upper bound on Hi which is approximately, 

Hi < 10~ 4 M P (4.42) 



where Hi is the Hubble parameter at the end of inflation _ This bound is very 
conservative as the primordial gravitational wave contribution to the CMB is expected 
to be sub-dominant. 

The scalar modes of the gravitational field are proportional to the perturbations of 
the inflaton field via Einsteins equation. The scalar modes may also be written in terms 
of the density perturbations ^ at the time of last scattering of the CMB photons, and 
these are related to the measured angular anisotropics of the CMB. Without going into 
these details, one obtains the COBE normalization of the inflaton's potential which 
is 



20 



- *' (4.43) 



y3/2 

- 1(T 5 (4.44) 



where the slow-roll approximation has been applied and where Hi and x are taken at 
the instant of 60 e-folds from the end of inflation which are to match the amplitude 
of the density perturbation on the largest observable scales today. The left-hand- 
side of this expression is proportional to the spectrum of perturbations, which again 
scales with powers of Hi, but the spectrum is also dependent on the shape of the 
potential which appears in the denominator of the second expression as V . The 



19 This argument can be found in |67j 

20 See [7| 69j or one of the many other references on inflation for a derivation of this relation. 
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COBE normalization (4.44) is a "normalization" of V(x) because, once the shape of 
the potential has been specified, then this equation fixes the value of Hj or equivalently 
fixes the energy scale. For example, the inflaton potential we consider is a quadratic 
potential V(x) = \m 2 x x 2 ■ Using the above formula, one obtains, 

m x ~ 1(T 6 M P (4.45) 

which is the value we will apply throughout. From this one determines the Hubble 
parameter at the end of inflation to be Hj ~ 10~ 6 M P . 

The implications of the COBE normalization cannot be underestimated (For a 
general discussion, see |5J). Once one chooses a form for the inflaton potential, the 
COBE normalization, sets the energy scale for all of the field dynamics and particle 
decays which follow inflation. In setting the energy scale it also suggests the particle 
theories which one may choose to study. For example, the quadratic inflaton potential 
considered is the simplest one can imagine in the chaotic inflation paradigm. This 
falls into the class of polynomial potentials of which many models of inflation can be 
mapped into [64, 7J. The COBE normalization then yields approximately the same 
energy scale Hi ~ 10 16 GeV for all these potentials [7] which is also the energy scale 



Mqut, at which the three gauge couplings appear to unify (see Figure 5.2). It is thus 
natural to consider supersymmetry and Grand Unified theories in combination with 
these potentials because one can then make predictions for the dynamics of these 
theories. 

It is worth noting that the quadratic potential is consistent with CMB observations, 
though other potentials such as a quartic potential are now disfavored [19J. We note 
finally that some models of inflation can allow for a lower energy scale of inflation. 
Many of these are well motivated and consistent with observations, but do not easily 
fit into the context of our discussion which assumes supersymmetry, B — L violating 
effects from the GUT scale, and large vacuum expectation values of flat direction 
fields which are generic to supersymmetric theories, but may not be present in other 
theories. We return to the discussion of the dynamics and decay of the inflaton and 
flat direction fields. 
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Chapter 5 

Reheating and Preheating 



In this section we describe reheating, which is the process by which the large entropy 
of the universe is generated through the decay of the inflaton and through the sub- 
sequent thermalization of the inflaton decay products. Reheating is thus the particle 
genesis for the universe. A target quantity in the study of reheating is the reheat 
temperature Trh which is the temperature of the inflaton's decay products once they 
have thermalized (assuming the inflaton energy density dominates). Another target 
quantity is the baryon asymmetry quantified by the difference of the baryon and lepton 
number densities divided by the entropy density, ng ~ £ . This quantity is expected to 
be conserved after the universe has reached a thermal state, and is related to the ob- 
served baryon to photon ratio today [76J. There are other possible observables related 
to reheating which will be mentioned at the conclusion of this section, but the focus 
here will be on Trh and t]b both of which are dependent upon the dynamics of the 
flat directions!]] 

To begin, a framework for determining T RH and i] B is outlined, based on the 
Hubble parameter during inflation, Hj the decay rate of the inflaton T x , and the 
rate of thermalization of the inflaton decay products T t herm- We then discuss three 
perturbative decay scenarios of the inflaton; decay into fermions, decay into scalars, 

1 While upper and lower bounds on Trh may be established, a reliable determination is not 
currently feasible due to uncertainties in cosmology and particle theory. In the future, if these 
uncertainties are overcome, a reliable determination of Trh might be possible 
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and a Planck suppressed decay. We describe the subsequent thermalization process 
of the inflaton decay products in the absence of flat direction VEVs, and then in their 
presence, explaining how in the latter case the thermalization may be delayed due 
to the effective vacuum established by the flat directions. However, this presumes 
the flat directions persist in a coherent state for a relatively long time. We note 
that the perturbative decay channels of the flat direction are suppressed by powers 
of the VEV [20J, so this presumption is not unwarranted, but indications are that 
a fast nonperturbative decay of the flat direction VEV may be possible [26J. To 
motivate the analysis of nonperturbative flat direction decay, we discuss the now well 
studied nonperturbative decay of the inflaton commonly referred to in the literature 
as "preheating" [56J. Based on this analysis we present a preliminary discussion of the 



nonperturbative decay of the flat directions. Figure [57T] shows how the elements in 
the above scenarios are related. The reader is also referred to the following books and 
reviews from which the background material here was gathered [5l IT6l 1671 1631 1641 170] . 



5.1 Determining the Reheating Temperature 

Here we continue to develop the picture of the post-inflation evolution of the inflaton 



and flat directions which was begun in Section 4.4 The emphasis in the previous 
discussion was on the classical evolution of these fields. The emphasis here is on 
quantum effects, namely the decay rates of these fields T x , and T$, as well as the 
rate of thermalization of the inflaton decay products, T therm . While the basic picture 
is developed, a simple framework for calculating the reheating temperature T RH and 
the baryon asymmetry is also established. The main difficulty and the main source 
of uncertainty is the reliable computation of the decay rates and thermalization rates 
and this is postponed to the next section. 

To begin, consider the following simple upper bound for the reheat temperature 
which assumes that at the end of inflation, the inflaton decays in less than one Hubble 
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Figure 5.1: Schematic of the model of reheating involving flat directions. For the i nf la- 
ton , decay may proceed via couplings to fermions and/or bosons through perturbative 
or nonperturbative mechanisms. The inflaton may also decay through gravitationally 
suppressed coupling (perturbative decay). Also shown is the flat direction evolution. 
The flat directions establish an effective vacuum into which the inflaton will decay 
and in which its products will thermalize, but the vacuum is effective only as long as 
the flat direction VEV remains coherent. Thermalizing interactions can be suppressed 
due to exchanged particles having very large effective masses. Additionally, the flat 



direction VEVs perturbative decay rate is strongly suppressed (see equation (5.24)) 
for the same reason that the mediating particles are very massive [20J. If the VEV 
is sufficiently long lived, its energy density may come to dominate over that of the 
inflaton. If this happens, the reheating temperature T RH will be determined by flat 
direction decay and thermalization of its decay products instead of that of the infla- 
ton. The figure shows the case of inflaton domination because our belief is that flat 
directions will decay quickly by nonperturbative mechanisms and thus removing its 
induced effective vacuum. Such nonperturbative decay mechanisms are introduced in 
this section and further developed in the remaining sections. 
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time. Specifically, all of the inflaton condensate decays and thermalizes with negligi- 
ble decrease in the total energy density in less than one Hubble time. Under these 
assumptions, and applying the upper bound on the Hubble parameter during inflation, 



Hj < 10 M p discussed in Section 4.5 we equate Hj, to the Hubble parameter during 



a thermal radiation dominated epoch during which time the temperature is related to 



H as follows (see (E.l)), 

/ 4^ X 1 / 4 7 

T =U^) ^ *(n=»+ 5 w (5.i) 

where g*(T) is the effective number of relativistic (massless) degrees of freedom at the 
temperature T with and being the number of bosonic and fermionic degrees of 
freedom. In the standard model g* = 106.75, and for supersymmetric GUTs g* < 1000 
from which we may obtain a convenient formula for the reheating temperature 



T RH ~ 0.WH rh Mp (5.2) 

where H RH is the Hubble parameter at the instant of thermalization. From this 
expression and the upper bound Hj < 10~ 4 Mp, one determines, 

T RH < 10 16 GeV (5.3) 

which is thus an absolute upper bound on Trh independent of any assumptions on the 
inflaton's potential or other interactions. An absolute lower bound on the reheating 
temperature is obtained from the fact that nucleosynthesis occurs in a thermal bath 
at a temperature of approximately 1 MeV [77], so 

Trh > W^GeV (5.4) 

It is remarkable that there exists approximately 19 orders of magnitude of uncertainty 
in Trh- However, the above are phenomenological bounds on Trh- If one makes 
assumptions on the particle theory applicable during reheating, this range may be 
narrowed. We will make some minimal assumptions on our particle theory in order to 
estimate the reheating temperature, the first of which is an assumption of a quadratic 
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inflaton potential \m?x 2 whose mass m x ~ 10 6 M P is fixed by the COBE nor- 



malization (4.45). This assumption also fixes the Hubble Parameter at the end of 
inflation, 



Hi ~ 10- 6 M p with V{x) = -m\ X 2 (5.5) 



1 

2 

so if the inflaton were to decay and thermalize instantly (within one Hubble time) the 
reheating temperature in this case is ~ 10 15 GeV which is six orders of magnitude 
above the temperature on needs to avoid an overproduction of the gravitinos |^] In 
light of this, one may look for mechanisms which would suppress either the inflaton's 
decay or the thermalization of its products. In particular, the role of flat directions 
in suppressing the thermalization of the inflaton decay products will be discussed 
in the following subsections. However, before making such assumptions, we erect 
a simple framework for determining T RH as well as the baryon asymmetry i] B with 
model assumptions entering through decay and reaction rates lY The only further 
assumptions we make in this framework are the quadratic inflaton potential, and the 
Friedmann equation. 

We determined in Section 14.41 the classical evolution of both the inflaton and 
a flat direction after the inflationary stage. The solutions are damped oscillations 



of the form (4.33) for the inflaton with frequency m x ~ 10~ 6 Mp and (4.36) for 
a flat direction with frequency m ~ 10~ 16 Mp. For both of these homogeneous 
scalar fields, there is an equivalent quantum interpretation as zero temperature (and 
thus zero momentum) condensates of scalar particles of mass m x and m respectively. 
These particles can have decay mechanism(s) into lighter species which we parametrize 
with the decay rates T x and T^. The decays are quantum effects which modify the 
respective classical equations of motion. Consider for instance, the inflaton. One 
estimates the number of \ particles in the inflaton condensate as, 

Px ~ 1 (X 2 + m 2 x x 2 ) (5.6) 
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v m x 2m x 

2 



V ~''m x x 2 (5.7) 



m x 



2 Recall in the context of the Gravitino problem discussed in Section 4.1 that the upper bound 
on the reheat temperature is in the range Trh % 10 6 — 10 9 GeV determined in 
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where the virial theorem has been applied in the last two steps to approximate the 
number density. The reaction rate equation for the x particles is 

where n x is the number density of the inflaton condensate. Substituting the expres- 
sions for the number density n x , we obtain, 

X+(3H + T x )x + m 2 xX = (5.9) 

which is the equation of motion Q4.8 ) including a decay term which appears on the 
same footing as the Hubble parameter.^] A similar equation may be derived for the 
complex flat direction field. For the moment, we consider just the inflaton decay based 
on the equation (5.9). When the decay rate is small compared to the Hubble parameter 
T x < H, the solution to the equations of motion is the result (4.33). Though when 
the expansion has slowed to a point where H < T x , the decays soon dominate over 
the coherent classical oscillation, and the condensate is soon depleted. This happens 
at a value of the scale factor we define as Rd x - Since the inflaton mass is orders of 
magnitude above the mass of most of its expected decay products, the products will 
be relativistic. For R > R dx , that is after inflaton decay, the energy density evolves in 
a radiation dominated universe (see Figure 4.2). In this regime after the initial decay, 



the decay products are not necessarily in a thermal state however [33], and there 
will be scatterings, particle creation and particle destruction to bring the system to a 
thermal and chemical equilibrium [321 IBS] . This thermalization rate is parametrized 
with a second reaction rate T t herm- If, at the time of decay, the thermalization rate 
is faster than the Hubble expansion, which implies T t herm > T x , then thermalization 
can be thought of as happening instantly at Rd x . The reheat temperature is then 
determined with (5.2) by setting H ~ T x in this expression. If on the other hand the 
thermalization process is slow, and T t herm < T x , then the reheating temperature is 
determined by setting H ~ T t h erm in this expression. Summarizing, 

therm ] (5-10) 



3 There are also quantum corrections to the inflaton mass in the equation of motion, but we 
assume these effects are already taken into account. 
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Reheating is thus determined by the slowest reaction rate. Note that we nowhere 
assumed the reaction rates T x and Ttherm were perturbative results. The rate T x in 
particular can equivalently parametrize a nonperturbative decay. 

The presence of flat directions with large vacuum expectation values can modify 
the above history of reheating through the value of inputs T x and T therm . Furthermore, 
if the flat directions have large enough expectation values, and if they are sufficiently 
long lived, the flat directions can come to dominate the total energy density. For this 
to happen, the inflaton must decay before or during the flat directions oscillations so 
that the ratio p$/p x is no longer constant but scales as R until a significant portion of 



the flat direction particles decay or up-scatter into a relativistic species (see Section 4.4 



and Figure 4.2). Referring to figure 4.2 the value of the scale factor at which the flat 



direction energy density would overtake that of the inflaton is obtained in (F.l) to be 

(5.11) 



R eq ( Mp\ 2 f m x ^ 2 



R \$oJ \T X 

Then assuming that the flat direction decay rate r$ is a constant, the condition 
for the flat direction to dominate is that R eq < Rdq> which may be rewritten (see 
formula Q), 

^ < f#V (5-12) 



r x \M P 



The region of flat direction domination is shown in the plane ($ ,r$) in Figure 5.4 



One may also wish to compute the reheating temperature in these cases when 
the flat direction dominates, but they imply a baryon asymmetry of order one, and 
thus require either an extremely small £b-l or a later "additional" source of entropy 
to dilute the B-L charge to the required amount. As the first case is not obviously 
natural and the second case adds complications to the simple picture of inflation, flat 
direction domination is not pursued here any further than establishing the parameter 



boundaries in Figure 5.4 
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5.2 Determining the Baryon Asymmetry 



The next task is to estimate the baryon asymmetry within the above context. It 



was noted in Section |2.4| that the MSSM flat directions may carry a B — L charge. 
It was also mentioned that such a charge may be generated by the early motion 
of the flat direction through Planck suppressed B — L violating interactions, and 
afterwards approximately conserved. While only a baryonic charge is strictly necessary 
for baryogenesis, it may be that a nonzero B — L charge is required [79J, though see 
[50] , The reason is that the sum of the baryon and lepton numbers B + L is violated 
by Electroweak processes known as sphalerons^] Howver the quantity B — L is 
conserved by Electroweak interactions, so any baryon asymmetry generated may be 
erased by these sphaleron processes when they are in equilibrium [79"] . 

In the following we will estimate the quantity nB ~ L which is the ratio of the B — L 
number density to the entropy density. This quantity is an adiabatic invariant since 
the B — L number per co-moving volume R?n B _ L is conserved soon after the flat 
direction oscillations begin and the entropy per co-moving volume R?s is conserved 
from Trh to the present epoch. The conversion of B — L to B by sphalerons is 
quantified approximately by the relation 



28 

n B = ^ u b-l (5.13) 

The baryon-to-entropy ratio today is (^) 6 ~ 10~ 10 which is obtained from the 
observed baryon-to- photon ratio j^] A challenge for any fundamental particle theory 
is to obtain this baryon-to-entropy ratio from first principles. There are many proposed 
mechanisms of baryogenesis [82J, but the mechanism which is intimately connected 
with the study of flat directions is that of Affleck and Dine (22]- We next compute 
the baryon to entropy ratio in the Affleck Dine baryogenesis scenario. 



4 For an introduction to sphalerons, see [16J. The first work on the generation of the Baryon 
asymmetry through sphalerons was [78J 

5 What is often reported in the literature is the baryon number to photon number ratio today 
r) — ^f- w 5 x 1CP 10 . The photon number is related to the entropy by s w 7n 7 which accounts for 

the neutrinos in the universe. 
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Recall from Section 14.41 the number density n b—l is computed from the classical 



trajectory (4.40), and once generated the B — L number per co-moving volume is then 
nearly conserved. The resulting baryon to entropy ratio depends on how much this 
charge is diluted by the entropy generated through inflaton decay and thermalization. 
Specifically, the entropy per co-moving volume is not fixed until after a thermal state 
is reached at Trh, so it should be computed at or after this instant. However, the 
resulting expression is the same if it were computed at the time of inflaton decay. The 
reason is that we have assumed the inflaton energy density dominates throughout, 
and the inflaton decay products behave like radiation before and after thermalization. 
The entropy density is related to the energy density through s ~ pj^ d from (E.2), 
so the entropy per co-moving volume is largely insensitive to when the inflaton decay 
products thermalize. The main approximation here is that the flat direction energy 



density is negligible compared to the inflaton energy density. As shown in Figure 4.2 
the dominance of the inflaton is not guaranteed. The entropy density and the B — L 
number density are both computed at the time of inflaton decay Rd x to obtain the 
desired ratio, 



(riB-LR 3 )dx 



£b-l 



-,3/4 



dx 



(-B-L 



El) 

Px) 

J*o| 
Mp 



1/4 



m 



Mf 



1/2 



m 



(5.14) 
(5.15) 
(5.16) 



where in the last step, the relation (4.39) has been applied as well as the relation 
T?,Md which follows from H dx ~ T x and the Friedmann equation. The 



Pd X 



baryon-to-entropy ratio is then estimated using the sphaleron conversion (5.13) and 
assuming a gravitational decay T x ~ ^f- with m x ~ lCT 6 Mp, m ^ 
obtains, 



l(T 16 Mp. One 



n B 
s 



£b-l 



5>0 

Mf 



10 7 



(5.17) 



Recall from (4.40) that e^-L quantifies the B — L number per particle in the flat 
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direction condensate with e b—l = ±1 corresponding to a maximally charged conden- 
sate. Also, our constraint that the inflaton energy density dominates translates to 



the constraint $ < Mp as shown in (4.39). Using both considerations, the above 
baryon to entropy ratio may still be significantly greater than one if the condensate is 
maximally charged and $ ~ M p . It has been argued by Linde [21] that generically the 
baryon to entropy ratio should not exceed one. In the example given in [21J, it is the 
thermal corrections to the masses of the states involved which result in the ratio being 
bounded below unity, but the example assumes a case of flat direction domination. It 
is not immediately obvious how one would generalize such arguments to our situation 



and assumptions. We do not attempt such corrections here. The above ratio (5.17) is 
sufficient to establish the leading dependence on the parameters $ and £b-l- Specif- 
ically, one concludes that to have a successful Affeck-Dine baryogenesis, one requires 
either £b-l <C 1 and/or $ ^ Mp given our assumption of a gravitational decay 
(and precluding any additional contributions of entropy such as from moduli decay). 
We emphasize that the baryon asymmetry, while sourced by early motion of the flat 
direction, appears to be largely independent of the flat direction's later dynamics, so 
long as the inflaton dominates the energy density. This result for the baryon to entropy 



ratio (5.17 ) and the above brief qualifications are thus the conclusion of our discussion 
of baryogenesis. 



5.3 Perturbative Decay and Thermalization of 
the Inflaton 

The perturbative decay of the inflaton and thermalization of its decay products was 
studied shortly after the first models of inflation were developed [83J (see also the 
text-books [TBI [84]). In the following, we briefly present the generic scenarios as 
well as some recent developments regarding the thermalizing interactions which follow 
inflaton decay [85]. 

In the absence of a flat direction VEV, the inflaton decays into a vacuum in which 
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supersymmetry is assumed to be softly broken at a scale m susy ~ TeV. We assume 
that there exist light states m^,m^ <C m x in this vacuum to which the inflaton can 
decay, specifically, spin 1/2 fermions via x ~^ or spin-0 bosons via x ~^ £6- The 
corresponding terms in the Lagrangian would appear as, 

£int = -gxi 2 - hx^ (5.18) 

where the coupling h is dimensionless, and g has dimensions of mass. The decay 
rates for these interactions can be computed, and assuming that the final states have 
masses m^,m^ <C m x , the perturbative decay rates are, 

r(x-«) = ^— . T( X ^^) = —^ (5.19) 
oirm x o7T 

The form of these rates could also be obtained from dimensional analysis. In the 

case these rates are gravitationally suppressed, one has g ~ m x and h ~ 

resulting in, 

Tfi 

r G ™, ~ (5.20) 

p 

which could also have been obtained by dimensional analysis and independently of the 
spin of the final state particles. The gravitationally suppressed decay would proceed in 
approximately — ~ 10 13 oscillations of the inflaton. To determine T RH assuming 
the gravitational decay, the thermalization rate r t h cr m must be computed. We use the 
results of [85J in which it is shown that a thermalization process can be expanded in a 
series involving the coupling a as well as the small quantity which quantifies how 
effective the process is at transferring the high momenta inflaton products k ~ m x to 
the reheat temperature. It is shown in [85J that the dominant thermalization processes 
after perturbative inflaton decay are 2 — > 3 inelastic scattering with a rate, 

r^ 3 ~«(^) 2 r 2 _ 2 (5.21) 

\ 2 / \2 
a \ I R 



r ^~¥ , -={i) [itj ( 5 - 22 ) 
where T 2 ^2 is the elastic scattering rate, a is the corresponding elastic cross-section, 
and n x is the number density for the inflaton condensate (which should be of the same 
order as the number density of its decay products). 
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Note that the thermalization rate T 2 ^ 3 is suppressed by three powers of the gauge 
coupling a 3 , but then enhanced by the factor which quantifies the processes 

efficiency in momentum transfer. One may determine the the thermalization rate and 
reheat temperature by a self-consistent method by first assuming r thcrm > Y x so that 
T RH ~ ^JT X M P , and then determining the conditions on the parameters for this to be 



true. By this method, one requires the gauge coupling to satisfy or 1 > 100 (see (F.3) 
and reference |26J). Keeping in mind that this is an order of magnitude approximation, 
this condition on a may be compared against the expected values of the gauge cou- 



plings, shown in Figure 5.2 where one would presumably choose the renormalization 
scale at or near a momentum transfer of order Q ~ Trh- Assuming the condition on 
the gauge couplings is satisfied, the thermalization happens immediately after inflaton 
decay and the reheat temperature is 

T RH ~ 10 9 GeV gravitational decay, w/o flat direction VEVs (5.23) 



Again, this temperature corresponds to the gravitationally suppressed decay (5.20) in 



the absence of flat directions. It will be a reference case in the remaining discussion. 

5.4 Perturbative Decay and Thermalization of 
the Inflaton With Flat Directions 

Next we consider the effect of flat direction VEVs on the inflatons decay rate Y x and on 
the thermalization rate of the decay products Y therm- During inflation, a subset of the 
MSSM flat directions have presumably accumulated a large condensate of particles 
in their zero momentum state (these are the VEVs). It is useful to think of this 
condensate as an effective vacuum. From Section |2j it was noted that this vacuum 
typically breaks all the standard model gauge symmetries, and thus the spectrum of 
particles is quite different than the unbroken vacuum. In particular, the induced particle 
masses are of order of the obtained VEVs which may be nearly at the Planck Mass. 
One notes however that the VEVs evolve slowly in time with frequencies ~ ?n susy as 
discussed in Section 4.4 For the inflaton whose mass is much greater than m susy , 
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Figure 5.2: Running of the gauge couplings in the standard model and the MSSM. 
The inflaton decay products are expected to thermalize through gauge interactions. 
With our assumptions, one expects the typical momentum transfer of the interactions 
to be T RH ^ 10 9 GeV from (5.23). Thus the relevant gauge couplings can be read 
from the righthand plot at Q ~ 10 9 GeV. Note also the apparent gauge coupling 
unification in the MSSM at the scale M GUT ~ 10 16 GeV. Figure obtained from [76] 



the vacuum will look stationary. However, for the light state^J of mass < m susy 
the vacuum will be moving relatively fast assuming the particle momenta is small, 
k < m susy . There can be nonperturbative effects for processes involving these states, 
some of which will be focused on later. 

However for the inflaton and for its relativistic decay products, the vacuum appears 
nearly static (adiabatic) but it contains a large number of heavy MSSM states. We 
note two effects of these heavy states. The first effect is that the perturbative decay 
of the inflaton into these heavy states will be kinematically blocked when > m x . 
However, if there is at least one light state in the spectrum for the inflaton to decay to, 
and assuming the coupling is not suppressed, perturbative inflaton decay should not 
be suppressed. Even in the case with all gauge symmetries broken such light states 

6 The light states will be the quanta of the flat directions or their superpartners as well as any 
other MSSM quanta which did not acquire masses through the Higgs mechanism 
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must exist. These will be the flat direction states or their fermionic superpartners 
which could be the quarks and leptons of the standard model or the Higgsinos if the 
Higgs fields are the flat direction VEVs. The presence of these states is guaranteed by 
supersymmetry, and the supersymmetry breaking can only lift the masses to the TeV 
scaleQ The second effect is that the gauge bosons have acquired the same large 
masses < Mp and this will suppress the thermalization processes which are assumed 
to be mediated by gauge bosons. Of course, if there are any remaining unbroken 
gauge symmetries, one expects the remaining gauge boson(s) to efficiently mediate 
the thermalization of the particles they are coupled to. 

Hereafter, it is assumed (i) the inflaton decays gravitationally implying a pertur- 
bative decay, (ii) that all gauge symmetries have been broken and (iii) there do exist 
light MSSM states into which the inflaton may decay (namely the flat direction states 
or their fermionic superpartners). Note the 2 — > 3 thermalization rate determined by 
the arguments of [85] relied up the existence of massless (or light) states which could 
be emitted as the third particle in the reaction. Assuming the existence of such states, 



and under the above three assumptions, r x is given by (5.20), and the thermalization 
rate is obtained by replacing a ~ in (5.22), and recomputing. In particular, for 
reheating to occur at the instant of inflation decay, one requires T 2 ^3 > T x , which 
translates to the condition, 

, 1 fm y \ 5 ( M P \ 2 ( $ r 
a 6 > 1 x 1 ' 



100 \Mp I \ m J \M J 



which has been determined in the Appendix (F.4). This is rewritten as a condition on 



3/2 



$0< licPJ 10 ~ 3Mp 

We may also determine the reheating temperature in the case this $ does not satisfy 



the above constraint. The expression for T RH in this case is determined in (F.5) and 



shown as a function of $ in Figure 5.3 It is again emphasized that these results are 



7 Only the Higgsinos should have mass ~ TeV . The quark and lepton masses should be that 
measured in the standard model obtained by the Higgs mechanism. 



99 



order of magnitude estimates, but the supression shown in the figure is apparently a 
mild one given our assumptions. 



5.5 Perturbative Decay of Flat Directions 

The suppression of T therm through the presence of the flat direction condensate as well 
as the case of flat direction domination both assume the condensate remains coherent 
for many inflaton oscillations. This is partly justified by making a perturbative estimate 
for the decay of the flat direction quanta into the effective vacuum established by 
the condensate. One expects the decay to be mediated by particles with induced 
mass of order of the flat direction VEV, and to decay into light fermions and bosons 
m^,m^ <C m susy of the effective vacuum. From this, we estimate a tree level decay 
rate 

If one takes into account the time dependence of one may determine the decay rate 
of the flat direction T p ^ rt which is the above rate T^ rt at the instant when H ~ T^ ert . 
Assuming the decay happens after inflaton decay, the rate has the form, 

4/5 / \ 1/5 

Ki>,J Vr 



X 



rr*=(^ rn (5.25) 



which has been determining in (F.6). The flat directions perturbative decay rate as 



specified by this formula is plotted as a function of |$ | in Figure 5.4 One sees that 
for |$ | ~ lCT 2 Mp, the flat direction energy density will overtake the inflaton. 

One reason why the perturbative result should be viewed with skepticism was 
mentioned above and this is the vacuum (condensate) is evolving on the time scale 



^■susy 



which is also the mass of the flat directions^] The vacuum thus does 



8 A second reason the perturbative result may be incomplete is that one expects the inflaton 
decay products to scatter off the flat direction condensate, and this may deplete the condensate. 
The interactions would be suppressed by powers of <!>; since the mediating particles are heavy, but 



it is not obvious the suppression would be as severe as in ( 5.24 ) since the kinematics are different - 
the inflaton decay products have an energy on the order of Trh instead of m susy . We do not pursue 
this decay mechanism here beyond noting its presence. 
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Figure 5.3: Approximate behavior of Trh as a function of $ assuming gravitational 
decay of the inflaton and suppressed perturbative 2^3 interactions which control 
the thermalization. The curve is determined from formula (F.5) which is obtained 



by the methods of Section 5.3 The reduction of the reheating temperature is less 
than an order of magnitude. The temperature 10~ W M P corresponds to thermalization 
occurring instantly after inflaton decay. 
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Figure 5.4: Flat Direction Domination: The shaded region above defines the region of 
the parameter space ($ ,r$) in which the flat direction energy density will eventually 



dominate over that of the inflaton. The region is determined through formula (5.12). 
In the case of a perturbative decay of the flat direction condensate, the decay rate of 



the condensate T$ is shown as the solid line from formula (5.25). In the case the flat 



direction dominates, its decay-rate becomes independent of $ c 
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not appear static to the flat direction quanta. An adiabatic vacuum is an implicit 



assumption in the perturbative computation of the decay rate (5.24) for instance. 
Hence the perturbative results involving the flat direction quanta may be completely 
wrong. Nonperturbative effects may in fact dominate, and perhaps facilitate a rapid 
decay of the flat direction even before the inflaton decays. Nonperturbative effects 
are discussed next; first for the inflaton, and then for the flat directions. 



5.6 Nonperturbative Decay of the Inflaton: Pre- 
heating 

The above estimates for the reheat temperature assume that T x , T t herm, and are 
obtained via perturbative processes in which the coupling is assumed small and the 
number of quanta involved in any reaction is also small. This assumption on T x is 
indeed good in the case of a gravitational decay of the inflaton, and so the pertur- 
bative estimate is appropriate. In fact, except for the content of this subsection and 



Section 5.3 we universally assume a gravitational coupling for the inflaton and thus a 
perturbative inflaton decay. However, it is worthwhile to consider the non-gravitational 
decay mechanisms which may be nonperturbative. The nonperturbative decay of the 
inflaton will prove to be a useful example for later studying nonperturbative flat direc- 
tion decay. 

The perturbative estimate for the decay of a species of particles is not appropriate 
if the relevant coupling is large. However, even when the coupling is small, if the fields 
involved have acquired VEVs, the perturbation expansion in the coupling may still 
break down. Essentially, the probability amplitude for a scattering with a condensate 
particle is very high due to the large occupation number. In this situation, one splits 
the field x i nt0 a classical background (x) and quantum fluctuations 8x, 

X = (X)+ S X (5.26) 

The classical background, which describes the condensate, is assumed known, and the 
computations done using the fluctuations 5x will not suffer from the breakdown of 
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perturbation theory present in the original field description. One must then be aware 
that any observable computed in the new field description may in principle correspond 
to a large number of quanta in the original field description. For example, one must be 
wary of making kinematic arguments which rely on there being a small number of initial 
state particles involved in the decays. These kinematic arguments can easily fail when 
the coupling is large and/or the fields have acquired VEVs. The computations should 
then be done using the above new field description using the methods of Section [3] 
and the results may then be interpreted in both frameworks. 

It has been shown that nonperturbative effects can be relevant to inflaton decay. 
These effects are referred to in the literature as "preheating" [56]. There are at least 
two striking consequence of preheating which are (1) that particles with masses orders 
of magnitude larger than the inflaton mass can be produced through decay, and (2) 
that the production of quanta can occur through an exponential growth. We will work 
through one example of nonperturbative inflaton decay following the development of 
[To] in which both these features will emerge. In the next section, we discuss the 
possibility of a nonperturbative decay for the flat directions, as it is the objective 
of this thesis to show that nonperturbative effects can also be relevant for the flat 
directions. 



We consider inflaton decay through the scalar coupling shown in (5.18), and we 
perform a perturbative calculation first from which the exponential growth of quanta 
is apparent even in the case of a very small coupling. The physical effect is that of 
Bose-Einstein condensation. In the decay \ ~> ££. the momenta of the product £'s 
are each \k\ = ^e. If one neglects the redshift due to Hubble expansion, one can 
see the occupation numbers in this momentum mode will quickly become large, and 
Bose-Einstein effects should be taken into account. To quantify the production, we 
write down the rate equation describing the growth in the number density of the £ 
particles which is, 

1 d(H 3 n £ ) 

-p-r — - — : — = 21> x 5.27 

i? 3 dt x x K J 

which is to be compared to the corresponding equation describing the depletion of the 

I i 2 

inflaton field (5.8). We already presented the decay rate T x = in the absence 
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of Bose-Einstein effects (see (5.19)). Now we must obtain the correction due to 



Bose-Einstein effects. The matrix element for the process x — * 6£ is. 



l,nk + l,n^ k + l\a^a_ k a x \n x ,nk,n-k) = (n^(k) + l)' 2 n x S' 6 (k{) (5.28 



\2„ r3/ 



and the matrix element for the reverse process — > % is. 

|(n x + l,n k - l,n_ k - l|a k a_ k a x |n x ,n k ,n_ k )| 2 = n 5 (k) 2 (n x <5 3 (k f ) + 1) (5.29) 

Note that the inflaton condensate populates the k = momentum mode so the 
occupation number for this mode is n x S 3 (k^) where n x is the inflaton condensates 
number density. For the £ spectrum, n^(k) is a strongly peaked function at |k| = 
but not as narrow as the delta function. The net rate of decrease of inflaton condensate 
particles is proportional to the difference of the above two matrix elements which is 
proportional to n x (l + 2n k ), where we have applied n x 5 3 (k = 0) 3> 1. Defining the 
occupation number of the £ particles as n^(k) = n^(k) = n^(— k), the decay rate is 
approximately, 

r "" , ftfc( 1 + 2n < (* = t)) (5 ' 30) 

where the prefactor is necessary to be consistent with our previous result for T x . To 
determine the number density we must integrate n^(k) over all momenta. If we 
can determine the approximate width of the function n^(k), we should then be able 
to approximate this integration. To do this, express the potential as follows, 

- \ m \ - 9Xi 2 = -g [ m \ + 2 9Xo sm(m x t)] £ 2 (5.31) 

where the solution to the equation of motion for x nas been used, (x) = Xo s ^ n ( m x^)- 
We have also assumed the scale factor is kept constant and neglected the depletion 
of the condensate. In the above x is treated classically, and the potential is written as 
the effective mass for the £ field. The momenta of the produced £ particles are then, 



k = V (?) 2 ~ m " ~ 29X0 sin ( m ^) ( 5 - 32 ) 
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Additionally, if both gxo and m| are small compared to m 2 ,, then the particles are 
produced in narrow band of phase space, centered at k 



f with width 



Ak 



mi 



(5.33) 



We can then integrate n^{k) over momenta to approximate the number density 

mgx Q n k 



ATxkln k ^k 



(27T) 



2tt 2 
gn x n k 

2tt 2 Xo 



(5.34) 
(5.35) 



where the virial theorem has been applied in the last step. Using this last relation, the 



decay rate (5.30) may be written 

-L v 



1 + 



(5.36) 



Substituting this into the rate equation (5.27), and assuming n k > 1, one obtains, 



7i£ oc exp 



n 2 9Xo 
m 2 



N 



(5.37) 



where 



m x t 
2vr 



is the number of oscillations of the inflaton field. The strength of 



the effect can be read from the exponent of (5.37) by determining the number of 
oscillations of the inflaton at which the depletion of the inflaton condensate becomes 
appreciable. Note that when approximating the resonance, we assumed a small cou- 
pling g < so even in the case of small coupling the exponential growth is obtained. 
The above example of narrow resonance demonstrates that the perturbative results are 
not always adequate. We should note however, that the narrow resonance is sensitive 
to the expansion of the universe. When this is included, the effect is diminished as 
the occupation numbers of any particular momentum mode are constantly shifting to 
lower momenta and reducing the occupation number at k = m x /2. Additionally, one 
should take into account the depletion of the inflaton condensate and interactions in 
which produced £ particles back-scatter on the inflaton condensate. Inclusion of such 
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effects, which are also necessary for energy conservation, can also change the final 
result. 

Note that when the coupling is large, g > — , our above result is invalid as the 
resonance is no longer narrow. In this case, the exponential production may still be 
obtained but the resonance is now for a broad range of momenta. We must use the 
methods of Section [3] which are generally applicable. For the case of large coupling 
consider instead the quartic potential \gx 2 ^ 2 which is similar to the flat direction 
potentials later considered. Here, the coupling g is dimensionless, and the regime of 
strong coupling is g > —. The perturbative decay corresponds to the annihilation 
process XX ~^ ££■ We perform the nonperturbative calculation for the production of 
£ quanta treating x as a classically evolving field (x) = Xosin(m x t). Again, we are 
freezing the scale factor, and neglecting the depletion of the inflaton condensate. The 
effective potential for the £ field in this case is, 

-\{ml + X 2 o^ 2 (m x t))e (5.38) 
and the equation of motion for the momentum modes £ fc are, 

4 + uo(t) 2 £ k = , oo(t) 2 = k 2 + m| + Xo sin 2 (m x t) (5.39) 

which is of the form discussed in Section |3j so the methods Section [3] can thus be 
applied. We will not work through the details of the calculation which may be done 
either numerically or analytically (see [I5])- We know that £ particle production may 
occur when the adiabatic parameter ^ is greater than one. The adiabatic parameter 
is computed, 

\u\ = m x ~g 2 xl\ sin(ro x t) cos(m x t)| 
uj 2 [k 2 + m 2 + g 2 x\ sin 2 (m x t)) 3 / 2 

and notice that for very small masses and momenta, the adiabatic parameter is 

greater than one for short durations of time when sin 2 (m x t) <C 1, or expanding about 

a time m x t* = nn in the small time interval e rn' 1 , the adiabatic parameter is 

approximated, 



\u\ _ m 2 x g 2 xle 



J 2 [k 2 + g 2 Xo'm x € 2 ) 3 / 2 
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(5.41) 



The upper bound of the resonance band is then the maximum value of k 2 such that 
^4 > 1 is still true, or when 

k 2 < [m x ~g 2 xl(rn x e)} 2/3 - g 2 x 2 Q (m x e) 2 (5.42) 
< m x gxo (5.43) 

where the right-hand-side of the above expression has been maximized with respect to 
variation of e. We note immediately, that momenta greater than m x can be generated 
which implies that many \ particles are involved in the decay process. Specifically, if 
we take Xo = M p and m x = 10~ 6 M p , then the momenta of the produced quanta may 
be as large as k ~ 10 3 m x . Additionally, when the occupation numbers are computed 
using the methods of Section |3j one finds an exponential growth of particles during 
these short time intervals, and a quick depletion of the inflaton condensate (within 
25 oscillations) for a range of couplings 1CT 3 > g > 1CT 6 and with m x = 10~ 6 M p 
(results quoted from [16J). 



5.7 Nonperturbative Decay of the Flat Direc- 
tion VEVs: Preliminaries 

This section provides the motivation for our choice of flat direction models in Section [6] 
and is based on the arguments of [2S1I2Z]. The focus is on the nonperturbative decay 
mechanism and establishing criterion for an efficient nonperturbative decay. We begin 
by noting the above model of inflaton decay via coupling gx 2 ^ 2 ,s n °t appropriate to 
model flat direction decay as the flat directions are complex fields. The analogous 
coupling would be 

g\$\ 2 e (5.44) 

where $ is the flat direction VEV, and £ is a real light scalar field. The motion of 
the flat direction is typically an elliptical motion in the complex plane as shown in 



Figure 4.3 The induced mass ~ |$| will never pass through zero and thus will 



never result in a large adiabatic parameter (5.40) as happens for the inflaton. To 
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shown this, one may estimate the adiabatic parameter shortly after the onset of flat 
direction linear oscillations to be, 



< ^susy 

4 



IT i ( 5 - 45 ) 

Lot e B _ L |<I>o| 



where cu^ ~ ~ |$ | ar| d where £b-l is again the B — L number per particle in 
the flat direction condensate defined in (4.40). Taking $ ~ 10~ 2 M P One then re- 
quires £b-l ~ ~ 10~ 14 in order to obtain an order one adiabatic parameter. One 
might suggest that the early nonlinear evolution of $ prior to Rn~m may contribute 
to non-adiabatic evolution, but this evolution would be model dependent and also 



conditional on the higher order terms of the flat direction potential (4.28) dominating 
over the mass term (see Section [4]). Our goal is to look for a model independent 
decay mechanism so we focus on the behavior during the linear oscillations of the flat 
direction when R > Rh^™- 



There is another problem with the coupling (5.44), and it is that in the context of 
the MSSM with all gauge symmetries broken, the field £ should correspond to either the 
heavy Higgs-like scalars or to the orthogonal light states which are flat direction quanta 



5(f) a . For couplings to the heavy scalars, one naively expects the conclusion (5.45) to 
a pplyj^] For the flat direction perturbations, this conclusion does not necessarily 
apply. Additionally, as has been discussed in Section [2~4| one rather expects a situation 



in which multiple flat directions have acquired VEVs so the effective mass terms (5.44), 
should properly be represented as a matrix expression composed of the heavy and light 
states. The matrix elements generally acquire time dependent terms of the sort $ a . 
Because the time dependence is on the same scale m susy as the light eigenstates of the 
mass matrix, one expects nonperturbative particle production. Then through analysis 
and numerics, as discussed in Section [3J one may determine whether this is in fact 
the case. We note that the production in the multifield case is sourced by both the 
rapidly changing eigenvalues and rapidly changing eigenvectors (see Section [3j We 
also note, that for reasons discussed in Section 3A the heavy Higgs-like states of the 
theory should not immediately be discarded as they will generically mix with the light 
9 Although we will see these heavy states may in fact be produced 
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states and may be produced on shell. Thus the mass matrix above M 2 should be 
expanded to include these heavy states. Specific models which further illustrate the 
above sketch will be presented in detail in Section [6] which follows. 

Before concluding, we note some potential obstacles to the nonperturbative decay 
of the flat directions. The first potential obstacle is the kinematic availability of states 
for the flat direction to decay to. Specifically, the flat direction quanta have mass 
m susy , and there may be limited or zero phase space for decay of flat direction quanta 
into other states of the same mass if there are equal numbers in the initial and final 
state. For instance, if the only allowed reactions were 2^2 reactions such as an 
s-channel process 0^0i — > 0202, the resonance band would be limited by the difference 
of the masses k 2 ~ m\ — m\. However, if the process is nonperturbative there may 
be more particles in the initial state than in the final state {(p\(pi) N — > 0202 which 
broadens the resonance by the ratio N . As we saw at the conclusion of the previous 
section for the case of the inflaton, the ratio could be as large as a thousand. This 
can make kinematics less of an obstacle. However, the flat direction condensate may 
possess a certain amount of B — L charge at the start of linear oscillations, and one 
must confirm that B — L conservation will not constrain the decay channels of the flat 
direction. For instance in the annihilation example just mentioned, ((fi{(j)i) N — > 0202, 
if the initial state possesses nonzero B — L, then the final charge must match. This 
could restrict the allowed interactions to iV — > N, and potentially restrict the decay 
rate if mi ~ m 2 . Charge conservation issues such as this are raised for instance in [37J. 
However, a supression of the flat direction decay as envisaged in [37] relies heavily on 
the condensate being maximally charged, and as we will see in the Results, a sizeable 
charge asymmetry is still consistent with nonperturbative decay. 

A second potential obstacle concerns the heavy states should they happen to be 
produced in our calculations. By dimensional analysis, the subsequent decay rate of 
our heavy states should be, 

T heavy ~ aM (5.46) 
where a is a gauge coupling and M is the heavy scale. Since the gauge couplings 



are all less than one (see Figure 5.2) at the TeV scale and above, this suggests that 
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^heavy < M. This conclusion is also supported by the analogy to the standard models 
Electroweak symmetry breaking and the decay width of the heavy W and Z bosons 
or the standard model Higgs. These states are analogous to the heavy states in our 
Lagrangians, and their widths are measured (or expected in the case of the Higgs) to 
be less than their mass [76J. If this inequality is not obtained for our heavy states, 
meaning that the decay width of our heavy states is comparable to or larger than the 
mass of the states, then these states are virtual. This would make the computation of 
their occupation number nonsensical as the notion of a particle state is not available. 

From the above considerations, our expectations are that once we have chosen 
a realistic flat direction model, one may obtain a resonance band for momenta k < 
m susy . One must verify that interactions exist within the model to mediate the decay, 
and one must verify that the classical dynamics of the VEVs can lead to an adiabaticity 
matrix with elements greater than order one. The only way to verify these two points is 
to specify a model, solve the classical evolution, and apply the methods of Section [3j 
These issues are explored in the following section for some simple models of flat 
directions. 



5.8 Summary of Reheating 

The above discussion of reheating is only a sampling of the scenarios for the reheating 



process considering very generic decay mechanisms (see Figure 5.1). The full range 
of scenarios is vast and corresponds closely to the vast array of inflaton models |64|, 
the large uncertainty in the energy scale of inflation, and the unknown appropriate 
particle/field description during reheating. These are major open questions. It is 
worth noting that answers to the latter two questions may come in the not too distant 
future as observational cosmologists are seeking to measure gravitational waves (and 
thus the energy scale of inflation) and particle physicists apply the LHC to characterize 
the TeV scale. 

In the remaining material, we aim to answer questions related to the effects of 
the flat directions on a situation of gravitational decay of the inflaton. We saw that 
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the thermalization process can be potentially delayed by the presence of flat direction 
VEVs, but more importantly, we saw that the flat direction can in fact dominate the 
energy density of the universe if it does not decohere or decay quickly. The following 
sections are focused on the realistic modelling and simulation of flat direction dynamics 
with the goal of estimating the decay time of the flat directions, and thus resolving 
these questions. 
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Chapter 6 



Modeling Flat Direction Decay 



The previous four sections aimed to present the broad picture of inflation and reheating 
in the context of MSSM flat directions, and to present technical tools for studying the 
evolution of scalar fields in the early universe. In this section, the picture necessarily 
narrows in order to make a concrete calculation. Specifically, the full SU(3) xSU(2) x 
£7(1) gauge symmetry of the standard model is here truncated to its U(l) part, 
and the full spectrum of scalar fields in the MSSM is truncated to a case of two 
complex scalars and a case of four complex scalars. The two cases correspond to one 
and multiple flat directions respectively. To describe the effective vacuum in these 
models, the fundamental scalar fields are split into a classical background (which 
are the flat direction VEVs) plus quantized fluctuations away from this background. 
It is then necessary to fix a gauge, and we fix both the background fields and the 
fluctuations to the unitary gauge. The backgrounds will be chosen to possess only 
time dependence and the perturbations to the background will possess the full space- 
time dependence. However, by evolving the perturbations one may determine how 
quickly spatial fluctuations develop (flat direction decay). The goal of this section is 
thus to express the two Lagrangians of our models in a canonical form that may be 
solved numerically in Section [7} 
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6.1 Gauge Fixing U(l) Symmetric Models 



The actions we consider consist of N complex scalars (pi which are similar to be the 
scalar particles of the MSSM with the main difference being the gauge group here 
contains only the U(l) part. We develop two simple U(l) models which have the 
generic form, 

S = j d 4 a;^|-^(S) 2 + ^|^| 2 -F(0i,0 2) ..A)J (6-1) 

where is the U(l) gauge field, qi are the charge assignments of the scalars and 
= (dfj, — ieqiB^) (pi is the covariant derivative. The fields thus transform under 
the gauge symmetry as follows, 

Bp - B, + ^ 

e 

& - 

We fix the fields to the unitary gauge where the degrees of freedom of the fields match 
the physically propagating degrees of freedom^ This gauge choice provides the 
most straightforward and intuitive field description in our canonical quantization. A 
simple way to fix to the unitary gauge is to combine the fields in the Lagrangian into 
gauge invariant combinations (for example see [16J). To do this, we first decompose 
the complex scalars into the representation 0, = ^e tai with fi > and with the 
scale factor appearing in the denominator. Of course, this decomposition will lead 
to a coordinate singularity if any one of the fi ever becomes zero. However, our 
(pi acquire nonzero vacuum expectation values, and by assumption will not cross the 
origin except in the specific limiting case where any of the phases o { are constants 
of the evolution. The decomposition is thus acceptable. To economize notation we 
will also write many expressions in this section setting the scale factor to a constant, 
R — 1. The time dependence of the scale factor will be restored in the final results 

1 Equivalently, there is no extra gauge-fixing constraint to enforce on the classical equations of 
motion or Gupta-Bleuler constraint to enforce on the quantized fluctuations. 
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using the prescription derived in Appendix |B| Using this representation, 4>i — jj e *' 
the gauge transformation is, 



The gauge fixing procedure is to construct linear combinations of the phases a, ; among 
themselves and with which are individually gauge invariant. The unique gauge 
invariant combination for the vector is, 



B^ = B^ — , a = ' 1 2 * which transforms a — > a + A (6.2) 



It is quick to verify that a, which will be the goldstone mode, transforms under the 
gauge symmetry as specified, and also quick to verify that the vector is gauge 
invariant. One may then combine the phases a,i into a basis of gauge invariant combi- 
nations 9i. This will be done concretely for the two examples below. We simply note 
here that our basis of gauge invariant phases 0j will be N — 1 dimensional and unique 
up to a time independent 0(N — 1) rotation. We note finally that the above gauge 
fixing procedure may be equivalently interpreted as performing a specific gauge trans- 
formation on the fields in which the combination of phases a is transformed explicitly 
to zero. 

With the decomposition = fie lai , the kinetic terms of the scalars may be 
written, 



W,? = E + f^d^f] - B,r + -M 2 (f t )B,B» (6.3) 



where the gauge current J /t and the effective mass of the vector M 2 are defined, 




fi 



fi 



9 




(6.4) 
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In the unitary gauge, the gauge kinetic terms (6.3 6.4) retain their form except for the 
unphysical phase a being fixed to zero. Similarly, the gauge current term J M is now 
a "physical" quantity and this offers some computational convenience. Specifically, 
because our fields are assumed only time dependent, then gauge charge conservation 
requires d {J ) = or equivalently ( J ) = J™ lt . Here, J™ 11 is the initial gauge charge 
density, and it is assumed to be exactly zero. This is imposed for phenomenological 
reasons, although from a mathematical stand-point, there is no obstacle to assuming 
a uniform charge density for the background. In the following material it will be 
apparent that gauge charge conservation has nontrivial implications to the dynamics. 



One final aspect to the gauge charge (6.4) to be noted is that when all phases are all 
frozen, d^ai, the charge is automatically zero. In our later numerical simulations, the 
initial charge will be made zero by this method. 



6.2 U(l) Model With a Single Flat Direction 

We now perform the above gauge fixing for the simple U(l) example introduced in 



Section |2.2| consisting of two oppositely charged scalar fields <\>\ and <p2, and one 
complex degree of freedom in the flat direction VEV. The potential for the model 
includes mass terms and the D-term, 

Vfa, 02) = m?|0i| 2 + m^l 2 + \<?{\<t>i? - |0 2 | 2 ) 2 (6.5) 

where it is assumed both masses are of order the SUSY breaking scale, m susy . One 
also reads from the D-term that <pi = f\e l<Tl and 2 = /^e" 72 are charged +1/2 and 
— 1/2 respectively. The phase a = (o- 1 — a 2 ) is then the gauge variant goldstone mode 
while the combination (01 + 02) is a physical gauge invariant phase. We transform the 
vector as specified above in (6.2) and we write the Lagrangian in terms of the gauge 
invariant "physical" fields B^, fx, f 2 , and (01 +02), 

{d,hf + {d,f 2 f + (^^) [9>i + V2)} 2 -V(f u h) 

~\f (b) 2 + ie 2 (J 2 + f 2 ) - \eW [f\ - ff) d^a, + o 2 ) (6.6) 
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which is the unitary gauge action in which the goldstone mode a has necessarily 
dropped out. The effective vector mass as well as the gauge current are easily read 
from the above Lagrangian, 

M 2 = \e\fl + f 2 ) , J» = le(fi-fl)d»(a 1 + a 2 ) (6.7) 

We make the following demands on the classical solutions: (i) they are only time 
dependent, (ii) they result in ( J ) = and (iii) the phase (ax + cr 2 ) is dynamic so that 
the fields may rotate in the complex plane (see Figure 4.3). These assumptions 



and the above form of the gauge charge density (6.7) restrict the solutions to the case 
(fi) = (f 2 ) = F which puts the fields precisely on the flat direction. However, one can 
show that the classical equations of motion will force F = given our assumptions. 
To show this, one writes the equations for (/x) and (f 2 ), and takes the difference of 
the two equations to obtain, 

[m\-m\)F = (6.8) 

Thus, in order to obtain a solution with nonzero (fi), we must relax one of our 
assumptions. A simple course is to take the special case m x = m 2 , but before doing 
this, we consider alternatives. We must retain a nonzero VEV F as well as a dynamic 
phase (ctx + ^2) as these define precisely the case we wish to study. It is also difficult to 
explain a case with a uniform nonzero gauge charge density from a phenomenological 
standpoint, so we do not relax this assumption^] We could relax our assumption of 
only time dependence on the fields - it may be that our model becomes consistent only 
when classical spatial fluctuations to the fields (which would be small after inflation 
but present) are included and local current conservation <9 M J M = is enforced. This 
case seems interesting, and it implies a classical instability (and thus decay) of the 
flat direction, but it is not explored herej^] We take the special case mi = m 2 , 
and study this theory because it is well defined, not completely trivial and also helpful 



2 Assuming the charge is one of the standard model gauge charges, it would have to be extremely 
small or its coupling extremely weak in order to have not been observed. 

3 One final possibility is that our model is simply "sick" as a physical theory due to our assump- 
tions on the gauge symmetry and the field content having been too narrow. However, one should 
not conclude this until the other cases have been explored. 
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for understanding the multi-flat direction case considered in the following subsection. 
Although there is a fine tuning involved, one expects that the model will still be 
predictive in the case of slightly different masses for a time duration of approximately 
— - — which is the time required for the field oscillations to acquire a significant 
phase difference. As there are arguments for some of the MSSM scalar fields being 



(B^) = B^ , (a 1 + a 2 ) = S , (f 1 ) = -F , (f 2 ) — -F (6.9) 



nearly degenerate in mass (see Section 2.5), this model should not be ruled out on 
phenomenological grounds. 

We continue to determine the classical background for the case m\ = m 2 . The 
vacuum expectation values of the fields are then specified, 

\f , CA>4 

With regard to the vector field, because the background vector is assumed to be only 
time dependent, one can show the non-dynamical component satisfies £>o = = 
(see also Appendix [b|. Additionally, at the end of inflation it is assumed B = 0. 
The spatial component of the vector B can then develop a large nonzero value only 
from the effect of parametric resonance. Since M! C M, it is not immediately 
clear that a parametric resonance is possible. We take the background solution to 
be B = 0, and the presence of a parametric resonance for the vector will be checked 
by studying the Lagrangian for the perturbations in the following discussion. The 
background equations of motion then involve only the scalar fields F and E from the 
above Lagrangian, and have the form, 

F" + (m 2 R 2 - |S /2 - f) F = 
\d [F 2 £'] = 

f +{f) 2 = AnGR 2 (2V^)+p x ) (6.10) 

where the dependence on the scale factor has been restored and the energy density 
of the inflaton p x has been incorporated. Note that when written in the complex 
representation, our VEVs are defined = (0 2 ) = $ = ||e ,E and using physical 
time, the above equations of motion are simply the form (4.35) of Section 4.4 where 
the solutions to the equations have also been discussed. 
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The decay of the flat direction VEV may be quantified by studying the perturba- 
tions to the background, and these are parametrized as follows, 



( C r 1 - cr 2 ) = S + ° 



f 



F 

F + f + 6 
2 



h - ^ (6.11) 

Note that while the background quantities have only time dependence, the pertur- 
bations f, a, and S possess the full space-time dependence. Substituting this 
decomposition into the Lagrangian one obtains the quadratic Lagrangian for the per- 
turbations, 

\ [(dj) 2 + (<V) 2 + («V) 2 ] + Ijy (fa' - fa) 

-Um 2 - -V 2 ) f 2 --(m 2 - \i:' 2 ) a 2 --( -e 2 F 2 + m 2 - ^£' 2 ^ 6 2 
2 V 4 J J 2 V 4 / 2\4 4 J 

-\[F,u{b)} 2 + \ (^*pj - b (6.12) 

Let us examine the couplings above. Notice the vector perturbation b has already 



decoupled at quadratic order since the current (6.7) does not contain a part linear in 
the perturbations. The elimination of the non-dynamical component of the vector 6 
contributes a mass term for 5 but otherwise does not effect the quadratic Lagrangian 
(see Appendix[B]). The excitation S is the Heavy Higgs which to leading order has mass 
Ai 2 and also decouples. The perturbations / and a are apparently mixed, but these are 
simply the flat direction perturbations and they are degenerate with mass m to leading 
orderj^] The degeneracy of the eigenvalues means the corresponding I, J matrix 
elements are zero to leading order and thus will not lead to a parametric resonance. To 
conclude, the eigenstates and eigenvalues in this model are stationary and the model 



4 



To quickly see that / and a are the flat direction perturbations, look at the definition of the 



perturbations (6.11), and see that they correspond to a shift on the VEV F and X. 
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is incapable of leading to parametric resonance. This conclusion is expected as there 
are no states in the above model lighter than m for the flat direction quanta to decay 
into. This shortcoming is corrected in the multi-flat direction model presented nextj^] 



6.3 U(l) Model With Multiple Flat Directions 

One suspects that introduction of more degrees of freedom to our U(l) model may 
change the outcome. We thus perform the same analysis on the slightly more com- 
plicated model involving four complex scalar fields specified by the potential, 

V(0 1; 2 ,0 3 ,04) = m 2 (|0 1 | 2 + |0 2 | 2 )+m 2 (|0 3 r + |04r) 

+ ^e 2 (|0 1 | 2 -|0 2 | 2 + |0 3 | 2 -|04| 2 ) 2 (6.13) 

where the charge assignments can be read from the D-term and where the mass 
choices will be discussed shortly. Again the fields are decomposed fa = fit mi . The 
physical phases and the goldstone mode are, 

01 = (<7l - 2 - a 3 + <7 4 )/2 , 6 2 = ((?! + <7 2 )/2 6 3 = ((7 3 + <x 4 )/2 

a = (cri - cr 2 + 0-3 - cr 4 )/2 
The Lagrangian in the unitary gauge is then, 

E(^) 2 + i(A 2 + fl + fl + f!)(dA) 2 + (A 2 + /2XW 2 
i=i 

+ (fl + fi)(d»e 3 ) 2 + (A 2 - ^)(W(W - (fl - fi)(9A)(9A) 

- V (fi, h) - \F (B) 2 + ie 2 (fl + f 2 + f 3 + fl) B,W 



- eB^ 



\(fl + fl - fl - fl)(dA) + (fl - fl)(d,o 2 ) + (fl - tf) WW 



(6.14) 



5 In principle a decay may occur in the single flat direction model with mi ^ m^, though 
as discussed above there are subtleties to this case. These subtleties are mentioned again in the 
Conclusions. 



120 



where the gauge current J M is the coefficient of the term linear in and it is equiv- 
alent^ written in terms of the Oi fields as, 

J„ = e [fld^ai - fid ^a 2 + / 3 2 <9 M a 3 - fid^t] (6.15) 

where, 

(Ti = e 2 + ^9x , a 2 = e 2 - l -e 1 , a3 = 9 3 -^9x , ^ = 9 3 + ^9x (6.16) 

Next we determine the background. The background for the vector is again zero. 
This leaves seven scalar fields fx, f 2 , f 3 , f±, 9 1} 9 2 , 9 3 participating in the flat direction 
background dynamics. We wish to select solutions to the corresponding equations of 
motion for which the fields will satisfy the gauge currrent conservation (J ) = 0. By 



examination of the current term in (6.14), one finds a convenient solution in which 



(9i) = const., fx = f 2 and / 3 = / 4 . This solution is convenient because it evolves 
exactly on the flat direction as the previous model did. The equations of motion are 



thus two copies of the system (6.10 ) each copy evolving independently except for weak 



coupling through the Friedmann equations. The VEV's are defined as follows 

(fi) = {f 2) = F , (9 2 ) = E 
(f 3 ) = (f A ) = G , (9 3 ) = t 
(9x) = 

or in the complex representation, this is written 

<£ 

71 



(</>!> = <0 2 >^=Fe ffi (6.17) 
(0 3 ) = (0 4 ) = 4= = Ge iS (6.18) 
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where $ and $ have been defined with a factor v2 to make the Lagrangian for these 
complex VEVs canonical. The background and perturbations are then specified, 



fl = F + (f + 5)/2 
f 2 = F + (J-S)/2 
f 3 = G + (g + 5)/2 
f 4 = G+(g-~5)/2 



0i 

03 

B„ 



9 E 

T - 



Vf 2 +g 2 

2FG 

a 

IF 
a 

2G 



C + O{^ 



(6.19) 



where 6 M is the vector field perturbation, and where the definitions have been con- 
structed again to make the fields canonical. The elimination of the time-like compo- 
nent of the vector b is performed again in Appendix[B| The flat direction perturbations 
{/, a} and {g, a} decouple from eachother and from the remaining fields, and they 
are degenerate with masses m and m respectively. The transverse and longitudinal 
vector degrees of freedom decouple to leading order in the expansion The remain- 
ing three degrees of freedom 5,5 and 9 C are strongly mixed. Specifically, forming the 
vector X = (5, 5, 9 C ), the quadratic Lagrangian is written in the canonical form (3.28) 
with matrices (ti 2 + K T K) and K which are 



(n 2 + K T K) n 
(Q 2 + K T K) 

(n 2 + k t k) 
(q 2 + k t k) 12 
(n 2 + k t k) 



22 



k 2 + m 2 R 2 - 
k 2 + rh2 R 2_ 



IT 
~R 

~R 



+ 



+ 



e 2 + 



e 2 + 



3E 



12 



F 2 + G 2 

3E' 2 
F 2 + G 2 



G 2 



k 2 + 



FG 



(m 2 G 2 + m 2 F 2 )R 
F 2 + G 2 

3E' E' 



R" 3 (FG' — G F') 
~R + <™ < ^ 2 



(F 2 + G 2 Y 



e 2 + 



13 



F 2 + G 2 
3E' (FG' - G F') 



(n 2 + k t k) 23 



K 



12 



K 



13 



VWTg 2 

G 

VF 2 + G 2 
-GT! 
~~ VF 2 + G 2 



F 2 + G 2 

3E' (FG' - GF') 
F 2 + G 2 



K 



23 



FT! 
VF 2 + G 2 



(6.20) 
(6.21) 
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and with uj 2 and r,-„ which are 



uj 3 

r i2 

Tl3 
T23 



+ 



R 2 (m 2 |$| 2 + m 2 |$| 2 ) 3(lm ||$| 2 + ||$| 2 



~R 



|$|2 + |$|2 (|$| 2 + |$| 2 ) 2 

R 2 fm 2 |$| 2 + m 2 |l>| 2 ' 



+ 



3|$|2|$|2 



(|$| 2 + |$| 



2^2 



+ 



R 2 (m 2 |$| 2 + m 2 |$| 2 
|$| 2 + |<3>| 2 







m 
M 2 



|*||*| 



|$|2 + |$|2 



R 2 ( 



lm 
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lm 
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I M 2 
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(6.22) 



Where the expressions for uji and Tj have been written in the complex representation 
of the fields. The mixing between the three fields is a strong mixing that can not be 
removed by any small rotation. The system is diagonalized by an 0(1) time dependent 
rotation, and is expected to lead to parametric resonance. From the eigenvalues, one 
sees the system includes a heavy Higgs excitation as well as two light excitations^ 



These results (6.22) will be applied in the following section to verify the parametric 
resonance and to compute the decay of this flat direction model. 



6 The light degrees of freedom are the third flat direction of this model. The gauge current 
conservation on the background and our mass assignments however prevented this flat direction from 
obtaining a VEV 
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Chapter 7 



Results 



Here we perform the numerical implementation of the four field U(l) model presented 
in section 6.3 The model is converted into a form appropriate for numerical solution, 
the initial conditions of the system are discussed, and some useful measures of the 
decay are also presented. Using arguments based on the scaling of the Heisenberg 
equations (or equivalently equipartition of energy), it is then shown that of the four 
parameters m, m, |$ |, |$ | that could be varied, only the ratios fh/m and |$ |/|$ | 
need be varied with the remaining cases being obtained through rescaling. The pro- 
duction rate of quanta is reported for a few values of fh/m and for a continuous range 
of values |$o|/|^o| i n which the particle production is exponential. Additionally, the 
decay time of the toy flat direction model is determined in this resonance band with 
all data having been rescaled to phenomenologically relevant TeV scale. 

7.1 Numerical Setup of the U(l) Model With 
Multiple Flat Directions 

We study the nonperturbative decay of the multi-flat direction model presented in 



section 



6.3 The background equations are of the form specified in (6.10) for the 



two separately evolving flat directions. The Heisenberg Equations of motion (3.21) 
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will be solved with the driving terms determined in (6.22). It is necessary to give the 



flat directions initial charges, and rather than impose the charge by hand, the charge 
is introduced dynamically according to the susy breaking potential terms motivated 



in (2.39). The effect on the solutions is to give the phase of the flat direction an 



initial kick so as to induce a rotation in the complex plane (see Figure 4.3). These 
quartic terms become sub-dominant to the mass term after a short period of time, 



so their modification to the eigenvalues and Gamma matrix elements (6.22) may be 



neglected after this time. We additionally introduce a matter fluid, to represent the 
inflaton energy density which dominates throughout. The background fields defined 



in (6.18) are again, 



(<h) 

(03) 



(02) 
(04) 



$ F 
<£ _ G 



(7.1) 
(7.2) 



where the scale factor has been restored. The scale M = 10 _2 Mp is introduced 
which represents the approximate scale of the VEVs, and the following dimensionless 
quantities which will be used for numerics are then defined, 



V* 
m* ee 

F, EE 



eMrj 
rn 



eM 
F 



m* ee 



m 
eM 



M 



G. t . 



G 

M 



Px* 



k 
eM 
Px 



A* 



A 



A* = 



A 



e 2 M 2 M p 2 



(7.3) 



With these definitions, the coupling e becomes an overall constant in the Lagrangian 
which will not enter the equations of motion. Also the scale M will now appear only 
in the equation for the scale factor. Our choice of quartic couplings 



A 



m 



10|$ | 



A 



~ 2 

m 



10|$ | 



(7.4) 



which are of the form described in (2.39). The above the couplings A and A are 



assumed real with complex phases having been absorbed by a redefinition of the fields 
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(and thus appearing in the initial values of the phases). The additional factor of ^ 
above was chosen to be as large as possible without introducing instability in the 
background equations. Using our dimensionless representation, the complete set of 
background equations are, 

R 



f'; + 
g: + 



mlR 2 



A* 



F„ 4 sin (4E) = 



JDII 

~ 2 D 2 lh w2 
m *R — - E 

R 



Gz E' ) -^Gfsinf I E 



G» + y G» cos f 1 E 



M 2 
M 2 



™* F 2 + 



Tip" cos (4E) 



+ m 2 G 2 + T^cos(4E 



+ it> 2 p x ^ (7.5) 



where prime now denotes derivatives with respect to 77* , and where the equation for 
the scale factor was expressed in a way so that the time derivatives of the fields $ and 
$ do not appear. In the above background equations, the evolution of the inflaton 
energy density will simply be that of a dominating matter fluid p x * = p x *(0) (^f^ ■ 
In the dimensionless units, the initial conditions of the background equations are, 



R(0) = 1 , p x *(0) = l (7.6) 
F*(0) = F , E(0)=0.25 , G,(0) = G , S(0) = 0.25 (7.7) 

F*e* E V _ (G^ 



R 







R 







(7.8) 



where the initial condition on the scale factor is chosen R = 1 for convenience. The 
initial condition for R' is then set by the Friedmann Equation. The phases E and E 
must be nonzero in order for an initial charge to develop through the A and A terms. 
Also in the above, we have chosen p x * large enough initially so that H ^> m,rh and 
thus it is consistent to set the fields $ and $ as frozen, thus fixing the initial conditions 
forF;,G;,E',E'. 
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The Heisenberg Equations in the dimensionless computer units are, 



a' 



Uu;* 



h)(3 + 



2u, 



(7.9) 



where we have introduced the dimensionless diagonal frequency matrix and dimen- 
sionless T, /, J matrices 



to 

~eM 
1 

2 



eM 

1 1 

± 



(7.10) 
(7.11) 



Recall that u 



diag(o; 1 , u>2, u> a ). The reader is also reminded of the adiabaticity 



matrix elements defined in (3.32) which again are 
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UJ CO 



A 



l = j 



no summation 



(7.12) 



where one notices these quantities are explicitly dimensionless and invariant under the 



redefinitions (7.3) 



The initial conditions for the Bogolyubov coefficients again are a = 1, /5 = 0. 
The evolution of the Heisenberg Equations is begun at the time when the smallest 



eigenvalue of (6.22) is positive and equal to twice the R" / R term of these eigenvalues. 
This ensures that any significant production of quanta will be due to the effects of 
the evolving flat direction VEVs and not due to the gravitational background which 
appears through the term R" /R. 

The occupations numbers are computed with rii = {(3* (3 T )u, and these along with 
the corresponding eigenfrequencies are used to compute the energy density of the 
produced quanta. In particular, the ratio of energy density of the produced quanta to 



127 



the energy density of the background flat directions is, 



prod 



= fH>**\ f dK k l uUK)ni{K) \ (7 .i3) 
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+ R 2 F* 
+ R 2 Gl 



m * + T§ cos(4S) 



^ 2 + y§cos(4S) 



This ratio is zero at the start of the evolution, and we evolve the equations until the 
ratio becomes unity at which point, we consider the flat directions to have decayed. 
Once r pro d = 1, the present analysis based on the quadratic Lagrangian is also invalid 
since we are ignoring the back-reaction of the produced quanta on the evolution of 
the flat directions. For example, we do not account for the decrease in the amplitudes 
of the flat directions due to the particle production. However the instant r pro d is a 
useful measure of the decay of the flat direction. As a measure of the time for the 
decay, we use the number of rotations of the first flat direction, N = (S — S ) /2n 
as this is somewhat easier to visualize than the physical (or conformal) time and in 
particular it is independent of the scale m. We thus denote by N^ecay the value of N 
at which the ratio r proc i equals to one. 



For numerical solution, the background equations (7.5) which are five second order 
equations are easily rewritten as a ten dimensional system of the form X' = f(X) 
by defining the state vector X = (F, F', E, a', G, G', S, R, R'). The numerical 
solutions are easily obtained using a typical desktop computer (1 Gigahertz Processor, 



1 Gigabyte memory) in approximately a minute or less. The Heisenberg equations (7.9) 
may also be written in this generic form with a state vector of length 36 composed 
of the two 3x3 complex matrices a, (3. However, because these equations require 
the background solutions for the computation of the driving terms, it is convenient 
to solve both the background equations and the Heisenberg equations as a single 
system in a vector X of length 46. Numeric solution of this system is demanding 
due to the increased number of dependent variables, and also due to the periods of 
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non-adiabatic evolution resulting in a numerically stiff differential equation for these 
periods. The solutions for a single value of momentum and for approximately five 
to twenty complete rotations of the flat directions can require up to few hours of 
computer time using a standard differential equation solver such as ODEPACK. 

There is one further aspect to the numerical solution which requires attention 
and this concerns the series expansion in the small parameter || which was used in 



obtaining the driving terms (6.22). The series expansion is in fact more properly an 



expansion based on the following hierarchy, 
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(7.14) 



This hierarchy is actually more restrictive than the condition S < 1 because the non 



adiabatic evolution can make terms such as ? sometimes larger than the masses 
m,rh. However, if m,rh are chosen small enough, the above hierarchy and the series 
expansion are valid for the whole numerical evolution. In the numerics, this has been 
arranged to be the case. 

To summarize, once solutions to the background and Heisenberg equations have 
been obtained through the above setup, the scaling arguments (discussed next) may 
be safely applied, and the numeric results obtained. 



7.2 Scaling of the Numeric Solutions 



In this section the scaling arguments of Section 3.4 are applied to show that of the 



l*o] 

the numerical simulations. The solutions of the remaining cases may be obtained 



four parameters $ , $ , m, and in, only the ratios ™ and L 5 ! nee d be varied in 



through the rescaling shown in Table 3.1 The heavy scale in our model is the scale 



of the VEVs and specifically it is taken to be M = y $0^0- The light scale is the 
scale of the masses and we use m for this scale. Recall the scaling is reliable only 
when both the known solution and the rescaled solution have sufficiently small ratios 
||. However, in practice, all four parameters $ , $ , m* and m* are varied in the 
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numerical simulations, and the expected scaling will be verified in the following section 
(see Figure 7.4). 





known sol. 


rescaled sol. 


Scales 




fim , \im , 7$o , 7$o 


Background Solutions 


$(t) , $(t) 





Table 7.1: With known solutions to the background equations for $ and $ of the 
form (7.15) at the scales m and M = (|$ | l^ol) 1 ^ 2 . one mav i n f er rescaled solutions 



at the scales /im and jM. The table shows how the solutions and their functional 
dependence change under the rescaling. 



The scaling of the background solutions was simply assumed in Section 3.4 We 
now briefly justify how the background rescales. We use the complex representation 
of the fields and physical time for which the equation for a flat direction is, 



m 



10|$r 



<T = 



(7.15) 



and where the assumed form for A in (7.4) has been substituted. Here the Hubble 
parameter is driven externally by the dominating inflaton so we may take H = which 
is the time evolution for a matter dominated universe. With these two substitutions, 
the rescaling of the solutions is straightforward to show and the results are outlined 



in Table 7.1 This Table may be compared to Table 3.1 constructed earlier, and 
note in particular, the time dependence has rescaled in the same way in both tables. 
What the above rescaling means is the trajectory of the background retains its shape 
under changes in m and M, with the only changes being the overall scale of the 



trajectory and how fast the trajectory is traversed. Also, from Table 3.1 one reads 
the occupation numbers will rescale as follows, 



{ni , n 2 , n 3 } 



rii , n 2 , n 3 



(7.16) 
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and applying this along with the rescaling of the background and the rescaling of the 



frequencies uoi, it is possible to show that r pro d defined in (7.13) will rescale as, 

(7.17) 

which is the main result of this section. This formula implies that r proa - may be 
expressed , 

mm „ / m |$ I \ 
|$ ||$o| \m \%\ J 

where the the dependence on the overall scales of the masses and the VEVs is made 
explicit. The function / may be obtained numerically, and in general it exhibits an 
exponential growth with the number of rotations N as well a strong dependence on its 
arguments |$ Q |/|<& | ar| d m/m. However both the exponential growth and the scaling 
will only be obtained within a band of parameter values for |$ |/|$o|, and outside this 



band, neither will be obtained. Within the resonance band (7.18) can be cast in the 
form 



Vo^T^n 10 CTAf , (7.19) 

where C and a are two time-independent quantities which are functions of the ratios 
m/m and |$ |/|$ |. One may also invert this expression to obtain, 



N decay ~ - log 10 . (7.20) 




And thus, so long as the parameters are within the resonance band, the decay time 
Ndecay is driven mostly by the the production rate a with other effects contributing 
small logarithmic corrections. 



7.3 Numerical Results 



The numerical model was presented in the first subsection along with some useful 



physical measures such as the adiabaticity matrix in (7.12), and the ratio of the 



energy densities r proc i in (7.13). The adiabaticity matrix elements are a quick method 
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for identifying a resonance band, but the ratio of the energy densities will provide a 
more precise method, and the number of rotations Ndecay at which the ratio becomes 
unity will be our primary measure of the decay time. In the previous section it was 
determined that only the ratios m/fh and |$ |/|^o| nave anv physical relevance in our 
numerics, so in the following, the results will be shown in terms of these two ratios. 



We begin by showing the background evolution in Figure |77T] for a sample set of 
initial conditions and parameters. As discussed in the previous section, the background 
trajectories retain their shape under a rescaling of the masses m and the VEVs. Thus, 
the trajectory shapes shown are relevant to any mass scale including the Electroweak 



scale. The dimensionless charge asymmetry (4.40) may also be determined, and 



for instance the flat direction $ shown in Figure |7.1| is computed to have a charge 
asymmetry £b-l ~ 0.26. In the following, the relative size of the trajectories is varied 
simply by adjusting the initial value of the ratio 1^4. 

Now with just information of the above background evolution, one may compute 



the eigenvalues (6.22), Gamma matrix elements (6.22) and also the adiabaticity ma- 



trix elements (7.12). The adiabaticity matrix elements are a good predictor for when 



particle production may occur. In Figure 7.2 the Root-mean-square leading adiabatic- 
ity matrix elements are computed over the course of five rotations of the flat direction 
$ while varying the ratio of the initial VEVs. The RMS values of A+j suggest 

there is indeed a resonance band. One notices that the predicted band is weighted 
more towards larger values of |$ |/|$ | in the figure. This may be explained by the 
fact that the heavier field $ becomes unfrozen at a later value of the scale factor 
than the lighter field $ and this results in an effectively smaller ratio |$|/|$| once 
both fields are oscillating. For example, although the ratio is approximately 15 for the 



trajectory shown in Figure 7.1 the effective ratio when both fields are unfrozen and 
oscillating is clearly closer to unity. 

One may now solve the Heisenberg Equations of motion for the same range of 
parameters to confirm the existence of a resonance band. To obtain some further in- 
tuition for the particle production, the total occupation number n tot = n 1 + n 2 + n 3 is 
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Figure 7.1: The evolution of the flat directions for the choice of masses m* = 
10~ 6 , fh/m = 3.72 and initial phases S = 0.25 , S — 0.156. The initial values 
of the VEVs are chosen such that |f» |/|<I>o| = 15 and |$ | = 10" 2 M p . We show 
the real and imaginary parts of the flat directions for the first three rotations of |<£> | 
(solid red) and |$| (solid green); the right panel shows the same evolution as the left 
panel, but magnified. The trajectory shown for the field $ is chosen to match that of 
Figure 14.31 hence the charge asymmetry is approximately the same e b—l ~ 0.26. 
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Figure 7.2: The root mean square (over physical time) of the four leading adiabatic 
matrix elements A 2 2, A 33 , Ai 3 and A23 computed according to (7.12). The compu- 
tation is performed for the modes with momentum k — 0, for five rotations of the 
flat direction $. While the RMS values are not directly correlated with the growth 
exponent a shown 



in Figure [776 
production may be possible. 



they indicate the parameter space where particle 
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shown in Figure 7.3 compared to the above background evolution. Notice particle pro- 
duction occurs for short periods of time when the fields have comparable magnitude. 



One may see this feature in the expressions for the elements of the T matrix (6.22). 



The appearance of factors such as |$|/ y |$| 2 + |$| 2 and |$|/y |$| 2 + |$| 2 in these 
expressions indicate a suppression when the ratio $/$ is very large or very small, but 
the T matrix, and thus particle production, is unsuppressed for those instants when 
$ ~ $ is obtained. Finally, we should mention that the particle production does 
not always occur in brief bursts. Specifically, when the fields are continually of a 
comparable magnitude, the production can then be a more continuous process. 

Next we check the rescaling discussed in the previous section. The Heisenberg 
equations are solved for different values of the parameters m*,m*,$ , and $ , and the 



occupation numbers are rescaled according to (7.16) in order to compare the spectra 



at a given time. This procedure is illustrated in Figure 7.4 where we show the spectrum 
of the first (heavy) eigenstate for three distinct cases in which the ratio || is varied 
over the range 1CT 7 < jj < 1CT 5 . The three spectra coincide once ri\ and k are 



rescaled according to Table |3.1| Although we do not show it here, we also verified 
that the occupation numbers for the light eigenstates do not change with 

Next the time evolution of r pro d is computed, and this is done for three values 
of the ratio m/fh and a sampling of the values |$o|/|^o| over the resonance bands. 



The scaling relation (7.17) is applied to determine the number of rotations at which 
r P rod = 1 is satisfied with m = m susy and M = 10~ 2 Mp, and these values are 



recorded. The results are shown in Figure 7.5 



The data of Figure [775] exhibits a resonance band for each of the three selected 
mass ratios. Note also that different values of |r were used to obtain this data. Thus 
the figure is a second verification of the scaling relation. Before rescaling there was 
noticeable scatter in neighboring data points, and after rescaling, the data collapsed 



to the smooth curves shown in the figure. The growth exponent defined in (7.19) was 



also computed for the same parameter set and these are shown in figure 7.6 



The data shown in the Figures 7.5 and 7.6 are the main quantitative results of 
this thesis. They explicitly show the decay of our toy model flat direction for realistic 
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Figure 7.3: Upper panel: Number density of the produced quanta for a given 
momentum = 1CT 7 . Lower panel: amplitudes of the two flat directions. We 
note that particle production occurs whenever the two amplitudes are compara- 
ble, as discussed in the main text. The parameters chosen for this evolution are 

m* = 10" 6 , F* = 2, fh/m = 3.72, |$ |/|$o| = 15 
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Figure 7.4: Occupation number for the first (heavy) eigenstate after three rotations 
of the first direction, for a specific set of parameters: |$ | = 10~ 3 M P , l^ol/l^ol = 4, 
rh/m = 3.72, S = 0.25, S — 0.156. The larger figure shows the occupation number 
as a function of the momentum for three specific choices of the flat direction mass m*. 
The insert shows the same result, but plotting rii/m* vs. k*/m*. In terms of these 
variables, the three curves overlap. Since m* oc ?f, this confirms the scaling (7.16) 
for the occupation number, and the fact that the momenta of the quanta produced 
also scale as ^. 
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Figure 7.5: Number of rotations of the first VEV where the production criterion 
fprod = 1 is satisfied and its dependence on the initial ratio of the VEVs |$ |/|$o|. 
for three different mass ratios. The numerical analysis was made for the first twenty 
rotations of the first flat direction and the results were rescaled to m/e = lO^GeV 

and ^/|$o||$o| = 10- 2 M p . 
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Figure 7.6: The dependence of the growth exponent a defined in (7.19) on the initial 
ratio of the VEVs |$o|/|^o|. f° r three different mass ratios. 
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values of the parameter space. The degree to which these results may carry over to 
MSSM flat directions is discussed in the Conclusions which follow. 
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Chapter 8 



Conclusions 



Inflation provides a strong motivation for the study of the dynamics of coherent fields 
in the early universe. The inflaton itself is a coherent field, and as we have shown, 
other fields such as supersymmetric flat directions may be made coherent through the 
process of inflation. The motivation for this thesis and the attempt was to determine 
the evolution and decay of flat direction VEVs in the early universe. As discussed, the 
VEVs will break at least some of the standard model gauge symmetries and establish 
an effective vacuum into which the inflaton would decay and in which the inflaton 
decay products would thermalize. However, evidence was presented that flat direction 
fields may decay quickly and nonperturbatively, and thus restore the vacuum to a gauge 
symmetric state perhaps even before the inflaton has decayed. The final thermalization 
of the inflaton decay products was not studied, but once the gauge symmetric vacuum 
is restored, existing techniques such as those of [85] would be applicable. 

Our study began with a very brief introduction to supersymmetric theories and 
the scalar potential of these theories. It was shown that the scalar potential - which 
is a functional of all the scalar fields of ones model - can possess special nonzero 
combinations of fields for which the potential is exactly zero. These flat directions 
are the light scalar fields described above. We thus studied the growth of the fields 
during inflation, and established reasonable upper bounds on their VEVs as shown in 



Figure 4.1 After the conclusion of inflation while the universe is behaving as a matter 
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dominated universe, a special time was noted when the Hubble parameter became 
comparable to the mass of the light scalar. At this time, the complex scalar would 
begin to execute well defined oscillation. The motion is generally an elliptic-spiral 



like trajectory in the complex plane as shown in Figure |4~3] or |7.1| with the frequency 
of rotation being simply the mass of the field. The orientation of rotation and the 
ellipticity of the trajectory will additionally quantify a conserved charge for the field 
which in the case of the MSSM flat directions is typically the difference of baryon 
number and lepton number or B — L. During this time both the inflaton and the flat 
directions both oscillate and reduce in amplitude as matter constituents and as shown 



in Figure |4T2] Both will eventually decay however, and the order of their decays is 
somewhat important. If the inflaton decays first, and into a relativistic species, then 
the flat direction may overtake the inflaton's decay products later. However if there 
are multiple flat directions which have obtained obtained large expectation values, 
then the flat direction VEVs may decay first through a rapid nonperturbative process. 

An example of the nonperturbative decay was shown here. It is shown in a toy 
U (1) gauge model that the decay can occur in approximately ten or twenty oscillations 
of the flat direction fields for a range of parameter choices. This is in contrast to the 
perturbative estimates which would put the flat direction decay to be long after the 



decay of the inflaton. The Figures [775] and |7T6] are thus the main quantitative results 
of this thesis. 

In the following we examine some notable aspects of our models. To begin, we 
note our study [27J was the first concrete calculation of a nonperturbative decay in a 
gauge model. One consequence of the gauge symmetry on our results has been the 
gauge current constraint. For instance, in the model of a single flat direction, the 



gauge current constraint (6.7) essentially would not admit a time dependent solution 
except in the circumstance when the masses of the two fields involved were degenerate. 
We were thus led to select this case in order to consistently apply our formalism. The 
general case appears to only admit a space-time dependent solution, and to pursue this 
model apparently would require lattice simulation. Similarly, in the model of multiple 



flat directions of Section 6.3 which possesses three flat directions, the gauge current 
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constraint forced one of three flat directions from obtaining a VEV in the solution. A 
simple lesson from our study is that the current conservation constraint can often have 
a non-trivial effect on the dynamics. We note that similar conclusions were drawn in 
the work of [87J. 

In the above, the conservation of current was considered on just the background. It 
should be possible to consider gauge current conservation at the level of the perturba- 
tions as well. Specifically, the conservation should appear in the Heisenberg Equations 
of motion. These equations do not conserve current on their own however, so it must 
be that the driving terms which are derived from the U(l) symmetric Lagrangian will 
contain this information. Verifying this conservation on the perturbations is another 
direction which may be pursued. 

Having studied the U(l) models, a logical next step would be to pursue the same 
analysis on an MSSM flat direction. As an example consider again the work of [S7] 
in which the background evolution of the H U L multi-flat direction was solved. These 
authors obtained only time-dependent solutions to this model. However, using the 
methods presented in this thesis, the spatial stability may in principle be studied^] In 
particular, for the H U L directions or some other choice of MSSM directions, obtaining 
the driving terms to the Heisenberg equations would be more involved due to the added 
gauge symmetry, but the series expansion in the hierarchy jj for example would greatly 
facilitate these calculations. In general, the techniques learned by studying the U(l) 
models here can make these more complicated problems feasible. 

Putting aside the above modelling considerations, we consider again the connec- 
tion between flat direction decay and the equipartition principle which was made in 



Section |3.4| and in the Results. The most significant consequence of equipartition is 
the production of heavy quanta. Our analysis and numerics both revealed that heavy 
quanta may be produced in the flat direction decay. The production of heavy quanta 
is not uncommon in preheating scenarios as was discussed for a toy model of inflaton 



1 The classical development of spatial instability in an Affleck-Dine condensate has been studied 
for example in [88J where the connection is made to non-topological solitons in gauge theories known 
as Q-balls [89J (see also |74|). Due to limitations of time, the relation of this thesis work to the work 
of [88] has not been explored. 
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preheating in Section [5] and has also been shown in more detailed studies (for example 
[861 l6l|). One might speculate that equipartition is a property that all preheating 
models which involve heavy states might share in common. 

However, it is unsettling that states with nearly Planck mass could be produced 
in such large quantities. One must look to the expected decay time of the heavy 
quanta to resolve the difficulty. A preliminary estimate of the decay time was made at 



the conclusion of Section |5.7| and the estimate indicates a decay time which may be 
comparable to the mass of the heavy quanta. In order to treat these states consistently 
as particles, we had to assume the decay time was somewhat longer than one period 
of oscillation of the heavy state. This time would still be very short compared to the 
period of a flat direction state, and in particular the ratio of the two periods is simply 
This suggests the heavy state occupation numbers would be depleted quickly 
on the time scale of our simulations. Equivalently, this means that the inclusion of 
interaction terms would be significant for the heavy states. Further study in this area 
is warranted. However, we note that the exponential production on our model did 
not greatly rely on the presence of large occupation numbers for the heavy state. 
This was checked by artificially setting to zero the mixing between the states via 
I = J = so that the only source of particle production could be from the evolving 
eigenfrequencies and not the evolving eigenvectors. After this change, an exponential 
particle production was still obtained for the light sector of the model. 

The above situation in which the final states are short-lived is an interesting one 
and to our knowledge has not been addressed in the literature. However, the effects 
of interactions and the back-reaction of produced particles on the condensate have 
been studied for example in [S5], so analysis of a similar nature may be possible. 
Alternatively, the issue may be resolved through classical lattice simulation. In gen- 
eral, classical lattice simulations may be useful to compare against the generic results 
obtained here. For instance, classical stability analysis such as the determination of 
lyapunov exponents may be applied in combination with a lattice simulation (5TJ. In 
principle, methods such as these may provide further insight. 
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Finally, we present some discussion of observables. Our model suffers in this re- 
spect in the same way as other models of preheating. The difficulty lies in the process 
of thermalization itself which erases most all information about what happened during 
the reheating process aside from gross thermodynamic quantities such as the baryon 
asymmetry and the reheating temperature. However other observables may survive 
the reheating process. These include such effects as modulation of the cosmological 
perturbations, gravitational waves, the creation of dark matter, primordial black holes 
and topological defects [63J. The most likely place one may look for a signature of the 
existence of a flat direction VEV may be in the statistics of the cosmological perturba- 
tions and in particular the imprint left on the Cosmic Microwave Background. In some 
specific models known as D-term inflation in which the inflaton is supersymmetric (see 
|7J for an introduction), it has been noted the flat direction VEVs may contribute to 
the cosmological perturbations [921122]. For instance, bounds may be placed on the 
magnitude of the flat direction VEV (based on the COBE normalization) which are in 
the vicinity of 10~ 3 Mp to 10~ 2 Mp and which correspond to a lowest non-vanishing 



terms in the superpotential n = 4 and n = 6 [92J. For comparison, see Figure 4.1 



More recently, it has been suggested that the effect of parametric resonance of a 
field coupled to the flat direction may generically amplify the perturbations of this 
field so they may be observable in the CMB as curvature perturbations [93J. Again, 
in the absence of a direct signature, one may constrain proposed models based on 
non-observation of such a signal [53]. We note this last work appears to have been 
somewhat motivated by our own work on which this thesis is based [27J. 

This completes the main body of this thesis on the dynamics of supersymmetric 
flat directions. While much has been determined in this study, we have emphasized in 
the above that there remains more work that can be done on this interesting subject. 
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Appendix A 



Notation 

The notation used in the thesis is here summarized for the convenience of both spe- 
cialist and non-specialist readers. 

The notation for fields is as follows: Fundamental complex scalar fields (specifically 
the fields of which the flat directions are composed) are denoted by 0j = /je * with i 
running over the whole set of fields. However, the inflaton field is denoted by \ (the 
author apologizes as this is an unusual convention). The U(l) vector field is denoted 
by B^. The SU(2) and SU(3) vector fields are A* where a is a group index, and it 
should be clear from the context which symmetry is being considered. Fermions are 
denoted by ip. 

The fundamental fields are typically decomposed into a classical background plus 
quantized fluctuations away from this background. The classical backgrounds are 
written in upper case and the perturbations in lower case. In the case of the vector, 
the classical background is written in calligraphic letters, for instance — + 
where are the quantized fluctuations. Note the background value of the vector £> M 
is always zero in the context of this thesis (See Section [6]). Finally, the notation for 
derivatives is, = (d^xY = {x,nY 

When discussing the classical background of the scalar fields $j, the terms Vac- 
uum Expectation Value (VEV) or scalar condensate are sometimes used as these are 
equivalent descriptions. When the term VEV is used, it is implied that it is the VEV 
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determined by a local observer in some causally connected patch of the universe at 
some time, (0i) pa tch = 3>i, and this is not to be confused with the VEV determined 
over "all" space which is zero (0j) a ii space = 0. We sometimes use the term "Flat 
Direction" when referring to the VEV, although strictly speaking the flat direction is a 
property of the scalar potential while the VEV is a realization of the field configuration 
on the potential (see Section [2]). Effective vacuum and Bose Einstein condensate are 
also used to describe the classical background. 

In describing the "decay" of the flat direction condensate, we sometimes refer 
to this as "decoherence" or "the development of spatial instability" as these are all 
equivalent statements for the decay in the context of our models. 

The notation for energy scales is as follows: The mass of the scalar fields which 
comprise the flat directions will be m or fh and these will always be of order m susy ~ 
TeV = 1000 GeV which is the assumed scale of the supersymmetric particle spectrum. 
The TeV energy scale (approximately 1000 times the mass of the proton) is heavy by 
the standards of nuclear and particle physics, but it is relatively light compared to other 
scales in Cosmology such as the Grand Unified Scale, M GU t ~ 10 16 GeV or the Planck 
scale, M P ~ 10 19 GeV. Note the Planck mass is defined by Gm = -j^r- This is not to 
be confused with the reduced Planck mass Mp = M P j\fM ps 2.4 x 10 18 GeV which 
is not used in the thesis. An arbitrary heavy mass scale will be written as M, and in 
the results section M is defined specifically as the scale of the VEVs |$j| ~ 10 _2 Mp. 
The ratio of the TeV scale states to one of these heavy scales is then generically 
written ||. A quantity is unsuppressed when it is zeroth order in the expansion in ||. 

The gauge couplings are denoted by "g" in Section [2] where supersymmetry is 
discussed, but in our model and results sections [6] and [7J we instead use e. This is not 
the electron's charge as it would be in Quantum Electrodynamics, but the notation is 
convenient. 

The terms parametric resonance, non-adiabatic evolution and nonperturbative dy- 
namics are often used interchangeably. 

Finally, the general metric tensor is g^ u and the Minkowski metric tensor is rj^. 
The cosmological scale factor is R. The conformal time is 77 with time derivatives 
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written as a prime (such as R'), and physical time is t with time derivatives written 
with a dot (such as R). 
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Appendix B 

Details Regarding U(l) Gauge 
Symmetric Actions 



The action for a U(l) gauge field, and n complex scalars charged under this gauge 



group is presented in the main text (6.1). The procedure for gauge fixing such actions 
to the unitary gauge was also presented there. Some details which were not explained 
in Section [6] are explained here. In particular, it is shown here (i) how to incorporate 
the expansion of the universe, (ii) how to decompose the vector field Lagrangian and 
eliminate the nondynamical time-like component B and (iii) how to diagonalize the 
Lagrangians for the systems of two fields and four fields in the expansion parameter ||. 
Incorporating the Expansion of the Universe: With a change of spatial co- 



ordinates and a rescaling of fields one may put the Lagrangian (6.1) in a form which 
is more easily quantized. The background metric is taken to be the Friedmann- 
Robertson-Walker metric with the spatial curvature constant set to zero. The line 
element for this metric is ds 2 = dt 2 — R(t) 2 dx 2 , where R(t) is the scale factor. It 
is convenient to use the conformal time coordinate i] which is related to the physical 
time by Rdr/ = dt. The line element in these coordinates is, 

g^dx^dx" = R{t) 2 r]^dx^dx" = R(tf(dt 2 - dx 2 ) 



where r/^ is the Minkowski metric. After choosing the conformal time, the action (6.1 ) 
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is written in a way that the raising and lowering of the space-time indeces is simply 
done with the minkowski metric as follows, 

S = J d 4 xR 4 {- J^/^F^v + ^(D^diD^y - Vfa, ..4r 

Additionally, if we redefine the scalar fields <\> = \, the action rewrites, 



s = Jd A x + p^PW + - R A v ( 



R 1 '" R 



where primes above denote where the scale factor now appears through the effective 
mass term and also may appear in the potential, though the potentials we consider 
such as that from the D-term are quartic in the scalar fields, and so i? 4 4 = y? 4 and the 



scale factor drops from these terms as well. To summarize, for the Lagrangians (6.1), 
the prescription for incorporating the scale factor is to replace the mass terms as 
follows 



and when later computing physical quantities to convert back to comoving coordinates 
via Rdr] = dt and to physical fields via 4>i = %■ 

Decomposing the Vector Field and Eliminating the Nondynamical Time- 



like Component: We now decompose the vector field in our U(l) actions (6.1). 
We note that is coupled to the scalars through the current as well as through 
the effective mass term A4 2 both defined in (6.3). The part of the lagrangian which 
contains the vector is thus, 

where the field strength tensor is F^ v = d^B u — d u B^. We split B^ into it's space 
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and time components B and B { as follows 

^massless ^F^F^" 

— 2 ~ B o,i f - | (Bj,i -Bi,j ) 2 

= Vif' + ±B°,i B° n +B*'B° n ~B*,i B\, +±B* , t B\ 3 

= X -B' -B' + 1 V5 • V5 + 5' .V5 -^|V xB\ 2 

where the identity e^e^ = (5ji5 km — 5j m 5 k i) has been applied. The Lagrangian for 
the massive vector field is then 

\B'\ 2 + (Vfiol 2 + 2B' ■ VB - |V x Bf 
+ M 2 [B 2 - \B\ 2 ) - 2J B + 2 J • B 



(B.l) 



We split the vector into a transverse and longitudinal part B = B T + VL where 

— * 

V • B T — 0. The transverse part of the vector decouples from the Lagrangian since 
it transforms under the Lorentz group as a spin-1 field, while the scalar longitudinal 
mode of the vector as well as the fundamental scalar fields transform as Lorentz scalars. 
Hence the transverse mode will decouple at quadratic order in the perturbations. We 
thus only consider the longitudinal mode V L as well as B which is non dynamical and 
will be replaced by substitution of its equation of motion. The equations of motion 
for B^ are 

F^+M 2 B„ = j„ 
UB U -{B^) )V +M 2 B U = j, 

and from these, the equations of motion for B and Bi are, 

- V 2 B -d t (V -B)+M 2 B = J (B.2) 
-UB - V(d B - V • B) - M 2 B = J (B.3) 
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We solve for B by inverting the equation, 



'-V 2 + M 2 )B 



Jo + (V 2 L') 



The part of the Lagrangian (B.l) corresponding to the longitudinal mode is then 
determined to be, 



- {|VJ'| 2 - B [J + V 2 L'] - M 2 \VL\ 2 + 2J ■ Vj} 



(B.4) 



The remaining aspects of the calculation will be done in momentum space for which 
the expression for B is, 



J - k 2 L' 
k 2 + M 2 



Bo 



(B.5) 



where k = \k\, The longitudinal part of the Lagrangian (B.4) in momentum space is, 

1 U k 2 M 2 \ , , k 2 (J*L' + J L*') J* Jo 

±j 1j ~\~ ~ 



'o J o 



k 2 + M 2 k 2 + M 2 

- M 2 k 2 L*L + ( J*L + JL*) ■ k 



k 2 + M< 



The longitudinal mode can then be made canonical with the field redefinition, 

Vk 2 + M 2 
kM 



(B.6) 



(B.7) 



Notice in the simple case where the current is exactly zero, the Lagrangian for the 
longitudinal mode of the vector after applying the field redefinition is, 

1 



VU - (AT + k 2 )L*L 



(B.8) 



which is the expected result for a freely propagating longitudinal mode of the massive 



vector field. The above normalization (B.7) will be performed on the specific examples 



which follow in the case of nonzero current Jo. The longitudinal mode should still 
obtain the above dispersion relation tu 2 = k 2 + M 2 when the fields are evolving 
adiabatically. 

To this point the Lagrangian has been written in the fundamental field descrip- 
tion. However for the purpose of computing observables, the Lagrangians are always 
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expressed in background fields plus perturbations. The models we are interested in 
have backgrounds which satisfy (B^) = and also ( J ) = 0, so we must only consider 
the linear perturbations in these quantities, 6 M and in order to obtain a quadratic 



Lagrangian from the above result (B.6). Additionally, we will employ the expansion 
parameter || which was also applied in Section 3.4 Our objective is to obtain the 



quantities and coi to leading order in this expansion. 

Diagonalization of the Two field Model: The Lagrangian for the perturbations 



is (6.12) from which we read the current density and effective mass of the vector to 
be, 



Jo 



\eFT!8 



j = 



M 



\eF 



The Lagrangian (6.12) after the elimination of b is, 



I [(<V) 2 + (<V) 2 + (<V) 2 ] + (iV - fa) 



1 

2 

1 (\ 



1, 



_ e 2 F 2 m 2 3 E/2 ^ 2 

2 V 4 + 4 



(B.9) 



where the vector field is not shown since it has been decoupled. The remaining coupled 
system involves / and a and these fields are obtained to be degenerate in mass m. 
They are simply the flat direction perturbations. 

Diagonalization of the Four Field Model: In the four field model, the current 



density of (6.14) and the effective vector mass are obtained in terms of the perturba- 



tions (6.19) to be 



Jo = e 



(F 2 - G 2 )6' + 2FZ'5 + 2Gt'~8 
J = e {F 2 - G 2 )V6 
M = eVF 2 + G 2 
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Substituting these into (B.6), one notices there will be coupling between L and 9 in 
their kinetic terms at leading order. Our procedure will be to first find the transforma- 
tion which makes these kinetic terms canonical (temporarily discarding mass terms and 
mixed terms). We then apply this transformation to obtain the complete Lagrangian 
in the form Q3.28 ). The kinetic terms are, 

/ bnl> 2 

k 



VU + 



L*'V + 9*'B' 



Vk 2 + M 2 \M 
k 



Jo f, + 3o_ L *, , _ 



M 
M 



M 

F 2 -G 2 



l + (£) 2 



M 



k \2 



+ (F + G 



'H + e'L* 



AF 2 G 2 + (F 2 + G 2 ) 2 (A) 



where j has been substituted and the kinetic term for 6 has been normalized. This 
last expression is compactly written in matrix notation in terms of the fields L and 9 
as follows, 



1 



(L*,9*) 



n L 
n e 



1 Q 
Q 1 



n L 
n e 




where 



Q 



G 2 



M 



k \2 



AF 2 G 2 + (F 2 + G 2 ) 2 (A) 

4F 2 G 2 + {F 2 + G 2 ) 2 (^) 2 
\ (F 2 + G 2 )(l + (±) 2 ) 

kM 
Vk 2 + M 2 

One may apply the following additional rotations and rescalings 

L, 



VT+Q 

o V^Q 



Rq 


n L 









(') 










(B.10) 
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where R is a time independent rotation and Ri is a small time dependent rotation, 



Rt) 



Ri 



q 



1 










I) 


" i 


k 2 q 2 
2 


kg 




kq 


1 - 




-F 2 + 


G 2 



2 



AeFGVF 2 + G 2 

and these rotations have been chosen specifically to decouple the longitudinal mode 



of the vector from the other scalar degrees of freedom. The expression (B.10) may 
be inverted and expanded finally to obtain, 



9 



2FG 

k 



G 2 



2eFGVF 2 + G 2 



e c + 



k ( m 2 \ 

2e 2 (F 2 + G 2 ) c+ yM 2 ) 



(B.H) 
(B.12) 



This final result is substituted into the Lagrangian (B.6). The remaining mixed scalar 
degrees of freedom are arranged in a vector X = (8, S,9 C ), and the Lagrangian for 



these has the desired form (3.27). Specifically, the K matrix, the matrix (Q 2 + K T K), 



the eigenvalues uji and the matrix Y are all shown in the main text (6.20 6.22) 
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Appendix C 



Calculating the Variance of a 
Scalar Test Field During Inflation 



We determine the variance of test scalar field during an asymptotically de Sitter 
expansion. The sampling is done inside a patch of space of width L. We are not 
interested in wavelengths smaller than L, so we smear out these fast fluctuations by 
use of a window function Wl(x), 

1 



L (t)= / d 6 xW L {x)4>{x,t) 



W L (x) 



-e 2z? 



(C.l) 



(2tt) 3 / 2 L 3 ' 

in which Wl(x) is any reasonable normalized distribution of order 1 inside the patch, 
and dropping to zero quickly outside of the patch. In the above, it is taken to be a 
normalized Gaussian. The quantity (fi L is then the quantity which we must determine 
the variance of. The calculation is performed using conformal time dt = Rdr], and a 
redefined scalar field <\> = -|. The fields are decomposed as in (4.18), and the variance 
is thus obtained to be, 



<0| 

1 

2 
1 



d 3 x W L (x)if 
d 3 k 



(2tt) 3 
L ~ l d 3 k 
(2tt) : 



10} 

\h k \ 2 \w(kL)\ 2 



■N 2 



(C.2) 
(C.3) 
(C.4) 
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where w(kL) is the Fourier Transform of the window function, and we've used the 
fact that w(kL) is order 1 for \k\ < L v and drops to zero quickly for \k\ > Lr 1 . 
The effect of the window function is simply to impose a momentum cutoff. In the 



de Sitter metric, the equations for the modes is listed in the main text (4.23) which 
again is 



hl + 



k 2 + 



m 



- 2 



rf 



h k = 



(C.5) 



and may be put in the form of Bessel's equation with the definitions x = —kr], and 



x: 



dx 2 

the solutions are, f(x) 
form, 



dx 



9 
4 



, =~-:-— (C.6) 
i z 

BhYJx), and so the mode functions take the 



h k {rj) = yj-kr) (A k J n (-kr)) + B k Y n (-kr))) 



n 



m 
IP 



(C.7) 



The constants A k and B k must now be determined from the initial conditions for the 
modes which are fixed by the quantum fluctuations. Technically this means we must 



canonically normalize the modes h k at early times (or large k). The equation (C.5) 



for the modes at sub-horizon scales, asympotes to h'l + k 2 h k = with solutions e ±tkv . 
Referring to (4.18), we want h k (r]) = ~^ e ~ lkr ' at eai "ly times for our fields to be 



canonically normalized. Our solutions (C.7) are then normalized as follows 



h k {fj) 



-TITX :( nir , ir 

—e l 2 +4 



4 



n (-kr)) + iY n (-krj)) 



(C. 



where one may apply the asymptotic form of the Bessell functions which is, 

(J n (-kri) + iY n (-kri)) r- 



7i(—kr]) 



n,7T 7T 

2 4, 



for large (-kr]) (C.9) 



from which it is quick to verify that the solutions (C.8) have the correct asymptotic 
form far within the horizon. We consider these solutions in the limit of small mass 



Hi 



< 1 where n 2 > 0, and large mass > 1 where n 2 < 0, and we note the 
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Bessell solutions will still hold in the latter case when n is pure imaginary. We are 
interested in calculating the variance, so we substitute in our general solution, and we 
approximate our solutions in the long wavelength limit as our integration will only be 
over scales larger than the horizon. The steps of the calculation are shown first for 
the case n 2 > 0, 

" L_1 d 3 k 



(<Pl(v) )™ 2 >o 



(2ny 
dk 



\h k \ 



k 



(2tt) 2 '~ 4 

phys dkpfoyg kphys 



J n {-kr,)+iY n {-kr,)\' 



J l) 



v phys 

~h7 



+ iY n 



phys 



'Hi 



(2tt) 2 4ifj 

{HiL p hys)~ l dx x 2 

\J n (x) + iY n (x)\ 2 



Hr 



(2vr) 2 4 



(HlLphys) 1 ^ 



X 



(2vr) 2 4 



2 n r(n) 



7T 



2 (2n-5) r ( n )2X r(H lLphys ) * ^ 



7T y 



-X 



(3-2n) 



for x < 1 



where we have put the expression in terms of the physical momenta, k = k phys R, 
and original fields = </?/ ' R, then put the integration over the dimensionless variable 



x 



^phys 



We may similarly calculate the case n 2 < 0, 



_d?k_> h ,2 

2./n (27T)3 |ftfc| 



2i 



2 <0 



a (27T) 

J phy S dk nh „ x k. 



phys ^phys -\n\ir 



J, 



i\n\ 



v phys 

~b7 



i\n\ 



(2tt) 2 AHj 

1 riHiLphys)- 1 j 2 



J phys 

w 

2 



1 



16|n|7r 3 



(2tt) 2 4 

(H I L phys )- 1 ^ x 



x for n 
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Appendix D 

Inflaton Oscillations After 
Inflation 



Here, the approximate evolution of the inflaton field at the conclusion of inflation is 
determined in the case of a quadratic inflaton potential. Note, the result has been 



used in (4.33) of the main text. Also note that the approximate evolution of the 



flat direction fields (4.36) used in the main text may be determined through a similar 
method as that shown below for the inflaton (see also |20J). 

To begin, the Friedmann equation and the equation of motion for the dominating 
isotropic and homogeneous inflaton field are, 

H 2 =(l) 2 = S -f(lx 2 + V( X )) (D.l) 
dv 

dx 



X + 3H X + f - = where V = \m 2 x 2 (D.2) 



At late times, after inflation as occurred, this field is oscillating around the minimum 
of its potential (V — 0, % — 0). The motion is determined by the above equations, 
and one assumes a solution of damped oscillations, 

X (t) = X(t) sin[mr] X(r) = ^ + + ... (D.3) 

where we've defined the time parameter r = (t — t ) with t being the onset of the 
inflaton oscillations. Specifically, x( T = 0) = and this constraint requires X(r) 
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must not contain even powered terms. We may solve for the series Ai,A 3 ,... by an 
iterative procedure. Substitute the trial solution into (D.l), and set At = A. Also set 
smmt = s and cosmt = c to simplify the notation, 

-Ar 2 + 0(r 4 ) 
2Ar 3 + o(r 5 ) 

Xs + mXc 

Ast~~ 2 + mcAr 1 (D.4) 
Xs + 2mXc - m 2 Xs 

-2mcAr 2 - m 2 Asr 1 (D.5) 



X 
X 

X 

X 



where only the lowest order terms have been retained. Then, from (D.l) 

AttG 



H 2 



3 

AnG 
AtiG 



(Xs + mXcf + m 2 X 2 s 2 



m 2 X 2 (s 2 + c A ) + 2mXXsc 



m 2 X 2 



AnG 2 A 2 

m —7T 

3 t 2 



2sc 



X 
mX 



1 - 



2sc 
mt 



Now substitute for x< X ar, d H into (D.2) and keep track of the order t n . 

X + 3Hx + m 2 x 








-2mct" 2 - m 2 Asr l + 0(t" 3 )) 



+ 3 



4 



nGm 2 



1/2 







Ar 1 + (hi 

+ m 2 s{Ar 1 + 0{r 3 )) 

(-m 2 As + rr^As^- 1 + I -2mcA + 3mcA 2 



+ mcAt~ 






l/2\ 


^nGm 2 






r 



(D.6) 



the first term vanishes, indicating that the guess for the solution was correct 
condition for the second term to vanish fixes the constant A, 

A = A 1 ' 



\J"iTxGm 2 
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and the 

(D.7) 



Considering higher order terms would yield the coefficients A 3 ,A 5 ,..., which will not 



be done. Substituting (D.7) into (D.6) yields the Hubble parameter to lowest order, 

/:/ ^ - (D.8) 



2 
3t 



which is the Hubble parameter of a perfect fluid of nonrelativistic and pressureless 
matter. 
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Appendix E 

Thermodynamics of a Relativistic 
Plasma 

Here we present some expressions used in the main text for a universe dominated by 
an isotropic and homogenous relativistic fluid in thermal equilibrium. First, the total 
energy density of the fluid is related to the temperature through the Stefan- Boltzman 
Law, 

7T 2 7 
P ~ Prad = 7^g*Tr ad , 9*(T) = 9b + -g f 

where g* is the effective number of degrees of freedom in the fluid and g b and gj 
are the numbers of bosonic and fermionic degrees of freedom respectively. Using the 
Stefan-Boltzman law in combination with the Friedmann equation, one determines 
the temperature in terms of the Hubble parameter, 
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The entropy density for this relativistic fluid may be determined to be, 

The entropy density may then be written in terms of the energy density using the 
Stefan Boltzman Law, 

27T 2 {30p rad 

45 V 9*^ 

Finally, the number density of a component relativistic species X is, 

nx = ^9xT 3 C(3) ~ 1.2 

where £ is the Riemann zeta function. 
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Appendix F 



Determination of Various Scale 
Factors and Events 



We determine some expressions used in the main text which involve the scale factor at 
different events during the evolution of the inflaton and flat direction. The reader is 



referred to Figure 4.2 as well as to the following table which describes the succession 
of events with R increasing as one moves down the table from the top. 



Event 


R 


Eq. of state 


inflaton begins oscillations 


Rq 








matter dom 


flat direction begins to oscillate 


Rm 




inflaton decays 


Rd X 








radiation dom 


flat direction overtakes inflaton 


Req 








matter dom 



Other scale factors not shown in the table are R^, the scale factor at the instant of 
flat direction decay and Rrh, the scale factor at the instant of thermalization. Only 
the ratios of scale factors are physically meaningful and to determine any ratio, we 
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can apply the following relations, 

R oc t 2 / 3 — > t oc — > i? oc H~ 2 / 3 during matter domination 
R oc t 1 ! 2 — > £ oc — > i? oc H~ x l 2 during radiation domination 

and we will also need to know the energy densities of the flat direction to the inflaton at 
some reference times which we will take to be the onset of their respective oscillations. 

x -m\Ml ' r *-m*& ' m -\ Px ) m r x \Mp) 
where the last expression is the ratio of the energy densities at the onset of flat 
direction oscillations^ For example, using the above tools, the energy density of 
the inflaton at the instant of inflaton decay is obtained, 

Pr x ~ r x {T x M P f 

Similarly, the scale factor when the flat direction overtakes the inflaton R eq is obtained, 

P* , 
— ~ 1 

Px 

/ \ 3 
_m / Rm 

P* Ue, 



,111 I Rrn\ ( R dX 



Px V Rdx ) V ^-eq 



Req ■'- 



R d X r >-- 

which is valid for r m < 1. The condition for which the flat direction dominates is 

R 



(F.l) 



■-eg 



Rd<S> 
Req R d\ 
Rd X R d<l> 

r 



< 



< 



r 



< r 2 m (F.2) 



1 For the inflaton, the energy density at the end of slow-roll inflation is p x w |m^ 2 where 
X = thus r x — g^. For the flat direction the energy density is simply m 2 |$ TO | 2 where \<5> m \ 

depends on the evolution of this field during inflation, thus r$ = 1. In order to trace the dependence 
on these quantities in the resulting expressions, the substitution are not made here. 
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After substitution for r m , this expression defines a boundary in the plane ($ ,r^) 
separating regions in which the flat direction dominates to regions where the inflaton 



dominates which is shown in Figure 5.4 



Next the reheat temperature is determined assuming the thermalization rate is con- 



trolled by 2 — > 3 interactions (with no assumptions on T x and ). From Section 5.3 
the approximate rate is, 

2 



2^2 
2 



r 2 ^ 3 ~«(^) r 

r 2 ^ 2 ~ n x a , a = 
and putting these expressions together, 

m X\ Pdx ( a V f R 



T2^3 ~ a 



~ a 



TrhJ m x \m x J \R dx 
m x V r x( T x M p) 2 ( a V ( R 



T RH J m x \m x J \R dXi 

Now considering the case with no flat directions, and assuming T therm > T x so 
reheating happens instantly at Rd x we obtain the following consistency requirement 

on a, 



« 



m x 



2 r x {T x M P f 



0.1 JT X M P m 



x 



T2-^3 


> 




\ 2 






ft \ 


> 








m x) 




ft 3 


> 





1 m x 
100r x M P 



(F.3) 



Next considering the case involving flat directions and supressed reaction rates, we 
simply replace C^~j ~~ > (§) i n tne expression for a above. We consider three situa- 
tions, Ttherm ~ {r x , T 2 -^3, T$} which can each be realized by tuning $ , ot and/or 
Considering the case where reheating is controlled by T x , one obtains the consistency 
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condition, 



Pd x f a \ 




r2^3 

2 



R 



Pd X 



rn 



x $ 



2 ( Rm\ 

\Rdx) 



r x F 2 x M 2 P a 2 (m 



*8 



100r x a 



m 

Wr 



Mr 



§0 

Mf 



(F.4) 

tp / \lVlp J \ i x - 

Next considering the case where reheating is controlled by T 2 ^3 which happens when 
r$ < < T x . One determines T RH self-consistently using the expresion for r 2 _»3, 

X V PdX ( a Y ( rrh\ 2 



T 2 ^3 ~ at 

T 2 ^3 ~ a 

1 2^3 ~ OL 

1 2 ^ 3 ~ a 



RH 



RH 



2 r x r x M P 



( n \ 



T, 



RH 



rn 



x 



10 7 



Trh 

~Mp~ 

Trh 
~Mp~ 



Trh 
m 



m x 

2 r 2 
; x 1 x 

m x 

2 r 2 
'x L X 



a 



CT)2 / Rr.. , 



V R d x J 



m 



x 



~ r x cr 



M f 
M f 

2 



Rr, 

m 



7? 



r 



(r x a 3 l(T 7 ) 



7 ,i/9 / m x 



m 



M 



p 



r 



X 



Mr 



1/9 



$ 



2/9 



2^3 
5/2 



5/2 



2/9 



Mr 



5/18 



(F.5) 



This result is also continuous and consistent with the previous result. In the final case, 
reheating is controlled by r$, which occurrs when T2-+3 < 

Now specifing to the case when the flat direction decays perturbatively, we deter- 
mine the expressions for T$ first assuming the decay happens during matter domination 
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before inflaton decay. Note that V p ^ rt indicates the instant of decay of the whole flat 
direction condensate which is the single particle decay rate at a specific value of the 
scale factor defined as R d $, 

Tl ert a 2 m 2 

f-NJ 

m $ 2 

m ~ $ 2 \R m J 
Tl ert \ 3 a 2 m 2 



m J §1 

rpert \ / \ 2/3 

%, \ I am \ 



m J V $o / 

Next the case when the flat direction decays during radiation domination and for which 
the above calculation gets modified on the third step which now becomes 

l # ami K d <s> \ I K dx 



m $g \Rd x J \Rm 

-r^pert 2 2 / T-i \ 3/2 / 

1 $ am I 1 v \ / m 



$ 2 ^pperty ^ 
^eri 5/2 ft 2 m 2 / m y/2 

From this last result, we can determine the condition on $ f° r flat direction domina- 
tion, 

R 



Tl ert m 



< 1 



$ " L 2 

— < r , 

m r x 



\ 4/5 / \ 1/5 / \ 2 / ^ 



177 



where we have substituted for T x . We may determine the flat direction decay rate in 
the case the flat direction dominates by, 
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^X 
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Note this result is also obtained by substituting (F.7) into (F.6). Also in the above 



we have computed H eq from (F.l) as follows 
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Appendix G 

Bogolyubov Transformations for 
Multiple Scalar Fields 



The material of this appendix is supplementary to that of Section |3.1| Some calcula- 
tions which were not performed in this section are performed here, and the derivation 



of the equations (3.21) is here made from an alternative perspective. As a reminder, 



the results of Section 3.1 and this appendix were first obtained in [41J, and the notation 
and conventions of this paper are used in the following. 

The starting point is a set of creation and annihilation operators a it a\ which we 
wish to transform to a different set b it b\ through a Bogolyubov transformation, 



OLijbj + j; ; b] 



V 3 



a 



(G.l) 
(G.2) 



where and are complex matrices. We require canonical commutation relations 
on both sets of operators, and this gives constraints on the and /3y coefficients. 
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First calculate the commutator \a,i,a 



ti 



[a h a]} 

[ocikbk + (3* k b{, a* m b ] m + (3 jm b m ] 
a ik a* jm [b k , &jj + f3* ik l3 jm [b{, b m ] + a ik (3 jm [b k , b m ] + !3* k a* m [b\, &jj 

and then the commutator [a^Oj], 

[a i; %] 
, cxj m b m -\- /3j- m 6 m ] 

[b k , b m } + /3* fc /3* m [&|., &]„] + ajfc/5* m [&fc, &Jj + /3* fe aj m [6^,, & m ] 

&ikf3j m 3 k m ftikPtjvrfikm 
®-ik@j k Pik^jk 



'i.l 



'10 



(G.3) 











(G.4) 



so that by requiring both sets of operators to satisfy canonical commutation relations, 



we have determined two constraints (G.3) and (G.4), which in matrix notation are, 

= 1 (G.5) 



aa' 



*fi - (3* a 1 







(G.6) 



Note that in the one-dimensional case where a and (3 are complex numbers, the second 
constraint is satisfied trivially, and the remaining constraint is \a\ 2 — \(3\ 2 = 1. 

Now, we wish to apply the above Bogolyubov transformations to a system of TV 
evolving real scalar fields (f>i(x,t) described by a Lagrangian with canonical kinetic 
terms and a time dependent and non-diagonal mass matrix M 2 



S 



d A x- (cj) im (j> im -fcMfjfa) 



(G.7) 



Note that here the physical time t will be used instead of the conformal time as is used 
in the main text. We define a basis of time dependent Heisenberg operators, di(k,t) 



180 



and a\(k,t) as well as a basis of fixed time Schroedinger operators a,i(k) = 0) 
and a\(k) = a\(k,0) with initial time t — 0. The operators are indexed i = 1,...N 
with N being the number of fields. The transformation between the two bases is, 

Oi(M) = t) a j(k) + P*j(k, t)a\j(k) (G.8) 

at(M) = a *j(k,t)atj(k) + (3 ij (k,t)a j (k) (G.9) 

where now ctij(h, t) and i) are momentum and time dependent. We will suppress 
the momentum dependence and time dependence below where possible in order to 
simplify notation. The real scalar fields are written in terms of our creation and 
annihilation operators as follows, 

1 - <* -ikx 1 



<i>i 



-I 
-I 
-I 
-I 



n = 



d 3 k 
d 3 k 
d 3 k 
d 3 k 
d 3 k 



Akx 



Akx 



ikx 



e^h^aj + e-"*hUt)a] 



1 
1 



di(k, t) + e 



4(k,t) 



(G.10) 



(eeijCLj + /3*aj) + e 
(oij + Pij)aj + e~ ikx 



ikx 



l%Ji 



ikxu* 



e lkx h tJ (t) aj + e- lkx h*Jt)a 



(G.ll) 



where it is implied in the first line above that our set of Heisenberg operators di- 
agonalize the Hamiltonian at the time t. We've also defined hij(t) = a ^^-' J and 
parametrized the conjugate momentum IT via the function hij which must be deter- 
mined. The commutation relation for the fields will fix the form of II and thus ha. 



1 



d 3 kd 3 k 



[ipi(x),U m (y)} 

I eikx^k'y^^^^^^ \ 

+ e ikx - ik 'yh ij hl n [a j (k),al l (k')} 
+ e- Wfc '^*^ rnn [a](A;),o n (A; / )] 
\ + e"' y h*r h * mn [a](k),ai(k')] 

hijh m j — h^hm^j 



= iS 3 (x - y)5y 

= i5 3 {x - y)Si, 

= i5 3 (x - y)6i, 

= ?S- 
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This is a nonlinear algebraic constraint from which we can solve for h In matrix 
notation, this constraint is written, 



hh ] - h*~h T = li 

The solution may be guessed, and it is (in matrix notation), 

1 



h 
h 



"2w 



[a + 13) 
'-a + /3) 



(G.12) 

(G.13) 
(G.14) 



where to is the diagonal matrix of eigenfrequencies, and these solutions are verified on 



(G.12) by application of the constraints (G.5 G.6). Finally, the equations of motion 
will determine the evolution of the ctj/s and /%'s if we are given their initial values. 



With this goal in mind, we assume the Lagrangian introduced at the beginning (G.7) 



The mass matrix of this Lagrangian is diagonalized with a time dependent rotation 
matrix C, 

C T (t)M 2 (t)C(t) = m 2 d (t) (G.15) 

where m 2 d is diagonal. The conjugate momenta are EL = J£ = <f>, and the Hamiltonian 
written in momentum space is, 



Hamilton's Equations of motion are 



2 



0* = {Ml + k%) 



(G.16) 



n 



dC 

d<j) 



n 



-Q 2 



(G.17) 



However, Hamilton's equations are written in fields <pi and IT which do not yield 



a diagonal Hamiltonian. From our Lagrangian (G.7), we see that the fields for the 
diagonal basis will instead be, 



C 



T , 



(G.18) 
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and we thus identify <pi with the fields defined above in (G.10) which we assumed to 



be a diagonal basis. Now we may use (G.18) and (G.10 G.ll) to write Hamilton's 
equations (Heisenberg's equations of motion) in terms of the a^'s and fy's, 



a 



(-iw - I) a + (& - J)P 
P=( iuJ -I)P+(^-j)a 



(G.19) 



where matrix notation is used again, and the antisymmetric T and I matrices and the 
symmetric J matrix are defined, 



I, J 



1 



r , ± —=Ty/uj 

uj Juj 



r = c T c 



(G.20) 
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